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HETRAN - A COMPUTER PROGRAM TO SOLVE THE TWO-DIMENSIONAL 
STEADY -STATE HEAT CONDUCTION ON A CLADDED TUBE 
WITH A CONNECTING FIN 
by Frank B. Molls 
Lewis Research Center 

SUMMARY 

HETRAN, a computer program, has been written in FORTRAN IV, version 13, for 
the IBM 7094/7044 DCS to solve the two-dimensional equation V • [K(T)VT] = 0, with 
convective boundary conditions on a cladded tube with a connecting fin. The tube wall, 
cladding, and fin can be made of different materials. The required input consists of an 
initial guess to the solution, the parameters defining the region, the heat transfer coef- 
ficients on the boundaries, the coolant and gas temperatures, and the relaxation param- 
eters. Temperature dependent thermal conductivity is handled by a user supplied sub- 
routine. 


INTRODUCTION 

There are several programs available for solving the two-dimensional steady-state 
heat conduction on arbitrary regions. However, they require as input a large amount of 
detail concerning the grid that is superimposed on the r^ion. This leads to an exces- 
sive amount of manual preparation, especially if solutions are desired for a varied set 
of configurations. HETRAN is designed to overcome this difficulty for a pq.rticular type 
of geometry since it is a computer program which generates the grid. 

HETRAN was originally written to determine design criteria for rocket nozzle tubes. 
A detailed discussion of this problem may be found in reference 1. 

Following the presentation of mathematical development, the FORTRAN IV program, 
HETRAN, is presented with a complete description of the required input and the output 
obtained. The input and output for an example case is also given. 



SYMBOLS 


A 

B 

D 

D, 


g 

MAT: 


Q 

Ql 

^1 

^2 


^3 


^Ib 

^2b 

T 


T(I,j) 



T 


refg 


one of subregions of tube 
one of subregions of cladding 
one of subregions of tube 
one of subregions of fin 
one of subregions of tube 
one of subregions of cladding 

thickness of fin measured between sides of rectangle, cm 

p 

heat "transfer coefficient on coolant side of tube, W/m 

2 

heat transfer coefficient on cladded, heated gas side of tube, W/m 

conductivity of i material (i = 1, 2, 3), J/(m)(sec)(K) 

i*^ material (i = 1, 2, 3) 

one of rectangular subregions of tube 

rectangular subregion of fin 

one of subregions of tube 

one of subregions of fin 

radius to inner wall of tube (top arc), cm 

radius to interface between outer wall of tube and inner arc of cladding, cm 

radius to outer boundary of cladding, cm 

radius to inner wall of tube (bottom arc), cm 

radius to outer wall of tube (bottom arc), cm 

temperature, K 

temperature at r ,j‘’" grid point, K 
temperature of coolant, K 
temperature of gas, K 

reference temperature used as lower limit in Kirchoff Transform for mater 
ial 1, K 

reference temperature used as lower limit in Kirchoff Transform for mater 
ial 2, K 
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■^refg 

"i 

u.(I,J) 

yi 

yib 

ZL 

ZR 

^TOP 

^Z 


reference temperature used as lower limit in Kirchoff Transform for mater- 
ial 3, K 

dependent variable for i*^" material obtained from Kirchoff Transform 
dependent variable for i*"^ material at grid point 

length of bottom rectangular section of tube, cm 
length of middle rectangular section of tube and fin, cm 

perpendicular distance from top of middle rectangular section of fin to interface 
of cladding and fin, cm 

perpendicular distance from bottom of middle rectangular section of fin to 
lower end of fin, cm 

one of subregions of tube 

one of rectangular subregions of tube 

i boundary (i = 1, 2, . . . , 8) of tube, cladding, and fin 

angle corresponding to subregions A and Aj^, rad 

angle corresponding to subregions B and Bj^, rad 

angle corresponding to subregions D and D^, rad 

angle by which top circular arc of tube is less than tt/ 2 radians, rad 

angle corresponding to subregions Q and Qj^, rad 

angle corresponding to subregions ZL, rad 


STATEMENT OF THE PROBLEM 

HETRAN solves the two-dimensional equation 

V . [K(T)VT] = 0 


d) 


with boundary conditions 


-K(T) ^ = h„(T„ - T) on 3, 
0n 

-K(T) ^ = h (T^ - T) on 3 
3n » ^ 
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on 9j^, 02, 9^, 9g, 00, drj 


0n 

for the geometry shown in figure 1. The geometry represents a tube with cladding and a 
connectii^ fin. The tube is convectively cooled on its interior boundary, convectively 



heated on the outer boundary of the cladding, and insulated on the other boundaries. The 
tube wall, claddii^, and connecting fin can be made of materials with different conduc- 
tivities in which case equation (1) has discontinuous coefficients. We assume perfect 
thermal contact and continuity of heat flow at the interface of differir^ materials. 
Therefore, each point of the interface between different materials, say material 1 and 
material 2, satisfies 
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"^MAT. “ "^MAT, 

1 t 




9n 


= Kr 


MAT, 


dn 


MAT. 


METHOD OF SOLUTION 

Equation (1) and its associated boundary conditions are transformed by means of the 
Kirchoff transformation (ref. 2, sec. 2. 1) leading to Laplace’s equation on the interior 
of the region. Defining the following transformations: 



Applying these transformations to equation (1) and its boundary conditions results in a 
set of equations, one for each material region. 

For material 1: 

V^Uj = 0 (3) 


- - T) 

Bn 


on Bg 


3u. 

— = 0 on B2, 84, 9g 
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For material 2: 


= 0 

0U9 

- if “ 

9uo 

— s = 0 on a. (4) 

an 

For material 3: 

I 

V®U3 = 0 

auo 

— = 0 on 9g, drj 

an 

These equations are coupled along their interfaces by the assumptions of perfect heat 
contact and continuity of heat flow. 

Next, the composite region of figure 1 is divided into 12 subregions (fig, 2). A 
grid, part of which is radial and part rectaftgular, is superimposed on these subregions 
in such a fashion that each of the boundaries and material interfaces lies on a grid line. 

The Laplacian on the interior of each of the subregions is replaced by the regular 
five point difference formula and the resultii^ matrix solved by successive line over- 
relaxation (ref. 3, sec. 6.4). The points on the boundaries and along material inter- 
faces are solved by point overrelaxation. Since the equations at these points are non- 
linear, an iterative method for solution is required. 

The derivation of the difference equations for these special points will be illustrated 
by taking typical points on the radial grid of regions A and A^. 

Let (I, J) be an ordered pair of integers where I refers to the radial position on 
the grid and J the angular position. Then the point (1, J) is on the inner bovuidary of 
region A at the J angular subdivision. Let a cell (cross-hatched area in fig, 3) be 
taken about this point in the manner of reference 3, page 11. 

The difference equation at (1, J) is obtained in the following steps: 



/ 


around cell 


au, 

— ds 

an 
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/ 9ui f auj p 

_ids+ / — ds+ / 

X A 

a-D u-c *'c-a 


9u< 

ids 



9uw 

—ids 

9r 


^ ^ u^(l, J) - u^(l, J-1) Ar j u^(l, J) - Ug(2, J) ^ Arj 


2rj^ sin 


A«A 2 


Ar. 


-l* — 


Ui(l, J) - Ui(l, J+1) Ar. r , 

+ -i i i - h T - T(1,J) r. A0 

A0^ 2 c J i 


2rj^ sin 


(5) 


The temperature T(l, J) is not known. But assuming a previous estimate of Uj at 
(1, J), say Uj^(l, J) and a corresponding temperature T(l, J), it is possible to approxi- 
mate T(l, J). Let 


f(T) = u^ - 



— = -K.(T) 
dT ^ 

Expanding f(T) in a Taylor’s series about T(l, J) and truncating after the second 
term of the series result in 


f(T) w f(T) + (T - T) — 

dT 


T=T 


( 6 ) 


0 



Kj(z)dz - (T - T)K^(T) 
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Substituting equation (7) into equation (5), rearranging, and solving for Uj^(l,J)give 


Ari / Ar,\A0* 

u.(l, J) = 1 [u (l, J-1) + u.(l, J+D] + r. + Uj(2, J) + h 

^ A^^^- 1 ^ \ ^ 2 / Arj ^ 


c^l 


4rj sin 



where w is the relaxation parameter. The temperature T(1,J) corresponding to 
Uj(l, J) is obtained from equation (7). 

The difference equations for a material interface point are obtained in a similar 
fashion. The cells about the point (N, J) of region A and the correspondii^ point (1,J) 
of region Aj^ have a common boundary (path b — c of fig. 4). The cell for (N, J), say 



Figure 4. - Detail of grid for material interface point regions A and Aj, 
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cell^, is cross-hatched area in figure 4, and the cell for (1, J), say cellg, is the stippled 
area in figure 4. 

Again using the divergence theorem results in 


0 = 


/ 


cell. 


V^uj dA + 


/. 


cellc 


V^ug dA 


<£ !!i 

•ground cellj^ 


/ 

around cellg 


9Uft 

— ds 
9n 


0 = 



3u, 

ds + 

0r 





ds 


+ 






( 8 ) 


The assumption of continuity of heat flow implies 



0U, 0U« 

— i ds + / — = ds = 0 

0r / 0r 

*<;-b 


Dropping these integrals out of equation (8) and continuity with the approximation lead to 
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0 


Uj(N, J) - Uj(N, J-1) Arj Uj(N, J) - Uj(N, J+1) Ar^ u^(N, J) - u^(N-l, J) 

■■■■ ' ■ "I" “ •4' "■ ■ I 


A0. 2 A0* 2 

2rg sin 2rg sin — — 


Ar, 


.4 U2(l,J)-U2(l,J-l)Ar2 U2(1,J)-U2(2,J)/ Ar2\ 

2 o 


UgCl, J) - Ugll, J+1) Arg 


2r£ sin 


A0a 2 




Ar.\ 


+ r. 


i A0, 


2 x 2 sin 


A0. 


2/ 


Ui(N,J) 


^^2 / ^^2\ 


2rg sin 


U2(1,J) 


4rg sin 


Ar, / ArA A0* 

i— [uj(N, J-1) + u^(N, J+1)] + ^2 - ^ ^i(N-l, J) 


Arr 


4ro sin. 


A0, 


[ugCl, J-1) + U2(l, J+1)] + 1^2 + Ug(2, J) 


(9) 


Letting TI(J) be the temperature associated with U 2 (l, J) and T(J) the temperature 
associated with Uj^(N, J) and using the same approximation as in equation (7) give 


/'T(J) 

J Kj(z)dz 


u.|(N, J) refi _ 

T(J) w ■ V ^ pi — , + T(J) 

Kj[t(J)] Kj^Lt(J). 


( 10 ) 
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TI(J) w 


+ TI(J) 


( 11 ) 


U2(1,J) 

K2[tI(J)] K2[tI(jJ 


However, the assvimption of perfect thermal contact at a material interface implies that 
T(J) = TI(J). Subtractii^ equation (10) from equation (11) results in 


/•T(J) 

/ Ko(z)dz / 

" •/iTi 

refo refi 


T(J) 


0 = 


Uj(N,J) U2(1,J) ^ref2 


K^(z)dz 


k£t(J)] K2[t(J)] K2[t(J)] 


Ki[t(j)] 


/•T(J) /-T(J 

4 L 

■"■refo ^ref^ 


T(J) 


Uj(N,J)^ U2(1,J) 


K^(z)dz 


Kj[t(J)] K^^T(J)] K2[t(J)] 


Kj[t(J)] 


( 12 ) 


Equations (9) and (12) are solved simultaneously for u^(N,J) and U2(l, J). A new 
estimate for T(J) is then obtained by adding equations (10) and (11) to give 


T(J) = T(J) + 



r /'t(j) 

/'T(J) , , 1 


/ K,(z)dz 

I K^(z)dz 

U2(1,J) Ui(N,J) 


"^refi 

- - - V 

I 1 ““ 1| 

K2|t(J)] Ki|t(J)] 

K2[t(J)] 

Ki|t(j3 



J 


2 


Convergence of the iteration process is assumed to occur when 


heat in - heat out 
heat in 


< tol 


(13) 


where tol is some preassumed number and 




Trapezoidal integration is used to approximate these integrals so as to be consistent 
with the difference equations. 

When convergence has occurred, the temperatures corresponding to the u. func- 
tions are determined by iterating on equations of the same type as equation (7). 

The program was converted to double precision after it was found that the large as- 
pect ratios of the cells of the grid led to an appreciable loss of significance. The pro- 
gram takes less than 1 minute to do 175 iterations for 1030 grid points. 


PROGRAM PROCEDURE 


The main program is HETRAN. There are 17 subroutines: CONSA, CONSB, 

RINIT, COMPA, COMPB, RCOL, CIRCOL, CIRCOT, RLINE, CLINE, OUTPUT, CONI, 
CON2, CONS, SUBTl, SUBT2, and SUBT3. The subroutines CONI, CON2, CONS, 
SUBTl, SUBT2, SUBTS, and part of RINIT must be supplied by the user. For more in- 
formation on these subroutines, see the section labeled USER SUPPLIED SUBROUTINES. 
The calling relation of all subroutines is shown in figure 5. 




SUBT3 

SUBT2 

SUBTl 



C01\13 

CONI 

CLINE 

RUNE 



C0N3 

C0N2 

CONI 

CLINE 

RUNE 


CIRCOL 

RCOL 



CIRCOT 

CIRCOL 

RCOL 


Figure 5. - Calling relation of subroutines. 
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HETRAN reads and prints all input data and controls the calls to the subroutines. 

Subroutines CONSA calculates the grid constants for subregions A, Al, D, Dl, B, 
and Bl and the point relaxation coefficients for the interfaces and boundaries of these 
regions. Through calls to RCOL, CIRCOL, and CIRCOT, it obtains the line relaxation 
coefficients for the aforementioned regions. 

Subroutine CONSB calculates the grid constants for subregions P, Q, Q^, ZR, 
and ZL and the point relaxation coefficients for the interfaces and boundaries of these 
regions. Through calls to RCOL and CIRCOL it obtains the line relaxation coefficients 
for the aforementioned regions. 

Subroutine COMPA calculates the u.'s at each point of the grid for subregions A, 
Aj, D, Dj, B, and Bj^ and the temperatures at interface and boundary points of these 
regions. 

Subroutine COMPB calculates the u.’s at each point of the grid for subregions P, 
Pjj Q> Qjj ZL, and ZR and the temperatures at interface andboimdary points of these 
regions. 

Subroutine RCOL factors a tridiagonal matrix with constant off-diagonal elements 
into a product of an upper and lower triangular matrix. 

Subroutine CIRCOL computes the five point difference coefficients for Laplace's 
equation on an annulus where the increments in the radial direction are constant. Next, 
it factors the tridiagonal matrix of coefficients (coupling the variables along a radial) 
into a product of upper and lower triangular matrices. 

Subroutine CIRCOT computes the five point difference coefficients for Laplace's 
equation on an annulus where the increments in the radial direction are variable. Next 
it factors the tridiagonal matrix of coefficients (coupling the variables aloi^ a radial) 
into a product of upper and lower triai^ular matrices. 

Subroutine OUTPUT prints the points of the machine generated grid and the temper- 
atures at each point of the grid* The temperatures are obtained through calls to SUBTl, 
SUBT2, and SUBT3. In addition, if IPUNCH is set equal to one in the input data, OUT- 
PUT will punch 414 BCD cards which contain the last computed u^’s at each point of 
the grid and temperatures at each point of the interfaces and boundaries. The number 
of cards punched is a fixed quantity that can only be charged by chai^ii^ the BCD dump 
statements in the section OUTPUT. 


INPUT 

Figure 6 shows the input variables as they are pxmched on data cards. The first 
input card is for a title which serves as problem identification. The remaining cards 
are for input variables. Fixed point variables must be punched in a five column field 
and floating point variables (decimal point must be pmiched) in an eight column field. In 
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Figure 6. - Input form. 


addition to the data cards, the user must supply a subroutine for an initial guess. For 
an explanation see the section labeled USER SUPPLIED SUBROUTINES. 


Input Variables and Definitions 

NRA number of division points for radius, T 2 - (see fig. 2), <15 

NTA number of division points for angle (see fig. 2), <24 

NRB number of division points for length fj^ (see fig. 2), 4 < NRB < 5 

NTB number of division points for angle (see fig. 2), 

NRB + NTB + NTD - 2 < 35 

NTD number of division points for angle Ojy (see fig. 2) 

NXP mrniber of division points for middle rectangular length x^^ (see fig. 2), 

<20 

NTQ number of division points for angle 0 q (see fig. 2), <20 

NZ number of division points for angle 9^ (see fig. 2), <20 

NZX number of division points for lower rectangular length Xj^ (see fig. 2), <19 

IPUNCH trigger to OUTPUT subroutine: 

= 0 No cards punched 

= 1 414 BCD cards punched containing u’s at each point of grid and 

temperatures along boundaries and interfaces. The number of 
cards punched is fixed by the dimension statements in the 
program 
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IRE AD 


trigger to RINIT subroutine: 

= 0 User must set initial guess for u’s 

= 1 RINIT reads 414 BCD cards containing u's and temperatures 
aloi^ boundaries and interfaces 

NCOUNT maximum number of iterations program will be cycled. After NCOUNT 
iterations, output is printed and next case is processed. 

MCOUNT after MCOUNT iteration, heat in, heat out, and relative error are printed 
(see eqs. (13) and (14)) 

ICOS trigger to activate a cosine variation for the heat transfer coefficient on 

the claddii^: 

= 0 EGAS = EGAS 

= 1 EGAS = EGAS COS 0 

Rl radius to inner tube wall (top arc) (see fig. 2) 

R2 radius to interface of outer tube wall and inner radius of cladding (see 

fig. 2) 

R3 radius to outer boundary of cladding (see fig. 2) 

YFINT perpendicular distance from top of middle rectangle to interface of cladding 
and fin (see y^ in fig. 2) 

FINTEK fin thickness measured between sides of rectangle (see fj^ in fig. 2) 

ANGTOP angle in radians by which top circular arc of tube is less than 90° (see 

®TOP 

ECOOL heat transfer coefficient on inner radius of tube 

TCOOL temperature of environment of inner radius of tube 

EGAS heat transfer coefficient on outer boimdary of cladding 

TGAS temperature of enviroiunent on outer boundary of cladding 

XLONGM length of middle rectangular section of tube (see x^ in fig. 2) 

R2B radius to outer wall of tube (lower arc) (see rgb in fig. 2) 

TOL trigger to stop iterations. If heat balance (eq. (13)) :< TOL, iterations are 

stopped and output is printed. 

WMEGA relaxation factor for subregion A (see fig. 2) 

WMEGAl relaxation factor for subregion A1 (see fig. 2) 

WMEGD relaxation factor for subregion D (see fig. 2) 
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WMEGDl 

relaxation 

WMEGB 

relaxation 

WMEGBl 

relaxation 

WMEGP 

relaxation 

WMEGPl 

relaxation 

WMEGQ 

relaxation 

WMEGQl 

relaxation 


factor for subregion D1 (see fig. 2) 

factor for subregion B (see fig. 2) 

factor for subregion B1 (see fig. 2) 

factor for subregion P (see fig. 2) 

factor for subregion PI (see fig. 2) 

factors for subregions Q, ZL, and ZR (see fig. 

factors for subregion Q1 (see fig. 2) 


2) 


Instructions for Preparing Input 

Units of measurement. - Any compatible system of units may be used for linear 
and thermal input data. In the example problem the English system was used with 
lengths given in inches and temperatures in degrees Rankine. The one ai^le that is in- 
put, ANGTOP, must be in radians. 

Relaxation factors. - The optimum relaxation factors are dependent upon the geom- 
etry and material. They will generally have to be fovind by experimentation. Provisions 
have been made for different relaxation factors in the various subregions. 

Case of no cladding . - It is possible, by inputii^ R3 = R2, to rvm a case without 
cladding on the tube. The outer tube wall and the top boundary of the fin are then con- 
sidered to be the heated surfaces. 

Format for input data. - All fixed point variables are punched in a right-adjusted 
five column field; all floating point variables (punch decimal point) are in an eight 
column field. 


USER SUPPLIED SUBROUTINES 

The user must supply the following subroutines: CONI, CON2, CON3, SUBTl, 
SUBT2, and SUBT3. If an initial guess is not available on BCD cards, then an initial 
guess to the u.'s at each point of the grid and the corresponding temperatures at inter- 
face and boTmdary points must be supplied in RINIT. 

Subroutines CONI, CON2, and CON3 have the same variable list. The subroutines 
perform the same function, but for different materials; CONI refers to material 1, 
CON2 to material 2, and CON3 to material 3. Consequently, a description of CONI is 
sufficient for all tiiree routines. 
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The purpose of CONI is to retum to the calling program the conductivity of mater- 
ial 1 at a given temperature and the variable 



Kj(x)dx 


Kj(T) 


Usage; Subroutine CONl(TOLD, U, NTIM, COND) where 


TOLD 

U 


vector of temperatures supplied by calling prc^ram 

vector ^ Kj^(x)dy^j(T) computed in CONI 
^ref- / 


NTIM number of elements in vector TOLD supplied by calling prc^ram 

COND vector of conductivities for material 1, computed in CONI, corresponding to 
vector of TOLD’s 


TOLD, U, and COND must be declared double precision and dimensional variables of at 
least dimension one in CONI. 

Method; The subroutine must have a ’*DO IX)OP*’ to calculate from 1=1, NTIM 
the conductivities, COND(I), for material 1, correspondii^ to the given temperatures 
TOLD(I). It must also calculate the variables 


PT 

cAref. 

U(I) = 3 


Ki(x)dx 


TOLD(I) 


The subroutines SUBTl, SUBT2, and SUBT3 must return the temperatures corre- 
sponding to a given vector of u.'s for materials 1, 2, and 3, respectively. Since the 
variable list of these three subroutines is identical and their function is the same, a 
description of SUBTl is sufficient to describe all three subroutines. 

The purpose of SUBTl is to calculate a vector of temperatures for material 1, 
given the correspondiiig u.'s. The method which is useful for this calculation, is the 
Newton method of equation (7). 

Usage; Subroutine SUBT1(U, Tl, TEMP, Nl, N2) where 

U vector of Uj's, supplied by calling prc^ram 
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Tl initial guess to TEMP(Nl) supplied by calling program 

TEMP vector of temperatures, calculated by SUBTl, corresponding to vector of u's 
Nl initial point of ”DO LOOP” supplied by calling program 
N2 final point of ”DO LOOP” supplied by calling program 

TJ and TEMP must be declared dimensioned variables in SUBTl with dimension of at 
least one. 

Method; The subroutine must have a ”DO LOOP” to calculate from I = Nl, N2 
the temperatures, TEMP(I), for material 1, corresponding to the given U(I). Newton’s 
method as given by equation (7) is an effective method for performing this calculation. 

Subroutine RINIT sets the initial guess for the u. at each grid point and the initial 
guess for the temperatures at each grid point of the boundaries and interfaces. These 
temperatures must correspond to the u. at the given points. 

There are two ways to set the initial guess in RINIT. If IREAD is set equal to 1 on 
the input data, then RINIT will attempt to read 414 BCD cards. These cards are gener- 
ated by subroutine OUTPUT and contain the complete field of u. and the temperatures 
along the boundaries and interfaces from the last iteration. If IREAD is set equal to 0 
as the input data, then the user must supply the appropriate initial guess. 

A constant temperature field can be inserted as a first guess by changing four cards 
in RINIT. First, card 22 of the listing for RINIT must be changed to the given temper- 
ature. Then, cards 26, 28, and 40 must be charged to the Uj^, Ug, and Ug correspond- 
ing to the temperature on card 22. For any other type of initial guess, cards 21 to 66 
must be deleted and the appropriate guess inserted. 


OUTPUT 

Sample output is given in the appendix together with sample input. The units for 
linear and temperature output are the same as for input. The sample input used inches 
and degrees Rankine; the output \mits are also in these units. All ai^lar output is 
given in degrees. 

In addition to the sample output given in the appendix, 414 BCD cards can be ob- 
tained by setting IPUNCH equal to one. This output is suitable for restarting the pro- 
gram from the last iteration, or as an initial guess for a similar case in which the geo- 
metric parameters have been varied. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 21, 1971, 

132-80. 
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APPENDIX - COMPUTER LISTING WITH SAMPLE INPUT AND OUTPUT 


Sample Input 


S&MPLE CASE — STAINLESS TUBE, NICKEL CLAD D ING, COPPER FIN— NO COS VAR* 


NRA NTA NRB NTS NTD NXP 

9 20 4 5 5 10 


NTQ NZ NZX IPUNCH IREAO NCCUNT MCQUNT ICOS 

10 10 10 C 0 3000 2 50 -0 


R1 

0.40500CCOE-01 


R2 

R3 

Y-FIN-TOP 

FIN- THICK 

ANG L6-T0P 

0.52500000E-01 

0.725)00006-01 

0.12600000E-01 

0.95959599E-03 

0 


H-COOLANT 

0*14805000F-‘OI 


TEMP-COOLANT 
70. 0000000 


H-GAS 

0.1 942 9000 E-02 


TEM P-GAS 
5897.00000 


X-MID-RECT 

0.325nCQCOE-01 


R2-B0T 

0.40000000E-01 


Y-F IN- BOTTOM 
0.94999999E-02 


TOLERANCE 

0.999999996-04 


DMEGA(A) 

1.91999999 


OMEGA(Al) 

1. 91999999 


□MEGAIO) 

1.91999999 


OMEGA! Dl) 

1.91999999 


OMEGA(B) 

1.91999999 


OMEGA <B1) 

1.91999999 


OMEGA(P) 

1.91999999 


QMEGA(Pl) 

1.91999999 


OMEGA! QJ 

1.91999999 


OMEGA! Ql) 

1.91999999 


Sample Output 


ITERATION 

ITERATION 

ITERATION 

ITERATION 

ITERATION 

ITERATION 


2 50 
500 
750 
1000 
12 50 
1421 


HEAT IN = 
HEAT IN = 
HEAT IN = 
HEAT IN = 
HEAT IN = 
HEAT IN = 


0.52211515 
0. 52490184 
0.52544611 
0. 52555250 
0.52557338 
0.52557682 


HEAT OUT = 
HEAT OUT = 
HEAT OUT * 
HEAT OUT = 
HEAT OUT. = 
HEAT OUT * 


-0*60466251 

-0.54039369 

-0.52846774 

-0.52615677 

-0.52570375 

-0.52562915 


RELATIVE ERR « 0.15810259 

RELATIVE ERR = 0. 29513791E-01 

RELATIVE ERR * 0. 57506089E-02 

RELATIVE ERR = 0. 11497 699 E-02 

RELATIVE ERR * 0. 24812427E-03 

RELATIVE ERR * 0. 995 5 8095 E -04 


TEMPERATURES ON THE TOP RADIAL SECTION OF TUBE 


THBT - 0. 


R 


T 


R 


T 


R 


R 


T 


0.40500CE-01 

0.4650COE-C1 

0.525000E-01 

O. 525 OOOE-OI 

O.535096E-O1 

0.571359E-01 


0.770166E^‘C3 
0.102212E^-04 
0.122742E+04 
0.122742E+04 
0.124270E+04 
0.12 9236E 4-04 


0.420000E-01 

0.480000E-01 

0.528333E-01 

0.538598E-31 

0.605076E-01 


0.838 638E403 
0.107715E+04 

0.123 210 E+04 
0.124649 e+04 
0.133917E+04 


0.435000E-01 

0.495000E-01 

0.531667E-01 
0.543152E-01 
0.6 54288E-01 


0.903183E4-03 
0. 112956b 4-04 

0.123677E+04 
0.125288E *04 
C.140611E4-04 


0.450000E-01 

0.510000E-01 

0.535000E-01 
0.549637E-01 
0. 7250006-01 


0.964230E4-03 

0.117954E4-04 

0.124144E4-04 

0.126197E4-04 

0.149903E4-04 
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THET 


0.250292E+01 


TRTRTRT 

0.4050n0E“01 0.769639E+03 0.420000E-01 0.838079E<^3 0.435000E-01 0.902596E+03 0.450000E-01 0«9636l8E*03 

0.465000E-01 0.X02149E+04 0.480000E-01 0.107650E+04 0.495000E-01 0.112889E+04 0.510000E-01 0.117886E+04 

0.525000E-01 0.122672E+04 

O. 525000E-P1 0.122672E+-04 0.528333E-01 0.123139E*04 0.5316676-01 C.123607E+04 0.535000E-01 0.124074E+04 

P. 535896E-01 0. 124199E 4-04 0.538598E-01 0 . 124578 E4-04 0.543152E-01 0.125216E+04 0.549637E-01 0.126125E+04 

O. 571359E-Q1 C.129163E+04 0.605076E-01 0.133843E404 0.654288E-C1 C.140537E+04 0.725000E-01 0.149832E+04 

THET = 0.5005 84E4-01 

R T RTRTRT 

P. 405000E-01 0.768057E+03 0.420000E-01 0.836401E+03 0.435000E-01 0.9C0833E+03 0.450000E-01 0.961780E+03 

0.465C0CE-OI 0.101959E4-04 0.480000E-01 0.107454E+04 0.495000E-01 0.112688E+04 0.5100006-01 0.117681E+04 

C.5250C0E-01 0.122463E4-04 

0.525000E-01 0.122463E4-04 0.528333E-01 0. 122929 E*04 0.531667E-CI 0.123395E+04 0.535000E-01 0.123861E^04 

0.535896E-01 P. 1 23 98 7E 4-04 0.538598E-01 0.124365E+04 0.543152E-01 0.1250C2E+04 0.549637E-01 0.125909E+04 

0.571359E-01 0. 128944E4-04 0.605076E-01 0.133621E+04 0.654288E-C1 C.140317E4-04 0.725000E-01 0.149619E+04 

THET >= 0.750876E4-01 

R TRTRTRT 

0.405t>OOE-01 0. 76541 8E4-03 0.420000E-01 0.833601E4-03 0.435000E-01 0.897890E+03 0.450000E-01 0.958710E4-03 

0.465000E-01 0.101641E+04 0. 480D00E-01 3.107 L26E+Q4 0.495000E-01 0.112352E+04 0.510000E-01 0.117339E^04 

C.525000e-Pl 0-122112E+04 

0.525P00E-01 0.122112E + 04 0.528333E-01 0.122577E+04 0.531667E-01 C.123041E+04 0.535000E-01 0.123506E4-04 

0.535896E-01 0. 123631£4^04 0.538598E-01 0.124008E+04 0.543152E-01 0.124644E+04 0.549637E-01 0.125549E+04 

0.571359E-01 0.128577E4-04 0.605076E-01 0.133250E+04 0.654288E-C1 C.139948E+04 0.725000E-01 0.149262E+04 

THET = 0.100117E+02 

R TRTRTRT 

0.405000E-01 0. 761 716E 4-03 0.420000E-31 0.829673E+03 0.435000E-01 C.893760E+03 0.450000E-01 0.954401E4-03 

0.46500CE-01 0.101194E+04 0.480000E-01 0.1D6666E+04 0.495000E-01 0.111880E+04 0.510000E-01 0.116857E+04 

0.525000E-01 0.121 61 9E4^04 

0.525000E-01 0.121 619E4-04 0.528333E-01 0.122082E4.Q4 0.531667E-01 0.122544E+04 0.535000E-01 0.123007E+04 

0.5358966-01 0.12313264-04 0.538598E-01 0.123507E4-04 0.543152E-01 0.124140E+04 0.549637E-01 0.125042E+04 

O. 571359E-P1 0.128061E+04 0.605076E-01 0.132728E4^04 0.654288E-01 0.13S428E+04 0.7250006-01 0.148760E4-04 

THET = 0.125146E+02 

R TRTRTRT 

P, 405000E-C1 0.756947E+03 0.420000E-01 0. 824610 E+03 0.435000E-01 C.888^34E4-03 0.4500006-01 0.948841E4-03 

0.465000E-01 0.100617E+P4 0.<t80000E-31 0 .10607164-04 0.495000E-01 0-111270E+04 0.510000E-01 0.1 16234E4-04 

O. 525000E-01 0.12098364^04 

P. 525000E-01 0. 1209836+04 0.528333E-01 0.121442E+04 0.531667E-01 0.121902E+04 0.535000E-01 0.122362E*04 

0.535896E-01 0.122486E+04 0.538598E-01 0.122859E+04 0 . 543l52E-jGl 0.123^896+04 0.549637E-01 0. 1243866+04 

0.571359E-01 0.1273936+04 0.605076E-01 0.132051E+04 0.654288E-01 0.138754E+04 0.725000E-01 O.X48i08E+04 

THET = 0.150175E+02 

R TRTRTRT 

0.4050006-01 0.7511046+03 0.420000E-01 0.818404E+03 0.435000E-01 0.8819C2E+03 0.450000E-01 0.942017E+03 

0.465000E-01 0. 9990906+03 0.480000E-01 0.105340E+04 0.495000E-01 0. 11051 9E+04 0.5100006-01 0.115466.E+04 

0.525000E-01 0.120199E+04 

0.525000E-01 0. 1201996+04 0.528333E-01 0.120655E+04 0.531667E-01 0.121111E+04 0.535000E-01 0.121568E+04 

0.535896E-01 0.121691E+04 0.5385986-01 0.122061E+04 0.543152E-01 0.122686E+04 0.5496S7E-01 0.123578E+04 

0.571359E-01 0.126570E+04 0.605076E-D1 0.131214E+04 0.654288E-01 0.137921E+04 0.725000E-01 0.147303E+04 

THET = 0.175204E+02 

RTRTRTRT 

0.405000E-01 0.744179E+03 0.420000E-01 0.8H046E+03 0.435000E-01 0.874152E+03 0.450000E-01 0.933915E+03 

0.465000E-01 0.990675E+03 0.483000E-01 0.104471E+04 0.495C00E-01 C.109626E+04 0.510000E-01 0.1145526+04 

O.525C0OE-01 0.H9266E+O4 

O.525000E-C1 0.119266E+04 0.528333E-01 0.119717E+04 0.531667E-01 0.120169E+04 0.535000E-01 0.120622E+04 

0.535896E-01 0.1207456+04 0.538598E-01 0.121111E+04 0.543152E-01 0.1217306+04 0.549637E-01 0.122615E+04 

0.571359E-01 0.1255876+0^ 0.605076E-01 0.130215E+04 0.654288E-01 0.136924E+04 0.7250006-01 0.146338E+04 

THET = 0.200234E+02 

R TRTRTRT 

0.4050006-01 0.736163E+03 0.420090E-01 0.802527E+03 0.435000E-01 C.865172E+03 0.450QOOE-01 0.924521E+03 

0.4650COE-01 0.9809086+03 0. 480000E-01 0.103461E+04 0.495000E-01 0.108587E+04 0.510000E-01 0. 1134876+04 

0.525000E-01 0.1181806+04 

O.525C0OE-C1 0.1181806+04 0.528333E-01 0.118625E+04 0.5316676-01 0.119072E+04 0.535000E-01 0.119519E+04 

0.535896E-01 C.119639E+04 0.538598E-01 0.1230036+04 0.543152E-01 0.120616E+04 0. 549637E-01. 0.1214926+04 

0.571359E-01 0.124440E+04 0.605076E-01 0.129045E+04 0.654288E-01 0.135755E+04 0.725000E-01 0.145203E+04 
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TH€T 


0.2252636+02 


R TRTRTRT 

O.A050006-01 0.7270646+03 0,4200006-01 0.7928396+03 0.4350006-01 0.8549526+03 0.4500006-01 0.9138186+03 

0.4650006-01 0.9697716+03 0.4800106-01 0. 102 3D9 6+04 0.4950006-01 0.1C7399E+04 0. 5100006-01 0.112269E+-04 

O. 525CC06-C1 0,1169356+04 

P. 5250006-01 0. 1169356+04 0.5283336-01 0.1173756+04 0.5316676-01 0. 1178156+04 0.5350006-01 0.1182566+04 

0.5358966-01 0.1183756+04 0.5385986-11 0.1187336+04 0. 5431526-01 0.1193396+04 0.5496376-01 0-1202046+04 

0.5713596-01 0,1231236+04 0.6050766-01 0,1277016+04 0. 6542886-01 C. 1344076+04 0.725000E-01 0.1439046+04 

THBT = 0.2502926+02 

RTRTRTRT 

0,4050006-01 n. 7168636+03 0.420000E-D1 0.7819756+03 0. 4350006-01 C. 8434826+03 0.4500006-01 0.9017946+03 

0.465000E-C1 0.957244E+03 0.4800006-01 0. 1010116+04 0. 4950006-01 0-106C61E+04 0.5100006-01 0.1108956+04 

0.525000E-C1 0.1155296+04 

0.525rOOE-Ol 0,1155296+04 0. 5283336-01 0.1159616+04 0,5316676-01 C. 1163946+04 0.5350006-01 0.1168286+04 

0.5358966-01 0,1169456+04 0,5385986-01 0.1172986+04 0,5431526-01 C. 1178946+04 0.5496376-01 0.1187476+04 

0,5713596-01 0.121631E+04 0.6050766-01 0.1261736+04 0.6542886-01 0.1328706+04 0.7250006-01 0. 1424166+04 

THET = 0,2763216+02 

R TRTRTRT 

0,4050006-01 0.7055656+03 0,4200006-01 0.7699326+03 0.4350006-01 C. 8307526+03 0.45000CE-01 0.8884356+03 

O. 465000E-C1 0.9433105+03 0.4800006-01 0.9956496+03 0.4950006-01 0.1045686+04 0.5100G0E-01 0.1093596+04 

P. 5250006-01 0.1139566+04 

C. 5250006-01 0.113956E+04 0.5283336-01 0,1143796+04 0.5316676-01 C.1148C4E+04 0.5350006-01 0,1152306+04 

0.5358966-01 P.115344E+04 0,5385986-01 0.1156916+04 0,5431526-01 0.1162776+04 0.5496376-01 0.1171166+04 

C. 5713596-01 0.1199586+04 0.6050766-01 0.1244566+04 0.654288E-C1 C. 1311346+04 0. 7250006-01 0.1407336+04 

fHET = 0.3003506+02 

R TRTRTRT 

0.4050006-01 0.6931716+03 0.4200006-01 0.756 7086+03 0.4350006-01 0.8167596+03 0.4500006-01 0.8737326+03 

0.4650006-01 0.927955E+03 0.4800006-01 0.9796966+03 0.4950006-01 0.1029186+04 0.5100006-01 0.1076606+04 

0.5250006-01 0.1122136+04 

G. 5250006-01 0.1122136+04 0.5283336-01 0.1126266+04 0.5316676-01 0,1130416+04 0.535000E-01 0.1134586+04 

0. 5358966-01 0.1135706+04 0.5385986-01 0.1139096+04 0.5431526-01 0. 1144836+04 0.5496376-01 0.H5305E+04 

0.5713596-01 0.1180996+04 0.6050766-01 0.1225416+04 0.654288E-C1 0.1291696+04 0.7250006-01 0.1388426+04 

THET = 0.3253806+02 

R TRTRTRT 

0.4050006-01 0.679689E+03 0.4200006-01 0.7423076+03 0.4350006-01 C. 8015026+03 0.4500006-01 0.8576816+03 

0.4650006-01 0.9111676+03 0.4800006-01 0.9622306+03 0.4950006-01 0.1011C9E+04 0.5100006-01 0.1057946+04 

0.5257006-01 0.1102956+04 

0. 5250006-01 0.1102956+04 0.5283336-01 0.1106976+04 0.5316676-01 C.llllClE+04 0.5350006-01 0.1115076+04 

0. 5358966-01 0.1116166+04 0.5385986-01 0.1119476+04 0.5431526-01 0.1125076+04 0.5496376-01 0. 1133116+04 

0.5713596-01 0. 1160476+04 0.6050766-01 0.1204226+04 0.6542886-01 0.1270236+04 0.7250006-01 0.1367246+04 

THET = 0.350 409 6+02 

R TRTRTRT 

0. 4050006-PI 0.6651316+03 0.4200006-01 0.7267386+03 0.4350006-01 0.7849886+03 C.450000E-01 0.8402826+03 

0.4650006-01 0.8929455+03 0.480000E-OI 0.9432406+03 0.4950006-01 0.9913906+03 0.5100006-01 0.1037586+04 

0. 5250006-01 0. 1081986+04 

0.5250006-Cl 0.1081986+04 O.528333E-01 O.1O8580E+O4 0.5316676-01 0. 1089806+04 0.5350006-01 0.1093746+04 

0.5358966-01 0.1094806+04 0.5385986-01 0.1098026+04 0.5431526-01 0. 1103466+04 0.5496376-01 0.1111296+04 

0. 5713596-01 0.1137996+04 0.5050766-01 0.1180926+04 0.6542886-01 0.1246256+04 0.7250006-01 0.1343606+04 

TH6T * 0.3754386+02 

R TRTRTRT 

0.405000E-01 0.6495156+03 0.4200006-01 0.710017E+03 0.4350006-01 0. 7672286+03 0.4500006-01 0.8215456+03 

0.4650006-01 0.8732906+03 0.4800006-01 0.9227256+03 0. 4950006-01 C. 9700716+03 0.5100006-01 0.1015516+04 

0.5250006-01 0.1059226+04 

0.5250006-01 0.1059226+04 0.5283336-01 0.1052986+04 0.5316676-01 C. 1066766+04 0.5350006-01 0. 1070576+04 

0.5358966-01 0.1071606+04 0.538598E-01 0.1074706+04 0. 5431526-01 0.1079986+04 0.5496376-01 0.1087566+04 

0. 5713596-01 0.1113506+04 0.6050766-01 0.1155456+04 0.6542886-01 0.121984E+04 0. 725000E-01 0.1317226+04 

THET = 0.4004676+02 

R TRTRTRT 

0.4050006-01 0.6328676+03 0.4200006-01 0.692 1686 +03 0.4350006-01 G. 7462436+03 0.4500006-01 0.80 I486 6+03 

0.4650006-01 0.8522166+03 0.4800006-01 0.9006936+03 0.4950006-01 0.9471356+03 0.5100006-01 0.9917286+03 

0. 5250006-01 0.1034636+04 

0.5250006-01 0.1034636+04 0.5 283336-31 0.103824E+04 0.5316676-01 0.1041886+04 0.5350006-01 0.104554E+04 

0. 5358966-01 0.1046536+04 0.5385986-01 0.1049526+04 0.543152E-01 0.105460E+04 0.5496376-01 0.106 191E+04 

0.5713596-01 0.1087006+04 0.605076E-01 0. ’*'‘7796+04 0.654288E-01 C.119091E+04 0.7250006-01 0. 1287816+04 
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THET = 0.425496E^02 


R 

T 

R 

T 

R 

T 

R 

T 

0 ♦4050006-01 

0.6152206^-03 

0.420000E-01 

0.6732226+03 

0.4350006-01 

0. 7280626+03 

0.4500006-01 

0.7801326+03 

O.465QOOE-01 

C.525000E-01 

0.8297476+03 

0.1008236+04 

0.480000E-01 

0.8771646+03 

0.495000E-C1 

0.9225996+03 

0.5100006-01 

0.966235E+03 

0.5250006-01 

0.100 8236 +04 

0.5283336-01 

0.101168E+04 

0.5316676-01 

Q. 1015166+04 

0.5350006-01 

0.1018666+04 

0.5358966-01 

0.1019616+04 

0.5385986-01 

0.1022476+04 

0.5431526-01 

0,1027346 +04 

0.5496376-01 

0. 1034356+04 

0.5713596-01 

0.1058496+04 

0. 6050766-01 

0.1097936+04 0,6542886-01 

THET = 0.4505256+02 

0. 1155386+04 

O.725QQOE-0L 

0.1255006+04 

R 

T 

R 

T 

R 

T 

R 

T 

0.4050006-01 

0. 5966156+03 

0.4200006-01 

0.653219E+03 

0.4350006-01 

0.706724E+03 

0.4500006-01 

0.7575206+03 

0.4650006-01 

0.5250006-01 

0.805917E+03 

0.9800356+03 

0.483000E-01 

0.8521716+03 

O.495CO0E-O1 

0. 8964936+03 

0.5100006-01 

0.9390636+03 

0.5250006-01 

0.9800356+03 

0.528333E-01 

0.9833136+03 

0.5316676-01 

C. 9866196+03 

0.5350006-01 

0.9899546+03 

0.5358966-01 

0.9908546+03 

0.538598E-01 

0.9935836+03 

0.5431526-01 

0. 9982226+03 

0,5496376-01 

0.1004926+04 

0.5713596-01 

0.1028026+04 

0,6050766-01 

0.106 592 6+04 0.65428 86-01 

THET = 0.4755556+02 

C.112525E+04 

0,7250006-01 

0.1218206+04 

R 

T 

R 

T 

R 

T 

R 

T 

0.4050006-01 

0. 5771046+03 

0.420000 E-31 

0.632 207 6+03 

0,4350006-01 

0.6842776+03 

0.4500006-01 

0.7336966+03 

0.4650006-01 

0.5250006-01 

0.780773E+03 

0.9500896+03 

0,4800006-01 

0.8257596+03 

0.4950006-01 

0.8688636+03 

0.5100006-01 

0.91025SE+03 

0.5250006-01 

0. 9500896+03 

0.5283336-01 

0.9531846+03 

0,5316676-01 

0.9563096+03 

0.5350006-01 

0.9594646+03 

0.5358966-01 

0.9603166+03 

0.5385986-01 

0.9629016+03 

0.5431526-01 

0.9672996+03 

0,5496376-01 

0. 9736546+03 

0.5713596-01 

0.9956516+03 

0.6050766-01 

0.103189 6+04 0.6 542886-01 

THET = 0.5485356+02 

0.1088706+04 

0.7250006-01 

0. 1176506+04 

R 

T 

R 

T 

R 

T 

R 

T 

0.4050006-01 

0.5152886 +03 

0.4200006-01 

0.5654246+03 

0.4350006-01 

0.6127066+03 

0.4500006-01 

0. 6575056+03 

0. 4650006-01 
0.5250006-01 

0. 7001256+03 
0. 8532916+03 

0.480000E-01 

0 .740 8 20 6+03 

0.4950006-01 

0,7797996+03 

0.510000E-01 

0. 8172406+03 

0.5250006-01 

0. 8532916+03 

0.5283336-01 

0.8558376+03 

0.5316676-01 

0.8504176+03 

0.535000E-01 

0.8610316+03 

0.5358966-01 
C. 5713596-01 

0.861 7386+03 
0.891 5356+03 

0.538598E-0L 

0.6050766-01 

Q. 8638876+03 0.543152E-C1 

0.922 7246+03 0.6 542886-01 

THET = 0.6215156+02 

0.667558E+03 

0.9719136+03 

0.5496376-01 

0.8728886+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.405C006-01 

0.4481286+03 

0.420000 E-01 

0.4924396+03 

0.4350006-01 

0.5340COE+03 

0.4500006-01 

0.5731746+03 

0.4650006-01 

0.5250006-01 

0.61 02756+03 
0.7430086+03 

0. 4800006-01 

0.6455756+03 

0.4950006-01 

C. 6793196 +03 

0.5100006-01 

0.7117306+03 

0.5250006-01 

0.743008E+03 

0.5283336-01 

0.7449866+03 

0.5316676-01 

0.747C13E+03 

0.5350006-01 

0.7490896+03 

0.5358966-01 

0.5713596-01 

0.7496556+03 

0.7746066+03 

0.5385986-01 

0.6050766-01 

0.7513816+03 0.5431526-01 

0. 8023466+03 

THET = 0.6944956+02 

0. 7543616+03 

0.5496376-01 

0.7587546+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.4050006-01 

0.3791726+03 

0.4200006-01 

0.4167706+03 

0.4350006-01 

0.451640E+03 

0.4500006-01 

0*4839486+03 

0.4650006-01 

0.5250006-01 

0.5139766+03 

0.6169846+03 

0.4800006-01 

0.5419906+03 

0.4950006-01 

0.5682686 +03 

0. 5100006-01 

0.5931306+03 

0.5250006-01 

0.6169846+03 

0.5283336-01 

0.6183166+03 

0.5316676-01 

0.6197376 +03 

0.5350006-01 

0.6212466+03 

0.5358966-01 

0.571359E-ri 

0,6216656+03 

0.6432926+03 

0.5385986-01 

0.6229696+03 0,5431526-01 

THET = 0.7674756+02 

0.6252986+03 

0.5496376-01 

0. 6289056+03 

R 

T 

R 

T 

R 

7 

R 

T 

0.4050006-01 

0.3168216+03 

0.4200006-01 

0.3468946+03 

0.4350006-01 

0.3742936+03 

0.4500006-01 

0.3990526+03 

0.4650006-01 

0.5250006-01 

0.4209466+03 

0.4697656+03 

0.4800006-01 

0.4395656+03 

0.4950006-01 

C. 4544966 +03 

0. 5100006-01 

0.4649916+03 

C .5250006-01 

0.4597656+03 

0.528333E-01 

0.469840 6+03 

0.5316676-01 

0.470033E+U3 

0.5350006-01 

0.4703406+03 

0. 5358966-01 

0.4704396+03 

0.5385986-01 

0.4707886+03 0.5431526-01 

THET = 0.7674756+02 

0. 4715406+03 

0.5496376-01 

0.4729676+03 

R 

T 

R 

T 

R 

7 

R 

T 

0.4050006-01 

0.3168216+03 

0.4200006-01 

0. 3468946+03 

0.4350006-01 

0.3742936+03 

0.4500006-01 

0.3990526+03 

0.4650006-01 

0.5250006-01 

0.4209466+03 

0.469765E+03 

0.4800006-01 

0.4395656+03 

0.4950006-01 

0.4544966+03 

0.5100006-01 

0.4649916+03 

0.5250006-01 

0.4697656+03 

0.5283336-01 

0.4698406+03 

0.5316676-01 

C. 4700336+03 

0.535000E-01 

0.470340 E+0 3 

0.5358966-01 

0.4704396+03 

0.5335986-01 

0.4707886+03 0.5431526-01 

THET = 0. 8006076+02 

0.4715406 +03 

0.54963 76-01 

0.4729676+03 

R 

T 

R 

T 

R 

1 

R 

T 

0.40 50006-01 

0,2933926+03 

0.420000 E-01 

0.3238616+03 

0.4350006-01 

C. 3458866+03 

0.4500006-01 

0.3686066+03 

0,4650006-01 

0.5250006-01 

0.3891006+03 

0.4511976+03 

0.480000E-01 

0.407429E+03 

0.4950006-01 

0.4236036+03 

0.5100006-01 

0.4379376+03 

0.5250006-01 

0.5358966-01 

0. 4511976+03 
0.451834E+03 

0.528333E-01 

0.5385986-01 

0.451330 6+03 
0. 4520776+03 

0.5316676-01 

0.5431526-Cl 

C. 4515186+03 
0.452561E+03 

0.5350006-01 

0.4517616+03 
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THET 


0,833738E+02 


R 

T 

R 

T 

R 

T 

R 

T 

©•405000E-01 

0.272676E+03 

O.42P000E-O1 

0. 297927 B+03 

0.435000E-01 

C.321C61E+03 

0.450000E-01 

0.342284E+03 

0.465000E-01 

0.525CO0E-01 

0.361788E+03 

n.428020E^-03 

0.480900E-31 

0.379792E+03 

0.495000E-01 

0.396577E+03 

0.510000E-01 

0.412509E+03 

0,525000E-01 

0.535896E-01 

O.A^28020E+03 

0,^28661E<-03 

0.528333E-01 
0*538 598E-31 

0.428 163E+03 0.531667E-01 

0. 4288956+03 

THET = 0.866869E+02 

0.428353E+03 

0.5350COE-01 

0.428590 E+03 

R 

T 

R 

T 

R 

T 

R 

T 

0*405000E~01 

0.25^074E«-03 

0*420000E-01 

0.277344E+03 

0.435C00E-01 

0.2S8777E+03 

0.450000E-01 

0.318591E+03 

0.465000E-01 

0.525000E-01 

0.336990E+03 

0.400941E*03 

0.480000 E-Ol 

0.354183E+03 

0.495000E-01 

0.370398E+03 

0.510000E-01 

0.385891E+03 

•0.525000E-01 

0.535896E-01 

0.4C0941E^03 

0.4f)lA06E<-03 

0.528333E-01 

0.401062E+03 0.531667E-01 

THET = 0.900D00E+02 

0.401202E+03 

0.535000E-01 

0.401361E+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.405000E-01 

0,237058E+03 

0. 420000 EH)1 

0.258 486E+03 

0.435000E-01 

0-278285E+03 

0.450000E-01 

0.296657E+03 

O.A65000E-01 

0.525C00E-01 

0.313770E4-03 

0.372384E<-03 

0.480000E-01 

0.329777E+03 

0.495000E-01 

0.344813E+03 

0.510000E-01 

0.358992E+03 

0.525C00E-01 

0.37238!VE*03 

0.528333E-01 

0.372466E+03 

0.531667E-01 

C.372522E+03 

0.535000E-01 

0.372554E+03 



TEMPERATURES 

ON THE MIDDLE 

RECTANGULAR SECTION 







X = 0. 





Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.237058E+03 

0.150000E-02 

0.258486E+03 

0.300000E-02 

C.278285E +03 

0.450000E-02 

0.296657E+03 

0.600000E-02 

0.313770E+03 

0.750000E-02 

0.329777E+03 

0.900000E-02 

C.344813E +03 

0.105000E-01 

0.358992E+03 

C.120000E-01 

0.3723 84E+03 







0.1200006-01 

0.3723B4E+03 

0.123333E-01 

0.372 466 E+03 

0. 12666 7E-01 

0.372522E+03 

0.130000E-01 

0.372554E+03 




X = 0.361111E-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.21 5722c +03 

0.150000E-02 

0.2350 20 E+03 

0.300000E-02 

C.253149E+03 

0.450000E-02 

0.270184E+03 

0.600000E-02 

0. 2862016+03 

0.750000E-02 

D.3D1273E+03 

0.900000E-02 

0.315466E +03 

0.1050006-01 

0.328836E+03 

0.120000E-01 

0.341430E+03 







0.120CnOE-01 

0.341 430E +03 

0.123333E-01 

0.341 504 E+03 

0.126667E-01 

t. 34154 SE +03 

0.130000E-01 

0.341564E+03 




X * 0.7222226-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.198638E+03 

0.150000E-02 

0.216085E+03 

0.300000E-02 

C. 23260 5E +03 

0.450000E-02 

0.248240 E+03 

0.600000E-02 

0.263028E+03 

0. 750000 E-02 

0.277004E+03 

0.900000E-02 

0.29C196E+03 

0.105000E-01 

0.302626E+03 

0.120000E-01 

0.314304E+03 







0.120000E-01 

0.314304E+03 

0.123333E-01 

0.314371E+03 

0.126667E-01 

C. 31441 lE+03 

0.130000E-01 

0.314424E+03 




X = 0.1Q8333£-01 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.194523E+03 

0.150000E-02 

0.200 386 E+03 

0.300000E-02 

0.215473E+03 

0.450000E-02 

0.229814E+03 

C.600CO0E-02 

0.243434E+03 

0.750000E-02 

0. 2563536+03 

0.900000E-02 

0.268587E +03 

0.105000E-01 

0.280 148 E+03 

0.120000E-01 

0.291044E+03 







0.120000E-01 

0.291044E+03 

0.123333E-01 

0.291104E+03 

0.126667E-01 

0.291140E+03 

0.130000E-01 

0.29U52E+03 




X = 0.144444E-01 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.172691E+03 

0.150000E-Q2 

0.187192E+03 

0.300000E-C2 

C.201C25E+03 

0.450000E-02 

0.214215E+03 

0.600000E-02 

0.226780E+03 

0.750000E-02 

0.238737E+03 

0.900000E-02 

0.250094E+03 

0.105000E-01 

0.260863E+03 

0.1200006-01 

0.271051E+03 







0.120000E-01 

0.271051E+03 

0.123333E-01 

D.271106E+03 

0.126667E-01 

0.271138E+03 

0.130000E-01 

0.271149 E+03 




X = 0.180556E-01 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.1627C3E+03 

0.150000E-02 

0.176027E+03 

0.300000E-02 

0.188766E+03 

0.45000CE-02 

0.2C0946E+03 

0.600000E-02 

0.212576E+03 

0.750000E-02 

0.223669E+03 

0.900000E-02 

0.234233E +03 

0.105000E-01 

0.244275E+03 

0.120C00E-01 

0.253799E+03 







0.120000E-01 

0.253799E+03 

0.123333E-01 

0.253849E+03 

0.126667E-01 

0.253880E+03 

0.130000E-01 

0.253890E+03 




X = 0.216667E-01 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.154249E+03 

0.1500006-02 

0.166555E+03 

0.3000006-02 

0.178348E +03 

0.4500006-02 

0.189645E+03 

0.600P00E-02 

0.200455E+03 

0.7500006-02 

0.210787E+03 

0.9000006-02 

C. 2206446+03 

0.105000E-01 

0.230028E+03 

O.12OCOOE-01 

0.238941E+03 







0.120C006-01 

0.2389416+03 

0.1233336-01 

0.238988E+03 

0.126667E-01 

0.239016E+03 

0.130000E-01 

0.239025E+03 
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0,2527 78 E-01 


X = 


Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.l47088E<-03 

0.150000E-02 

3 .158517E+03 

0.300000E-Q2 

C.169496E +03 

0.450000E-02 

0.180038 E+03 

o.6onoroE-02 

o,i 2 acooe-oi 

0,1901476+03 

O.226262E+03 

0.7500006-02 

0.199826E+03 

0.900000E-02 

0.209073E+03 

0.105000E-01 

0.217886E+03 

0.120000E-01 

0.226262E+03 

0.123333E-01 

0.226305E+03 0. 12666 7E-01 

X = 0.288889E-01 

C.226331E+03 

0.130000E-01 

0.226340 E+03 

y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.141006E+03 

0.150000E-02 

0.151683E+03 

0.300000E-02 

0.161977E +03 

0.450C00E-02 

0.171894E+03 

0.60C000E-02 

0.120000E-01 

0.181433E+03 
0.215620E + 03 

0.750000E-Q2 

0.190587E+03 

0.900000E-02 

0.199344E+03 

0.105000E-01 

0.20 769 I £+03 

0,120000E-01 

0.215620E+03 

0.123333E-01 

0.215661E+03 0.126667E-01 

X = 0.325000E-01 

0.215685E+03 

0.130000E-01 

0.215693E+03 

Y 

T 

Y 

T 

Y 

T 

Y 

T 

0, 

0.135739E+03 

0.150000E-32 

0. 145778 E+03 

0.300000E-02 

0.155523E+03 

0.450000E-02 

0.164971E+03 

0,600000E-02 

0,12n00OE-01 

0.174113E+03 

0.206911E+03 

0.750000E-02 

0.182921E+03 

0.900000E-02 

0.I9I356E+Q3 

Q.105OQ0E-0L 

0.I99366E+03 

0,120C00E-01 

0.206911E+03 

0.123333E-31 

TEMPERATURES 

0.206949 E+03 0.126667E-01 

ON THE LOWER RADIAL SECTION OF 

THET * 0. 

0.206972E+03 

TUBE 

0.X30000E-01 

0.2U6979E+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.280000E-01 

0.135739E+03 

0.295000E-01 

0.145778E+03 

0.310000E-01 

C.155523E+03 

0.32S000E-01 

0.164971E+03 

0,34Q000E-0l 

C,400000E-ni 

0.1741136+03 

0.206911E+03 

0.355000E-01 

0.182921E+03 

0.370000E-01 

0.191356E +03 

0.385000E-01 

0.199366 E+03 

C.400000E-'C1 

0.206911E+03 

0.403333E-01 

0 . 206949 E+03 0 .406667E-01 

THET * 0.144952E+01 

C.206972E+03 

0.410000E-OI 

0.206979E+Q3 

R 

T 

R 

T 

R 

T 

R 

T 

0.280000E-C1 

0.13468PE+03 

0.295000E-01 

0.144559E+03 

0.310000E-01 

C.154145E+03 

0.325000E-01 

0.163449 E+03 

0.34CC00E-ni 

C.400000E-01 

0.172468E+03 

0.204789E+03 

0.3550006-01 

0.181176E+03 

0.370000E-01 

0.189522E+03 

0.385000E-01 

0.197429E+03 

o,4oorcoE“Oi 

0.204789E+03 

0.403333E-31 

0.204822E+03 0.406667E-01 

THET s 0.289903E+01 

0. 20483 6E +03 

0.410131E-01 

0 .20482964-03 

R 

T 

R 

T 

R 

T 

R 

T 

0.280rOOE-01 

0.133521E+03 

0.295000E-01 

0.143238E+03 

0.310000E-01 

0.152669E+03 

0.325000E-01 

0.161839E+03 

0.340fOOE-ri 

0.40P0G0E-01 

0.170752E+03 
0.202850E + 03 

0.355000E-01 

0.179387E+03 

0.370 00 OE-01 

0.187691E+03 

0.385000E-01 

0.195562E+03 

C.400000E~ni 

0.2028506+03 

0.403333E-01 

0 . 20 2 8 8 0 E +0 3 0 . 406 6 6 7E-0 1 

THET = 0.434855E+01 

C.2C2089E+O3 

0.410525E-01 

0.202872E+03 

R 

T 

R 

T 

R 

T 

R 

T 

O.20OCOOE-ni 

0.132269E+03 

0.295000E-01 

0.141818 E+03 

0.310000E-01 

0.151093E+03 

0.325000E-01 

0.160132E+03 

0.340PO0E-01 

0.400000E-C1 

0.168953E+03 

0.201124E+03 

0.355000E-01 

0.177546E+03 

0.370000E-01 

0.185860E+03 

0.385000E-0L 

0.193781E+03 

0.400000E-01 

0.201124E+03 

0.403333E-01 

0 . 20 1 1 55 E+03 0 .405 667E-0 1 

THET = 0.579806E+01 

0.2C1161E+03 

0.411184E-01 

0.201136E+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.280n00E-‘Cl 

0.13092 5E+03 

0.295000E-01 

0.140 298E+03 

0.310000E-01 

0.149412E+03 

0.325000E-01 

0.158319E+03 

0.340000E-01 

C.400000P-01 

0.167057E+C3 

Q.199625E+03 

0.355000E-01 

0.175633E+03 

0.370000E-01 

0.184012E+03 

0.385000E-01 

0.192081E+03 

0.400000E-CI 

0.199625E+03 

0.403333E-31 

0.199657E+03 0.406667E-01 

THET = Q.7247 58E+01 

C.199664E+03 

0.412108E-01 

0.199638E+03 

R 

T 

R 

T 

R 

T 

R 

T 

n.2800O0E-ni 

0.129492E+03 

0. 295000 E-01 

0.138677 E+03 

0.310000E-01 

0.147619E+03 

0.325000E-01 

0.156387E+03 

0.3400 00 E-*0l 
0.400000E-01 

0.165040E+03 

0.198350E+03 

0.355000E-31 

0.173617E+03 

0.370000E-01 

0.182113E+03 

0.385000E-01 

0.190436E+03 

0.400000E-01 

0.198350E+03 

0.403333E-01 

0.198384E+03 0.406667E-01 

THET = 0.869709E+01 

0.198397E+03 

0.413302E-01 

0.198377E+03 

R 

T 

R 

T 

R 

T 

R 

T 

O.2R000OE~01 

0.127970E+03 

0.295000E-01 

0.136952E+03 

0.310000E-01 

0.145707E+03 

0.325000E-01 

0.154320E+03 

0.34n000E-01 

0.400C00E-C1 

0.162875E+03 

0.197281E+03 

0.355000E-01 

0.171452E+03 

0.370000E-01 

C.180102E+03 

0.385000E-01 

0.1 88789 E+03 

0 .4000006-01 

0.197281E+03 

0.403333E-01 

0.197 321 E+03 

0.406667E-01 

0.157342E+03 

0.414769E-01 

0.197347E+03 
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THET * 0.101466E4-02 


R 

T 

R 

T 

R 

T 

R 

T 

0.2800r?1E-01 

0.126359E+03 

0.295000E-01 

0.135123E+03 

0.310000E-01 

0.142671E+03 

0.325000E-01 

0. 1521016+03 

O«3^P0OOE-Ol 

0,4000nnE-0i 

0*160527E*03 

0.196382E+O3 

0.355000E-01 

0.169080 E+03 

0.370000E-01 

C.177889E+03 

0.385000E-01 

0.1870376+03 

0.4000COE-01 

0.196382E+03 

0.4033 3 3 E-ai 

0.1964326*03 0.406667E-01 

THET * 0.1159616*02 

0.196472E+03 

0.416514E-01 

0.X96546E+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.28n000E-01 

■ 0.124661E+03 

O.Z95000E~01 

0.133191E*03 

0.310000E-01 

C.1415C5E+03 

0.325000E-01 

0.149716E+03 

0.34O0O0E-01 

n.400000E-01 

0.157963E^‘C3 

0.195563E+03 

0.355000E-01 

0.1664296*03 

0.370000E-01 

0.175346E +03 

0.385000E-01 

0.184994E+03 

O.AnnrooE-01 

0.195563E*03 

0.403333E~01 

0. 195648 E*03 0.406667E-01 

THET = 0.130456E+02 

C.195734E+03 

0.418543E-01 

O.195994E+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.280000E-01 

0.122882E*-03 

0.295000E-01 

0.131159E+03 

0.310000E-01 

0.1392C9E+O3 

0.32500 OE-01 

0.147154E+03 

0,3AOOOOE-Ol 

0.400000E-01 

0.155152E+03 

0-194382E4-03 

0.355000E-01 

0.163426E*03 

0.370000E-01 

0.172308E *03 

0.385000E-01 

0.182326E+03 

0.400000E-01 

0.194382E+03 

0.403333E-01 

0.194839E*03 0.406667E-01 

THET = 0.215961E+02 

C.195120E+03 

0.420862E-01 

0.195764E+03 

R 

T 

R 

T 

R 

7 

R 

T 

0.28f)0O0E-01 

0.110818E+O3 

0.295000E-01 

0.117222E*03 

0.310000E-01 

0.123C49E+03 

0.325000E-01 

0.128317E+03 

0.34f)000e-01 

0.40POOOE-01 

0.133007E«^03 

0.143583E+03 

0.3550006-01 

0.137053E*03 0.370000E-01 

THET * 0.301466E+02 

C.140330E +03 

0.385000E-01 

0.1426Z2E+03 

R 

T 

R 

T 

R 

T 

R 

T 

©•280P00E-01 

0.996175E*02 

0.2 950 00E-01 

0.104311E*03 

0.310000E-01 

C.108361E+03 

0.325000E-01 

0.U1794E+03 

0.340000E-01 

0.400COOE-01 

0.11461 5E+03 
0.119643E+03 

0.355000E-01 

0.1168 20E*03 0.370000E-01 

THET = 0.3869716*02 

0.118358E+03 

0.385000E-01 

0.U9337E+03 

R 

T 

R 

T 

R 

T 

R 

T 

0.280PPOE-01 

0.908152E^OZ 

0.295000 E-01 

0.9416956*02 

0.310COOE-01 

C.5699C3E+02 

0.325000E-01 

0.993073E+02 

0.340000E-ni 

0.400000E-01 

0.10U46E+03 

0.1042C0E+03 

0.3550006-01 

0. 102529 E*03 0.370000E-01 

THET * 0.4724766*02 

0.103481E+03 

0.385000E-01 

0. 104028 E+0 3 

R 

T 

R 

T 

R 

T 

R 

T 

0.280000E-01 

0.844081E^02 

0.295000 E-01 

0.8677106*02 

0.310000E-01 

0.887387E+02 

0.3250006-01 

0.903359E+02 

0.340000E-CI 

0.400000E-01 

0.915868E4^02 

0.936211E+02 

0.355000E-01 

0.925161 E+02 0.370000E-01 

THET s 0.557981E+02 

0.931482E+02 

0.385000E-01 

0.935082E+D2 

R 

T 

R 

T 

R 

T 

R 

T 

0.280000E-01 

0.799075E+C2 

0.295000E-01 

0.8155666*02 

0.310000E-01 

C.829264E+02 

0. 32500 OE-01 

0.840342E+02 

0.340000E-01 

0.400000E-01 

0.848983E+P2 
0.862 946E+02 

0.355000E-01 

0.8553776*02 0.370000E-01 

THET s 0.6434856*02 

C.859712E+02 

0.3850006-01 

0.862174E+02 

R 

T 

R 

T 

R 

T 

R 

T 

©•280OP0E-C1 

0.768027E+02 

0.295000E-01 

0.779480E+02 

0.310000E-01 

C.788994E+02 

0.3250006-01 

0.796683E+02 

0.340000E-01 

0.400000E-01 

0.802676E«^02 

0.812345E»02 

0.3550DOE-01 

0.8071066*02 0.370000E-01 

THET * 0.728990E+02 

C.8101C7E +02 

0.3850006-01 

0.811811E+02 

R 

T 

R 

T 

R 

T 

R 

T 

0.2800nOE-01 

0.746864E+02 

0.295000E-01 

0.7548206*02 

0.310000E-01 

C.761434E +02 

0.325000E-01 

0.766781E+02 

0.340POOE-01 

O.APOPOOE-n 

0.770949E4-02 

0.777674E<-02 

0.355000E-01 

0.7740 31E*02 0.370000E-01 

THET * 0.8144956*02 

C. 77611 8E +02 

0. 38500 OE -01 

0.777303E+02 

R 

T 

R 

T 

R 

T 

R 

T 

O.28PP00E-C1 

0. 73262 5E>02 

0.295000E-01 

0.738 195E+02 

0.310000E-C1 

C.74283CE+02 

0.3250006-01 

0.746580E+02 

0.34000PE-C1 

0.400000E-C1 

0.749505E«>02 

0.754221E4^02 

0.355000E-01 

0.7516676*02 0.370000E-01 

THET * 0.9000006*02 

0.753131E+02 

0.3850006-01 

0.753961E+02 

R 

T 

R 

T 

R 

T 

R 

T 

P,280POOE-Ol 

C.723253E+02 

0.295000E-01 

0.7272336*02 

0.310000E-C1 

0.730552E+02 

0.3250006-01 

0.733245E+02 

0.340000E-01 

0.40P0n0E-0l 

0.735348E*02 

0.738698E+02 

0.355000E-01 

0.7369006*02 

0.370000E-01 

C.737943E+02 

0.3850006-01 

0.738525E+02 
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TEMPERATURES ON THE LOWER RECTANGULAR SECTION 


X = 0. 


Y 

T 

Y 

T 

Y, 

T 

Y 

T 

0. 

C.60CC0OE-P2 
0-120C OOE-Cl 

n,723253E»02 

C.7353^8E<^02 

0.738698E+02 

0.150000E-02 

0.7500006-02 

0,7272336+02 
0. 736900 E+02 

0.300000E-02 

0.900000E-02 

0.730552E +02 
0.737943E +02 

0.450000E-02 

O.iOSOOOE-01 

0.7332456+02 

0,7385256+02 




X = 0.138889E-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.60PCOOE-02 

0.12000nE-Cl 

0.721021E^-02 

0,732428E+02 

0.735848E4'02 

0.150000E-02 

0.7500006-02 

0.724681E+02 

0.733971E+02 

0.300000E-C2 

0.900000E-02 

0.7278C2E +02 
0.735034E+02 

0.450000E-02 

0.105000E-01 

0.7303796+02 

0.7356506+02 




X = 0.2777 78 E-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

0.60C00CE-P2 

0.12OCPOE-0I 

0*719211E+02 

0.729939E^'0^ 

n.733335E4^02 

0.150POOE-02 

0.750000E-02 

0.722587E+02 

0.731449E+02 

0.300000E-02 

0.900000E-02 

0,7255116+02 
0.732507E +02 

0.4500006-02 
0.1 050006-01 

0.7279636+02 

0.7331306+02 




X * 0,4166676-02 




Y 

T 

Y 

T 

Y 

7 

Y 

T 

0. 

C.600000E-02 

P,1200n0E-Cl 

0.717741E+02 
0. 72784 7E <-02 
Q. 731171E+02 

0,150000E-02 

0,7500006-02 

0.720876E+02 

0.729312E+02 

0.300000E-02 

0.900000E-02 

0.723622E+02 
0,7303506 +02 

0.450000E-02 
0.1 050006-01 

0.7259506+02 

0.7309676+02 




X = 0.555556E-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

P. 

0.6n0O0OE~02 
0.120C OOE-Cl 

0.716558E+02 
0,726124= <-02 
0. 7293566+02 

0.150000E-02 

O.T50000E-02 

0.719495E+02 

0.727541E+02 

0.300000E-02 

0.9000006-02 

C. 7220876+02 
C.728551E+02 

0,4500006-02 

0,1050006-01 

0.7243036+02 

0,7291556+02 




X = 0.694444E-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0, 

0.600CO0E-02 

0.120f00E-01 

0.715628E+02 

0.724745E+02 

0,727885E+02 

0.150000E-02 

0.750000E-02 

0.718406E+02 

0.726116E+02 

0.300000E-02 

0.900000E-02 

0,7208726+02 

0.727099E+02 

0.4500006-02 

0.1050006-01 

0,7229926+02 

0.7276896+02 




X = 0.S33333E-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

O. 600000E-C2 

P. 120000E-01 

0,7149256+02 

0,7236926+02 

0,7267516+02 

0.150000E-02 

0.750000E-02 

0.717584E+02 

0.725024E+02 

0.300000E-02 

0.900000E-02 

0,7199516+02 
0,7.2 59 8 2£ +02 

0,4500006-02 

0.1050006-01 

0.7219956+02 

0.7265596+02 




X = 0.972222E-02 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0 . 

0.600000E-02 

O.12OC0OE-01 

0.714434E+02 

0,722949E+02 

0.725947E+02 

0.150000E-02 

0.750000E-02 

0.717008 E+02 
0.724253E+02 

0.300000E-02 

0.900000E-02 

0.7193C5E +02 
C.725192E +02 

0.4500006-02 

0.1050006-01 

0.7212946+02 

0.7257586+02 




X * O.llllllE-01 




Y 

T 

Y 

T 

Y 

T 

Y 

T 

0. 

C. 600 000 E -02 
0.120rOOE-01 

#■ 

0,7141446+02 

0.7225086+02 

0,7254666+02 

0.150GOOE-02 

0.750000E-02 

0,716667 E+02 0.300000E-02 

0.7237936+02 0.900000E-02 

X = 0.125000E-01 

0.7189236+02 
0.724720E +02 

0.4500006-02 

0.1050006-01 

0.7208786+02 

0.7252796+02 

Y 

T 

Y 

T 

Y 

T 

Y 

T 

0 . 

0.6C00P0E-02 

0.120DOOE-01 

0. 71404 7E +02 
0.7223626+02 
0.725307E+02 

0. 1500006-02 
0.750000E-02 

0.716554E+02 

0.723641E+02 

0,3000006-02 

0.900000E-C2 

C.718796E+02 

0.724563E+02 

0.4500006-02 

0.1050006-01 

0.7207406+02 

0.7251216+02 


REC= OOOCO FIL= 00002 
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Main Program 


COMMON /F0RMA/R1,R2,R3 ,YFINT,FINTHK,ANGTOP,HCOOL,TCOOLf HGAS»T€AS, 1 

1XLONGM,R2B,YFINB,ANGBOT,SUMAO»SUMAI 2 

DOUBLE PRECISION R1 ,R2 tR3 t YF I NT, FI NTHK, ANGTOP t HCOOL»TCOOL, HGAS, 3 

1TGAS,XL0NGH,R2B,YFIN3,ANGB0T,SUHA0,SUMAI ,TOL 4 

COMMON/FORAN/NRA,NRAM,NTA,NTAP,NRAl ,NRA1M,NRAIX,NRA1XM, 5 

INRB,NR8M,NTB,NTBM,NTD,NTDM 6 

COMMON /FORCN/ NXP ,NXP1 , NTQ, NTQl ,NZ * NZl ,NZX,NZXl 7 

COMMON /F0RWA/WME3A,WMEGA1 ,WMEG0,WMEGD1,WMEG8,WMEGB1 8 

DOUBLE PRECISION WME3 A ,WMEGA1 , WMEGO , WMEGOX ,WHEGB,H MEGBl S 

COMMON /FOR h:/ WMEGP ,HMEG PI , WMEGQ , WHEGQi 10 

DOUBLE PRECISION WMEGP tWMEGPl ,WMEGO,WMEGQl U 

DIMENSION 0ENT(12) 12 

EXTERNAL DUBIO 13 

14 

IPUNCH ^ 0 NO CARDS PUNCHED 15 

IPUNCH = 1 PUNCH 414 BCD CARDS 16 

TREAD = 0 SET INITIAL GUESS FOR U,S AND TEMPS IN RINIT 17 

IREAD = 1 RINIT READS 414 BCD CAROS FOR INITIAL GUESS 18 

NCOUNT MAXIMUM NUMBER OF ITERATIONS PROGRAM WILL CYCLE IS 

MCOUMT HEAT BALANCE PRINTED AFTER EVERY MCOUNT ITERATIONS 20 

ICOS = 0 HGAS = HGAS 21 

ICOS = 1 HGASm = HGAS ♦ COSITHETAdU 22 

R2 = R3 CASE WITH NO CLADDING 23 

24 

in READ! 5,103) (OENT(I) ,1=1 ,12) 25 

READ! 5,100) NRA ,NTA, NRB ,NTB , NTO , NXP, NTQ, NZ ,NZX , I PU NCH, IREAO, 26 

INCOUNT, MCOUNT, ICOS 27 

READ! 5,101) R1,R2,R3,YFINT,FINTHK,ANGTOP,HCOOL,TCOOL,HGAS, TGAS 28 

READ! 5,101) XLON3M,R2B,YFI N0 ,TOL 29 

READ! 5,101) WMEGA,WMEGA1 ,WMEG0,WMEGD1 ,WMEGB,WMEGB1 3C 

READ! 5,101) WMEGP, WMEGPl , WMEGQ, WMEGOl 31 

WRITE! 6,200) lOENt < I) , I =1 ,12 ) 32 

WRITE! 6,205) NR A,NTA ,NR8 ,NTB ,NTO ,NXP, NTQ, NZ ,N IX , IPUNCH, IRE AOtNCOUN 33 

IT, MCOUNT, ICOS 34 

WR ITE( 6, 201) R1,R2 ,R3 , YF I NT,F I NTHK ,ANGTQP 35 

WR ITE( 6,203) HCOOL ,TC OOL , HGAS , TGAS 36 

WR ITEI 6, 202) XLONGM ,R2B, YF I N8 , TOL 37 

WRITE! 6,206) WMEGA,WMEGA1 , WMEGD, WMEGOl , WMEG8, WMEGBl 38 

WR ITE( 6, 207) WMEGP, WMEGPl , WMEGQ , WMEG 01 39 

ANGBOT = ANGTOP 40 

ICOUNT = 1 41 

CALL CONSA(ICOS) 42 

CALL CONSB 43 

CALL R INIT(IREAD) 44 

1 SUMAI = 0.0 45 

SUMAO = 0.0 46 

CALL COMPAIICOS) 47 

CALL COMPB 48 

REL = ABSI (ABSISUMAI )-ABS<SUMAQ> )/ABS<SUMAI) ) 49 

IF(REL .LE. TOU GO TO 2 50 

IFIMOOCICOUNT, MCOUNT) .EQ. 0) WRITE(6,204) IC CUNT, SUMAI, SUMAO, REL 51 

ICOUNT * ICOUNT + I 52 

IF nCOUNT .ST. NCOUNT) GO TO 2 53 

G3 TO 1 54 

2 ICOUNT = ICOUNT - 1 55 

WR ITE( 6, 204) ICOUNT, SUMAI , SUMAO, REL 56 

CALL OUTPUTHPUNCH) 57 

GO TO 10 58 

100 FORMAT ( 1415) 59 

101 FORMAT ( 10D8. 5) 60 

1C3 F0RMAT(12A6) 61 

200 FDRMATUHK,25X12A6//1H ) 62 

2Cl FORMAT! IHK, 11 X2HRl,18X2HR2 ,1 8 X2H R3 , 15 X9HY-FI N-TOP, 11X9 HFIN-TH ICK, 63 

110X9HANGLE-T0P/6(5XG15.8) /IH ) 64 

20 2 FORMAT! IHK, 8X10HX-MID-RECT, 12 X6HR2-B0T, 11 X12HY-FIN-BQTT0H, 9X12HT0L 6 5 

lERANCE /4(5XG15.8)/1H ) 66 

203 FORMAT! 1HK,8X9HH-C00LANT,11X12HTEMP-CQ0LANT,12X5HH-GAS, 14X8HTEMP-G 67 

lAS/4! 5XG15.8) /IH ) 68 

204 F0RMAT!1HK,5X12HITERATI0N = I5,5X10HHEAT IN = G15. 8 ,5X 1 IHHEAT OUT 69 

1= G15.8,5X15HRELATIV= ERR = G15. 8//IH ) 70 

20 5 FORMAT ! IHK ,4X3HNRA ,5X3HNTA ,5 X3HNRB , 5 X3 HNT 8 ,5X 3HNT0, 5X3HNXP , 5X3HNTQ 71 

1, 5X2HNZ, 72 

16X3HNZX, 3X6HIPUNCH,3X5HIREA0 ,2 X6HNC0UNT ,2 X6HMC0UNT , 4X4HIC0S/ 14! 3X I 73 

25)/lH ) 74 

20 6 F3RMAT(1HK,8X8H0MEGA(A) ,12X9H0MEGA( All ,11 X8H0MEGA! 0 ), 12X9H0MEGA! 01 75 

1), llXaHOMEGAia ),12X9H0MESA(BU /6(5XG15.8)/1H ) 76 

20 7 F0RMAT(1HK,8X8H0MEGA!P) ,12X9H0M6GA!P1) ,11X8H0MEGA( Q ), 1 2X9H0MEGA! Q1 77 

l)/4! 5XG15.8) /IH ) 78 

END 79 
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SUBROUTINE CONSA(ICOS) I 

CO'<M3N/FORAS/NRA»NRAM,NTA,NTAP,NRA1 i NRAl M»NRA1X, NR AlXMf 2 

1MRB,NIRBM,NT8,NTBM,NTD»NT0M 3 

CD'^MON /FORCN/ NXP ,NXP1 , NTQ, NTQl »NZ t NZl ,N ZX,NZX1 A 

COMMON /FORA l/DRB,HCA,AlX,A2X,A3X,A4X,A5X,A6X,HGAfHCD» 5 

1HCAD,HGA0,HC3,HCAB,01X,D2X,D3X,AE,BE,DEX,DEY,EC#DENAD, 6 

20ENBZ 7 

DOUBLE PRECISION DRB,HCA , AIX , A2 X ,A3 X, A4, A5X ,A6X , HG A, HCO, 8 

IHCAD,HGAD,HCB,HCAB,OIX,02X,D3X,AE,BE,OEX,OEY,EC,DENAD»DENBZ 9 

COMMON /FORA2/AA a 5» ,AB (15) ♦ AC ( 15 ) » A1 AOS I t A1 8 ( 35 ) , AlC( 351, AD( 15», 10 

18D(15 ),CD{ 15) , AID (35) , BID (35 ), Cl 0(35 I , BAD (15) ,CAO( 15), 11 

2B1AD( 35),C1AD( 35) ,BA(15) ,BBU5) ,BC(15) ,B1A(20 ) tBlB<20l,81C(20)f 12 

3B8D( 15),CBD( 15) ,B1BDC20) ,C1 BO (20 ) , AP(15 ) , 8P( 15 ) , A1 P( 5 ) , 81P ( 5 ), 13 

A01BX< 15),02BX(15) lA 

DOUBLE PRECISION AA ,AB ,AC , A1 A , A1 8, AlC , AO, BD,C 0, A1 0, B1 D, CIO, 15 

1BAO,CAD,81AD,C1AD,BA,8B,3C,B1A,B1B,B1C,B80,CBO,B1BO,C180, 16 

2AP,BP , A1P,B1P,D1BX,02BX 17 

COMMON /FOUTA /THETA (25) ,RA(15) ,RAl X (35 ) ,R A1 (35 ) ,THET0(15), 18 

1THET8(15) 19 

DOUBLE PRECISION THETA ,RA ,RA1 X,RA1 ,THETO,THET B 20 

DIMENSION DRAX(35) 21 

DATA P 12/1.5707963268/ 22 

DOUBLE PRECISION DRAX,PI2 23 

COMMON /FORMA /R1,R2,RA2,Y1 ,RAH, PT ,HC , TC ,HG,TG, PX , R2 B,8Y 1, PT L, SUMAO, 2A 

ISUMAI 25 

DOUBLE PRECISION Rl ,R2 ,RA2 ,RAH ,P T,HC ,TC,HG ,TG ,PX,RZ8, PTL, SUM AO, 26 

ISUMAI, BYl 27 

DOUBLE PRECISION A4X, Y1 ,0R A,0T0, ST A, DTODR A,OT BORA, CUT Al, ENT A, AY, 28 

ICUTA2, THE8,X1,DRA2,OTORA,DTB2,YA,THA1,DTA,ENR,THEO,OTB,X2, DTA2, 29 

2R1DR,STB 3C 

COMMON /CLA/ UA (1 5) ,UA1 ( 35 ) , UD { 1 5) ,U01 <35 ) ,UB (15 ) , U8K 20 I 31 

DOUBLE PRECISION UA ,UAl , UD ,UOl ,U8,U8l 32 

COMMON /F0RANN/NRA11,NRAM1 ,NR8Ml 33 

DOUBLE PRECISION XA 34 

35 

GRID CONSTANTS FOR REGIONS A ,A1 ,D,D1 ,B,B1 36 

37 

XA := R2 + RAH 38 

CUTA2 * DATAN(Y1/XA) 39 

IF(RA2 .EQ. R2) GO TO 200 40 

YA s OS0RT(RA2*RA2 - XA*XA) 41 

CUTAl = OATANt YA/XA) 42 

GO TO 201 43 

20 0 THAI = PI2-PT-CUTA2 44 

CUTAl = CUTA2 + THAI ♦ .25 45 

201 THAI ^ PI2 - PT - CUTAl 46 

NTAl = NTA-1 47 

ENTA = NTAl 48 

NTAP = NTA+1 49 

DTA = THAl/ENTA 50 

DO II * UNTA 51 

AY = r-1 52 

I THETA (M = AY^DTA 53 

NRAM = NRA-l 54 

ENR = NRAM 55 

DRA = (R2-RD/ENR 56 

DO 20 1=1, NRA 57 

AY = I-1 58 

2C RA( I) = R1+AY*DRA 59 

THEO = CUTAl - CUTA2 60 

NTDM = NTO - 1 61 

ENTA = NTOM 62 

OTO = THEO /ENTA 63 

00 10 I = 1,NT0 64 

AY = I - 1 65 

10 THETD( I) = AY ♦ DTD 66 

THEB = CUTA2 67 

NR8M=NR8~1 68 

NTBM=NTB~1 69 

ENTA = NTBM 70 

ENR= NRBM 71 

DTB = THEB /ENTA 72 

ORB = RAH/ENR 73 

DO 30 I=l,NTB 74 

AY = I-1 75 

30 THFTB( I)=AY«OTB 76 

00 4C^ 1 = 1, NR8 77 

AY = I-l 78 
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AO RAK I)=R2 ^ AY^DRB 79 

SRAl = NTB + NR3i^ 8C 

NRA1M=NRA1-1 81 

DD 50 I=1,NTB 82 

J=NRA1+1-I 83 

50 RAl(J) = XA/DCOS ( THE 8 -THETBUH 8A 

RAIX { 1)=R2 0RB/2*0 85 

DO 60 I=2,NRBM 86 

60 RA1X( I)=RA1X(I-1J<-DRB 87 

Xl = THETB(NTB)-DT8/2.0 88 

DO 70 I=NRB,NRA1M 89 

RAIXU) = XA/DCnS(THEB -XII 90 

70 X1 = X1-DTB 91 

NRAIX = NRAl + NTDM 92 

IF<R2 -EO. RA2) NRAIX = NRAl 93 

IF(R2 .EQ. RA2) GO TO 202 94 

X2 = C0TA2 + DTO/2,0 95 

00 mo I == ItNTDM 96 

J = NRAl 4- I 97 

K = NTD-I 98 

RAKJ) = XA/DCOSfCUTAl - THETD(K)) 99 

J1 = J -I 100 

RAIX(JI) = XA/DC0StX2) 101 

iro X2 = X2 + DTD 102 

202 NRAIXM = NRAlX-1 103 

NRAll = NRAlXH-1 104 

NRAMl = NRAM-1 105 

NRBMI = NRBM - 1 106 

DRAX< 1 )=DRB/2.0 107 

DO 801 =2,NRA1XM 108 

80 DRAX( n^RAlXn )-RAlX U-1) 109 

DRAX(NRAIX) = RAUNRAIX) - RAlXtNRAlXM) 110 

111 

ITERATION CONSTANTS FOR REGIONS A, AI 112 

113 

CALL C IRCOUDRA,DTA,RA,AA,AB,AC,UA,NRA) 114 

HCA = HC Rl* OTA 115 

IF(R2 .eg, RA2J GO TO 203 116 

CALL C IRC0T(DRAX,DTA,RA1 , R Al X , Al A, Al B, AlC ,UAI , NRAIX ) 117 

A2X = AlCm 118 

203 AIX = AA(NRA) + AB(NRAM) 119 

A1C(NRA1X> = 2.0*A1A(NRAIXI + AIBINRAIXH) 120 

HGA = HG * RA2 ♦ OTA 121 

122 

ITERATION CONSTANTS FOR REGIONS D,01 123 

124 

CALL C IRCQUDRA,DTD,RA ,A0 ,B0 ,C0 , UO t NRA ) 125 

HCD = HC * R1 * OTO 126 

IF(R2 .EQ. RA2> GO TO 204 127 

CALL C IRCOT(ORAX,DTD,RAl ,R Al X»A1 0, B1 0, CIO ,UDl , NRAl XM) 128 

A6X = 2.0 * AlD(l) 4- BlOU) 129 

204 IFIR2 .EQ. RA2> A6X = HG ♦ R2 ♦ DTD 130 

A5X = AO(NRA» + BDtNRAM) 131 

AD(NRA) ^ AD(NRA)/2.0 132 

133 

ITERATION CONSTANTS FOR BOUNDARY BETWEEN REGIONS A,D AND Al, D1 134 

135 

DD 103 I = 1,NRA 136 

8AD(I» = <A3m ^ 8D(IM/2.0 137 

103 CADd) = AAm + AOin + RA(I) ♦ (OTA 4- OTO)/ORA 138 

CAD(NRA) = AACNRA)/2.0 + AO(NRA) + BAO(NRAM) 139 

CAOm ^ AA(1) ♦ AOU) + BAO(l) 140 

IF(R2 .EQ. RA2) GO TO 205 141 

BlAOdI = (AIBU) 4- B10(l))/2.0 142 

DO 104 I = 2,NRA1XM 143 

BIAD(I) (Aism + BlDUH/2.0 144 

1C4 ClADin = AlA(n + Atom + BlAOd) 4. BIAOU-U 145 

ClAOin = AlA(l) 4 - AlD(l) + 81ADU) 146 

20 5 HCAO = (HCA 4. HCD)/2.0 147 

HGAD = HGA/2.0 148 

IF(R2 .EQ. RA2) HGAD = HG AD4-A6 X/2. 0 149 

CIAD(NRAIX) = AIA(NRAIX) 4- B1 AO ( NRAl XM I 150 

151 

ITERATION CONSTANTS FOR REGIONS 6,81 152 

153 

CALL CIRCOL(ORA,OTB,RA ,8A , BB , BC , UB , NRAl 154 

CALL CIRC0T(DRAX,DTB,RA1 ,R Al X , B1 A, B1 B , BIC ,UB1 , NRAl Ml 155 

A4X * BIC (11 156 

XI = 2.0«DSIN(DTA/2.0) 157 
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A4X = 2,0* 81A(1)+ BIBU) 158 

A3X = BA(NRA) + BB(NRAM) 159 

BA(NRA ) = BA(NRA) /2.0 160 

HCB = HC*R1*)T8 161 

DO 92 1=1, NRA 162 

8BD(I)= (BDtn + 8B(n)/2.0 163 

92CBDCI) =AD(n +BA{I) + R A( I ) DTD+DTB) /DRA 164 

CBD(NRA) = A0(NRA)/2,D BA(NRA) + BBD(NRAM} 165 

CBDU) = A0(1) ^ BA(1) BBDdI 166 

167 

ITERATION CONSTANTS FOR BOUNDARY BETWEEN REGIONS 0,8 AND Dl, B1 168 

BIBDUI = (31D(l)+BlB(in /2*0 169 

DO 93 I=2,NRA1W 170 

81BO( n = (BIDU >+BlB tm /2.0 171 

93 ClBDtn = Alocn + B1D(I> 4- BlBOdH- BIBDU-IJ 172 

ClBDdl = A13(l) + BlOd) *■ BlBOCl) 173 

ClBD(NRAl) = AlO(NRAl) 4- BIBDINRAIM) 174 

HCAB = (HCD 4- HCB)/2.0 175 

DIX = AD(NRA) 4- BBO(NRAM) 4-BA{NRA) 176 

D2X = AlD(l) 4. B10(l)/2.0 177 

D3X = B1A(X» + BlB(l)/2.0 178 

OTD = DEY 179 

DO 110 1= 2,NRA1M 180 

OIBX (U =A1DU)+ (BlD(n+BlOU*in/2«0 101 

110 D2BX( n = BIA (n 4^(318(1)4- 81 B( I-l ) ) /2.0 182 

OIBX (NRAl) = AlO(NRAl) 4- B1 D ( NRAlM) /2«0 183 

028X CNRAl) * B18(NRAlM)/2.0 184 

IF(RA2 .NE* R2) GO TO 206 185 

D2X = A6X/2.0 186 

OlBX(l) = HG * DRAX(l) ♦ THAI 187 

IFdCOS .NE. 1) GO TO 207 188 

XI = OCOSCfHAl 4-THED) 189 

02X = D2X ♦Xl 190 

01BX( 1 1 = 018X11) * XI 191 

207 D3X * A4X/2.0 192 

DO 208 I = 2, NRAl 193 

D1BX( n = HG ♦ DRAXU) ♦ THAI 194 

208 IFdCOS *E0. 1) DIBXU) = OIBX(I) * XI 195 

196 

ITERATION CONSTANTS FOR BOUNDARY BETWEEN REGIONS B, P AND Bl,Pl 197 

198 

206 NXPl = NXP - 1 199 

ENTA = NXPl 200 

Xl= PX/ENTA 201 

AE = Xl/ORA 202 

BE = ORA/Xl 203 

OTBDRA = DT3/0RA 204 

DO 31 1=1, NRA 205 

APd) = (BB(I ) 4- AE)/2,0 206 

31 BP(I) = BA(I) 4- BE 4- AE 4- RA(I)*OTBDRA 207 

BPd) = BA(1)4- AP(1)4- BB/2.0 208 

DEX = BA(NRA) 4- 8E/2.04- APCNRAM) 209 

OENAD = Xl/DRB 210 

OEMS? = DRB/Xl 211 

DO 32 1=1, NR3 212 

AIP(I) = (BXB(n + OENA0)/2,O 213 

32 B1»(I) = BlAd) 4- DENBZ (RAKU^DTB 4- XD/ORB 214 

OEY = BlAd) + OENBZ/2.0 4- AlP(l) 215 

BIP(NRB) = B1A(NRB)4- DENBZ/2.0 4- AIP(NRBM) 216 

EC = (HCB4-Xl«HC)/2.0 217 

RETURN 218 

END 219 
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SUBROUTINE CONSB X 

common /FOR AN/NRA,NRAM,NTA,NTAP,NRA1 »NRAI M, NRAIX , NR AIX M* 2 

lNRB,NRBM,NTB,NT8M»NTDtNT0M 3 

common /FORCN/ NXP,NXPl,NTQtNTQl,NZ,NZl»N2X,N2Xl 4 

COMMON /FORCl /P X22 tHCP t PX2 f PXl t PX4 1 P2 X , PI X4 , PI X3 *HCPQt S XP»S 2P« 5 

1HCQ,01X,Q2X,HCQZ,QZ1,QZ2 ,HCZfZR3 ,HCRL,ZR1 ,ZR2 , ZR4, Z R5 , HCZR i PX3. 6 

iPlX2,PlXl 7 

DOUBLE PRECISION P X22 ,HCP t PX2 t PXl , PX4 , P2X , PIX4t PIX 3, HCPQt S IPt S2Pt 8 

lHC3,QlX,Q2X,HCQZfQZl,QZ2 ,HCZ ,ZR3 ,HCRL,ZR1 »ZR2 » ZR4, ZR5» HCZR. PX3 9 

2*P1X2*PIX1 10 

COMMON /F0RC2/AQU 5) ,6QU5> »C QU5 I , AX Q( 5 > » 8IQI 5 > tCX Q(5 ) » PQA( X5)» XX 

XP3B( 15>,PQA1(15) fPQBl (15) ,AQT (2U , BQT (2U ,COT (21) , Z A< X5 ),Z B ( X5 )• 12 

2ZC(X5), AOZ (15) ,CQZ(15) ,ZARL(15) tZ8RLU5) 13 

DOUBLE PRECISION AQ.BO ,CQ * AX Q, B1 Q,C1 Q, PQA .PQB , PQAl , PQBlt X4 

1AQT,8QT,CQT,ZA*ZB,ZC • A QZ »C QZ , ZARL, ZBRL X5 

COMMON /F0UTC/THETQ(20) f XP (20) , RQA( 15 ) »RQ(5) »THET Z ( 20 ) , 16 

1 XZ (20),QT (2(?),YP(20) ,YP1(5) 17 

DOUBLE PRECISION THETQ ,XP ,RQA,RQ,THETZ,XZ ,07, YP,YPX 18 

DATA P 12/1.570796328/ 19 

COMMON /FORMA /Rl,R2tRA2 »Y1 ,RAH .PTtHC, TC t HG t TGt PX , R2 B,BY It PT L, SUM AOt 20 

ISUMAI 21 

DOUBLE PRECISION Rl »R2 ,RA2 tRAH ,PTt HC ,TC.HGtTG,PX fR28t PTLt SUMAOt 22 

ISUMAItYl 23 

DOUBLE PRECISION BYl ,ORA ,DXP ,OQR ,X2 t THEZ , OELZ tX3,T HEQ. ORBt CTQ, OQLt 24 

UlDRfP I2 ,DZXtSUM 1,ENT,DRA2 tAY,STA,Xl ,R20R »ENZ t DPIX 2 ,X QR, RX 8t XQT » ZX 25 

COMMON /CLB/UQ(X5)»UZ<X5) ,CP(15) ,UP(15) ,CP1(5) ,U PX( 5 )♦ CZ ( 15 )t 26 

1UZR(15) 27 

DOUBLE PRECISION UQ,U Z ,CP »UP ,CP1 ,UP1 tCZfUZR 28 

COMMON /CLAB/ PXA ,PXAP t PXAl , PXAIP, ZRA» ZRP 29 

DOUBLE PRECISION PXA,PXAP, PXAl ,PXAXP fZRAf ZRP 30 

31 

GRID CONSTANTS FOR REGIONS P , PI t C» Q1 t ZL,ZR 32 

33 

R1B= R2B - R2 + R1 34 

THEQ= DATAN(8Y1/(R2B+RAH)| 35 

ENT = NR AM 36 

DRA =: (R2-RD/ENT 37 

ENT = NRBM 38 

ORB = RAH/ENT 39 

0RA2 = ORA/2.0 40 

NXPls* NXP-1 41 

ENT= NXPl 42 

OXP = PX/ENT 43 

00 7 I=1»NXP 44 

AY- I-l 45 

7 XP( AY*OXP 46 

DO 200 I = 1,NRA 47 

AY = I - 1 48 

200 YP( I) = AY * DRA 49 

DO 201 I = IfNRB 50 

AY = I 1 51 

201 YPl(I) = AY * DRB + YP(NRA) 52 

NTQ1= NTQ-l 53 

ENT= NTQ^l 54 

OTQ= THEQ/E.NT 55 

00 2 I=lfNTQ 56 

AY= I-l 57 

2 THETO( I)= AY*OTQ 58 

DO 3 1*1, NRA 59 

AY* I-l 60 

3 RQA( I )= R1B+ AY ♦ ORA 61 

R2DR = R2B + 0R8/2.0 62 

RQ(1) * RQA(NRA) 63 

DO 5 1 = 2, NR8 64 

5 RQ( I) * RQ(I-l) + DRB 65 

THEZ* PI2 + PTL - THEQ 66 

NZl* NZ-l 67 

ENZ* NZl 68 

DELZ* THEZ/ENZ 69 

DO 30 1=1, NZ 70 

AY* I-l 71 

3C THETZ( n* AY^DELZ 72 

ZX= Rl-RIB 73 

NZXl* NZX-1 74 

ENZ* NZXl 75 

DZX* ZX/ENZ 76 

DO 32 1*1, NZX 77 

AY* I-l 78 

32 XZ( I)= AY+OZX 79 

C 80 
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ITERATION CONSTANTS FOR REGIONS P,P1 81 

82 

PXl = dXP/DRA 83 

PX2 = DRA/DXP 84 

PlXl = DRB/DXP 85 

P1X2 = OXP/ORB 86 

PX22^ PX2/2.0 87 

HCP= HC+OXP 88 

PX3= PX1+PX2 84 

i>X4= 2.0*PX3 90 

P1X3= 2#0*(P1X1+P1X2I 91 

PlX4s= PlXl/2.0 92 

P2X= P 1X1+P1X2 93 

PXA = -PXI/PX4 94 

PXAP = PX2/PX4 95 

PXAl = -P1X2/P1X3 96 

PXAIP = P1X1/P1X3 97 

CALL RCOL(PXA,CP»UP,NRAM) 98 

CALL RC0UPXA1,CP1,UP1 »NR8M) 99 

100 

ITERATION CONSTANTS FOR REGIONS Q, Q1 101 

102 

CALL C IRCOL(ORA,OTQ,ROA,AQ,BQ,CQ,UQ,NRA) 103 

HCQ= HC^RIB^'DTQ 104 

OQR^ OTQ/DRa 105 

STA= DSIN<DTQ/2.0) 106 

X2= DRB/(2.0*S7A) 107 

AiaC n - X2/(2.0*R2B) 108 

81Qm= R20R*D0R 109 

C10(l):= X2/R2B + 81Q(1) 110 

00 6 I = 2,NR8 111 

AlQm= X2/RQ(n 112 

810(1 ) = 81Q(I-1) + DTQ 113 

6 C1Q<I)= 2.0»(A1QU) ♦ RQ(I)*OQR) 114 

01X=» AQ(NRA» + BQ(NRAM) 115 

32X= 2.0*A1Q(1> 4* 810(1) 116 

117 

ITERATION CONSTANTS FOR REGION ZL 118 

119 

CALL C IRCOL(DRA,OELZ,RQA,ZA,Z8,ZCfUZiNRA) 120 

ZC(NRA)= ZA(NRA) 4* ZB(NRAN) 121 

ZA(NRA)= ZA(NRA)/2,0 122 

HCZ= HC*R184'0ELZ 123 

124 

ITERATION CONSTANTS FOR REGION ZR 125 

126 

ZR1= ORA/(2*0=«'OZX) 127 

ZR2= OZX/ORA 128 

ZR3= 2.0*ZRl 129 

Z.R4= ZR3 4. ZR2 130 

ZR5= 2.0*ZR4 131 

ZRA = -ZR2/ZR5 132 

ZRP = ZP3/ZR5 133 

CALL RC0L( ZRA,CZfUZR,NRAM) 134 

HCZR* HC*DZX 135 

136 

ITERATION CONSTANTS FOR BOUNDARY BETWEEN REGIONS P,Q AND P1,Q1 137 

138 

00 55 I *1,NRA 139 

P3A(I) = (PXl+BQ(in/2.0 140 

55 PQBd) = PX24- AQ(I) 4* PX14- RQAdl^DTQ/ORA 141 

PQBd) * PX22 4. AQ(1) 4- PQA(l) 142 

00 56 1*1, NRB 143 

PQAKM * (P1X2 4- BIQ (I))/2.0 144 

56 PQBK n * PlXl 4- AIQ(I) 4. P1X2 + RQ (!)♦ OTQ/DRB 145 

X2 * OTO/2.0 146 

XI = R2B 4* RAH 147 

X3 = Xl-0R8/2*0 148 

SUMl = 0,0 149 

DO 90 1= 2,NTQ 150 

QT m * Xl/OCOS(THETQ(n ) 151 

XOT * QT(I) ♦ D$IN(THETQ(I n - SUMl 152 

SUNl = SUNl + XQT 153 

AQT(I) = XI -X3* OCOS(X2) 154 

AOT (I) = A3T(I)/XQT 155 

X2* X2 + DTQ 156 

90 B3T ( n = X3 ♦ DTQ/( QT(I) - RQ(NRBM)) 157 

AQT (NTQ 4 - 1 ) =0*0 158 

DO 91 I = 2, NTQ 159 
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91 COT(I) = AQTtn + AOT(H-U ♦ BQTdl 160 

CQTdJ s P1X4 +PQAUNRBM> + AQT(2) 161 

QT(1) == XI 162 

BOT(NTQ) * 8QT<NTQ)/2.0 163 

SIP = PX22 + A0<NRA)/2.0 + PQAINRAMI 164 

S2t> = P1X4 +A1Q(1) + PQAUl) 165 

HCPQ = <HCP +HCQ)/2,0 166 

167 

ITERATION CONSTANTS FOR BOUNDARY BETWEEN REGIONS QtZL 168 

03 92 I = 1,NRA 169, 

AOZ(I) = (BQ(I) + ZBm»/2*0 170 

92 C3Z(I) = AQ(I) + ZAin *■ RQAm*(0TQ+ DELZI^ORA 171 

C3zm = AQ(U + ZAU) + AQZm 172 

HC3Z = IHCQ +HCZ)/2.0 173 

QZl = AQ<NRAI/2.0 + ZAINRA) + AQZINRAMI 174 

QZ2 = AIQ (1) + B10(1>Z2.0 175 

176 

ITERATION CONSTANTS FOR BOUNDARY BETWEEN REGIONS ZL,ZR 177 

178 

03R = OELZ/ORA 179 

00 33 1= 1,NRA 180 

ZARL (I) = (ZB(I) 4^ ZR2)/2«0 181 

33 ZBRL in = ZAU) ZR4 4- RQAin* DQR 182 

ZBRUl) = ZARUl) +ZAa) ♦ ZRl 183 

ZBRL(NRA) = ZA(NRA) 4. 2R1 ^ ZARLINRAM) 184 

HCRL = (HCZ 4 - HCZR)/2.0 185 

RETURN 186 

END 187 


SUBROUTINE RCOL C A ,E L» U tN) 1 

2 

ROUTINE TO TRIANGULAR! Ze MATRIX CF COEFFICIENTS IN RECT. REGIONS 3 

4 

DIMENSION EL(1),U(1) 5 

DOUBLE PRECISION A,ELtU 6 

EL(2) = 1.0 7 

U( 2) = A 0 

DD 1 I = 3,N 9 

EL( n = 1.0 - A * UU-l) 10 

1 U( I) = A/ELU) 11 

RETURN 12 

END 13 


SUBROUTINE CIRCOUORA ,DTA , RA , AA , AB , AC , AO, N) 1 

2 

ROUTINE TO COMPUTE COEFFICIENTS AND TRI ANGULARIZE MATRIX FOR 3 

RADIAL REGIONS OF TUBE WALL 4 

5 

DIMENSION RA(l) ,AA(1) ,A8U) ,ACC1),AD(U 6 

DOUBLE PRECISION RA ,AA ,A8 ,AC ,AD,0RA2 ,DTA2 ,STA,DTORA,X 1, RIDR 7 

IfDRAfOTA 8 

0RA2 = DRA/2.0 9 

OTA 2 = DTA/2.0 10 

N1 = N-1 11 

STA = OSINOTA2) 12 

OTORA = OTA/DRA 13 

XI = DRA2/STA 14 

RIDR * RAID 4 - 0RA2 15 

AA(l) = X1/C2.0 ♦ RAUn 16 

AB( 1) = RIOR * DTDRA 17 

AC( 1) = Xl/RAU) + A8(l) 18 

DO 1 I = 2,N 19 

AA( I) = Xl/RAU) 20 

AB( I) = ABU-1) + OTA 21 

1 ACM) = 2.0 ♦ (AA(U 4 - RAID * OTORA) 22 

AD(2) =-AB(2)/AC(2) 23 

DD 2 I = 3,N1 24 

ACm = AC(I) ♦ ABU-1) ♦ AOU-1) 25 

2 AOm =-AB(n/ACm 26 

RETURN 27 

END 28 
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SUBROUTINE C I RCOT{ ORA X ,D TA ,RA1 ,RAX X ♦ A1 A, Al B, A1 Ct A1 0, N ) 1 

2 

ROUTINE TO COMPUTE COEFFICIENTS AND TRI ANGULAR IZE MATRIX FOR 3 

RADIAL REGIONS OF CLADDING AND FIN 4 

5 

DIMENSION ORAX(l),OMAX(n ,RAl (1) ,RA1 X(X ) , A1 A( 1 J , A1 BU > , A1C( 1 1, 6 

XAXDIXJ 7 

DOUBLE PRECISION ORAXtRAX t RAX X, AXA ,AX8, AXC,AX 0,XX* OT A 8 

NX = N-1 9 

XX = 2.C * 0SIN(DTA/2.O) XO 

AXA(X) = ORAX(X»/(XX * RAXCXH XX 

AXBIXI = RAXXIX) ♦ 0TA/(RAX(2) - RAX(XM X2 

AXCIX) = 2.0 * AXAIXJ + AXB(X) X3 

DO X I = 2,N X4 

AXAU) = ORAXm/IXX ♦ RAX (I M X5 

AlBCn = RAXxm ♦ DTA/(RAX(I+X) ~ RAX(D) X6 

X AlCm = 2,0 * AXAm + AXB(I-X) + AXBin X7 

AXD<2) =-AXB(2)/AXC(2) X8 

DO 2 I = 3, NX X9 

AXCdl = AX:m ♦ AXB(I“XI * AXO(l-XI 20 

2 A1D( I ) =~AXBU >/AXCU I 2X 

RET(JR.N 22 

END 23 


SUBROUTINE RINTIUREAD) X 

COMMON /FORAN/NRAfNRAM, NTA ,NTAP,NRAX ,NRAXM,NRAXX t NRAlXMt 2 

XNPB,NRBM,NTB,NTBM,NTO tNTDM 3 

COMMON /FORCN/ NXP »NXPX ,NTQ, NTQX #NZ f NZX t N ZX, NZXX 4 

COMMON /FORA/ A(X5,25I tAX (35,25) ,0(X5 ,X5 ) , BZ(20 ,X5) , 0( X5, X5 ), 5 

X DX(3 5,X5) ,T8A(25) ,TAI (2 5) ,TTA(2 5) ,TBB (X5 ) ,TB I( X5 ) ,TBD(X5) , 6 

2TDI(X5),T0D(X5) 7 

DOUBLE PRECISION A , AX , B ,BX ,D , DX , TBA , TAI ,TT A,F BB,T B I ,T BD, 0 

XTDItTDD 9 

COMMON /FORC/ P (20 ,X5) ,PX (20 ,5 ) , <3(X5 ,20 ) ,QX (X 5 ,20 ) , ZL( X5, 20 1, 10 

XZR ( 20,15) ,TP3(20) ,TPI (20) ,TB0(20) ,TQI(20 ) ,TZ 0( 20 ) , TZ BR( 20 ) XX 

DOUBLE PRECISION P,P1 ,Q,Q1 ,ZL,ZR ,TPB,TPI ,TBQ,TQI ,TZB,TZ8R 12 

DOUBLE PRECISION BZ X3 

IF (IREAO.EQ. X) GO TO 600 X4 

X5 

IREAD = 1 READ BCD CAROS CONTAINING U,S FOR ALL REGIONS 16 

AND TEMPS ON BOUNDARIES AND INTERFACES 17 

18 

INITIAL CONSTANT GUESS OF 1000 DEGREES FOR SAMPLE PROBLEM 19 

20 

DO X I = X,NTAP 21 

T8A( I ) = XOOO.O 22 

TAK I ) = XOOO.O 23 

TTA( I ) = XOOO.O 24 

DO 2 J = X,NRA 25 

2 A( J, I ) = .1802084 26 

DO X J = 1,NRA1X 27 

1 AX( J, I ) = .75271 28 

DO 4 J=X,NT3 29 

00 5 I=X,NRA 30 

5 D( I,J ) = A(I,X) 31 

DO 6 I=X,NRA1X 32 

6 Dl( I, J ) = Aid ,1) 33 

TBD( J ) = TBA( X) 34 

4 TDI(J) = TAI(X) 35 

DO 7 J=X,NTB 36 

DO 8 1=1, NRA 37 

8 B( I,J ) = A(I ,X) 38 

00 9 I=X,NRAX 39 

8Z( I, J ) = 4.7233244 40 

TBB( I ) = TOI ( X) 41 

9 TDD( M = TDUX) 42 

TBOU ) = TBA(X) 43 

7 TRKJ ) = TAI (X) 44 

DO XO J = X,NXP 45 

DO XX 1=1, NRA 46 

11 P(J,I) = A(I,1) 47 

DO 12 I=X,NRB 48 

12 Pl( J, I) = BZd ,X) 49 

Tf>B(J ) = TBA(X) 50 
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10 

TPl(J) = TAI(l) 

51 


03 13 

J=*1,NTQ 

52 


DO 14 

1 = 1, NRA 

53 

14 

3( I,J ) = A(I 

,1) 

54 


00 15 

1= 1,NR8 

55 

15 

31(1, 

J ) = BZ(I ,1) 

56 


TBQ(J ) = TBA(l) 

57 

13 

T3 I(J ) = TAI(l) 

58 


DO 16 

J=1,NZ 


59 


DO 17 

1=1, NRA 

60 

17 

ZL( If 

J ) = A( 1,1) 

61 

16 

TZB(J) = TBA(l) 

62 


no 18 

J=1,NZX 

63 


DO 19 

' 1=1, NRA 

64 

19 

ZR( J, 

I ) = A(I ,1) 

65 

18 

TZBR( J ) = TBA(l) 

66 


GD TO 

601 


67 





68 


READ 

INITIAL 

GUESS FROM BCD CAROS 

69 





70 

6C0 

CALL 

8CREAD 

(A(l,l) ,A(15,25) ) 

71 


CALL 

BCREAD 

(Al(l,l) ,A1(35,25) ) 

72 


CALL 

BCREAO 

(0(1,1), 0(15, 15) ) 

73 


CALL 

BCREAD 

(01(1,1) ,01(35 ,15) ) 

74 


CALL 

BCREAD 

(B(l,l), 8(15,15)) 

75 


CALL 

BCREAD 

(BZ(1,1) ,BZ(20,15) ) 

76 


CALL 

BCREAD 

(TBA(l) ,TBA(25) ) 

77 


CALL 

BCREAD 

(TAKl) ,TAI(25)) 

78 


CALL 

BCREAD 

(TTA(1),TTA(25)) 

79 


CALL 

BCREAD 

(T80(l ) ,T8D(15) ) 

80 


CALL 

BCR EAD 

(TDI(l) ,T0I(15) ) 

81 


CALL 

BCREAD 

(TDO(l) ,T00(15) ) 

82 


CALL 

BCREAD 

(TB8(1),T8B(15)) 

83 


CALL 

BCR EAD 

(TBI(l) ,TBI(15>) 

84 


CALL 

BCREAD 

(P(l,l) ,P(20,15) ) 

85 


CALL 

BCREAD 

(Pl(l,l) ,P1 (20,5) ) 

86 


CALL 

BCREAD 

(Q(l,l) ,Q(15,20) ) 

87 


CALL 

BCREAD 

(Q1(1,U ,Q1(15,20) ) 

88 


CALL 

BCREAD 

(ZL(1,1),ZL(15,20)) 

89 


CALL 

BCREAD 

(ZR(l,l),ZR(20,15n 

90 


CALL 

BCREAO 

(TPB(l),TPB(20n 

91 


CALL 

BCREAD 

(TPU1),TPI(20)) 

92 


CALL 

BCREAD 

(TBQ(l) ,TBQ(20) ) 

93 


CALL 

BCREAO 

(TQKl) ,TQI (20)) 

94 


CALL 

BCREAO 

(TZB(1),TZB(20) ) 

95 


CALL 

BCREAD 

(TZBR(1),TZ8R(20) ) 

96 

601 

continue 


97 


RETURN 


98 


END 
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SUBROUTINE COMPA(ICOS) 1 

CO'IMON /F0RAN/NRA,NRAM,NTA,NTAP,NRA1 ,NRA1M,NRA1X*NRA1XM, 2 

1NRB,NRBM,NTB,NTBM,NTD ,NTDM 3 

COMMON /FORA1/ORB,HCA,A1X,A2X,A3X,A^X,A5X,A6X, HGA,HCO* 4 

1HCA0,HGA0,HC3 ,HCAB ,01 X ,02 X ,03 X , AE, BE ,OEX,OEY, EC,OENAO, 5 

20ENBZ 6 

DOUBLE PRECISION ORB ,HCA , A1 X , A2 X ,A3 X , A4 , A5 X ,A6X , HG A, HCD, 7 

1HCAD,HGA0,H:3,HCAB ,Dl X,02X,03X,AE, be,dex,oey, EC,DENA0, OENBZ 8 

COMMON /F0RA2/AA(15) ,ABI15I ,AC(15) ,A1A(35) , AIB ( 35 ) , A1C( 35 ), A0( 15), 5 

1BD(15),CO<15),A10(35) ,B10(35 ),C10(35» ,BA0(15) ,CAO( 15), 10 

2B1AD( 35) ,C1A0(35) ,BA(15) ,B8<15) ,BC(15) ,B1A(27 ) , B1B( 20 ) , 81C ( 20 ) , 11 

3BBD( 15 ),CBO( 15 ),B 180(20) ,C18D (20) , AP ( 15 ) , BP(1 5 ) , A1 P( 5 ) , 81P ( 5 ), 12 

4D1BX( 15),D23X(15) 13 

DOUBLE PRECISION AA ,AB ,AC ,A1 A ,A1 8, AlC , AO, BO,C 0 ,A1D, B1 D, CIO, 14 

18AD,CAO,81AO,C1AO,BA,8B,3C ,B1 A ,B1B , B1 C , BBO ,CBO, BIBO, CIBO, 15 

2AP,BP ,AlP,BlP,0iaX,028X 16 

COMMON /FORA/ A ( 15 , 25) , A1 { 3 5 ,2 5 ) , B( 15 ,15 ) , BZ ( 20 ,15) ,D( 15, 15), 17 

1 01(35,15) ,TBA(25) ,TAI (25) ,TTA(25) ,T8B(15) ,TB I ( 15 ) ,T80( 15 ) , 18 

2T0K 15),TD0( 15) 19 

COMMON /FOUTA/THETA (25) ,RA(15) ,RA1X(35) ,RA1 (35 ) ,THET D( 15 ), 20 

1THETB(15) 21 

DOUBLE PRECISION A , A1 , B , B1 ,0 , 01 , TB A , TAI ,TT A,T 6B, T8 I ,T BO, 22 

1TOI,TDO 23 

DOUBLE PRECISION THETA ,RA , RAl X ,RA1 , THETO,THET B 24 

OIMEN SION COAA (25) ,CDNA1 (2 5) ,CONAI (25) ,C0NA1I (25), 25 
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ITBH 25 ), TTK2 5) ,TAI1 (2 5) , TAI2 (25 ) f TAI3 (25 ) fC0NBB(15 )i 26 

2CDNB U 15),C3N81Itl5) ,CONBC (15) ,C0NBD(15) fC0N0I(15) , CONDO 27 

3( 15)fCGN0lI(15)fX(35) ,F(35) 28 

DOUBLE PRECISION COAA ,C0NA1 ,CONA I ,C0NA1 I ,TBl r TT 1 fT A ll , 29 

ITA I2f TAI3,C3.M8B,C0NBI fCONBH fCONBCfCCNBDfCONDI f 30 

2C3ND0,C0ND1I fXfF 31 

C3MN0N /FORC/ P ( 20 ,1 5) ,Pl (20 ,5 ) , C(1 5 ,20 ) » Q1 (1 5 ,20 ) ,2L( 15t 20 ) * 32 

m ( 20,15) ,TP8 (20) ,TPI (20) ,T8Q(20) ,TQI (20) ,TZB(20 ) ,TZBR(20 ) 33 

D3U8LE PRECISION P ,P1 , Q, Q1 ,ZL , ZR,T PB ,TPI ,TBQ, TQI ,T Z B,TZ BR 3A 

DATA P 12/1.5707963268/ 35 

CO'IMON /FORMA /R1,R2,RA2 ,Y1 ,RAH , PT ,HC ,TC ,HG ,TG, PX , R2 8 , BY 1, PT L, SUM AD, 36 

ISUMAI 37 

DOUBLE PRECISION R1,R2 ,RA2 ,RAH ,PT,HC ,TC,HG ,TG ,PX,R2B, PTLfSUMAO, 38 

ISU'IAI, YIfBYI 39 

CO'IMDN /FORWA/WMEGAfWMEGAI , wmegd ,WMEG01 , wmegb , WMEG B1 40 

DOUBLE PRECISION WME3 A ,WMEGA1 , WMEGD , WMEGDl , WMEGB, W MEGBl 41 

DOUBLE PRECISION A4X , BZ, B2 ,HGA1 , H6AD1 ,BI X ,X1 , DTD, OET, X 2, B3, X3, P 1 2 42 

1,X4,X5,X6 43 

COMMON /CLA/ UA( 1 5 ) ,UA1 ( 35 ) , UD(15) t UDl (35 ) ,UB (15 ) , U8K 20 ) 44 

DOUBLE PRECISION UA ,UA1 ,U0 ,U01 ,UB, UBl 45 

COMMON /F0RANN/NRA11,NRAM1 ,NRBP1 46 

47 

COMPUTE U*S ON REGION A 48 

49 

DO 1 I = IfNRA 50 

1 A( 1,1) = A(I,3) 51 

DO 2 I = 1,NRA1X 52 

2 Al( I, 1) = Aid ,3) 53 

TAKl) = TAK3) 54 

TTA( 1) = TTA( 3) 55 

TBA(l) * TBA(3) 56 

CALL CON1(T3A,TB1,NTAP,COAA) 57 

CALL CQN2( TTA,TT1,NTAP,C0NA1) 58 

CALL C0N1(TAI,TAI1,NTAP,C0NAI) 59 

CALL CaN2(TAI,TAI2*NTAP,C0NAlI ) 60 

00 4 J 2,NTA 61 

JP = J+1 62 

J1 * J-1 63 

BIX = MCA » (TC - TBA(J) ♦ TBIU)) 64 

A( 1,J )*A(1,J ) + WMEGA»( (AAa)*{A(l ,Jl) + A(l,JP)M*AB(l)MrA( 2,J)+BIX)/ 65 

HAC(1)*HCA/:3AA( J) )-A(l, JH 66 

TBA(J) = TBA(J) - TBl(J) + A (1 , J )/COAA (J) 67 

SUMAO = SUMAO +HCA *(TC - TBA(J)) 68 

IF(J.EQ. 2) SUMAO = SUMAO/2.0 69 

F(2) =AB(l) *A(l,J) fAA(2)*(A(2,JPI +A(2,Jl)| 70 

03 5 I = 3,NRAM1 71 

5 F( n = AA(I) ♦ (Ad,JP)4- A(I,JD) 72 

F(NRAM) = AB(NRAM) ♦ A(NRA,J) + AA(NRAM) ♦ ( A ( NRAM , JP ) + A( NRAM, J 1 )) 73 

CALL CLINE (AB , AC, UA,F ,X,NRAM) 74 

OD 20 I = 2,NRAM 75 

20 A(I,J) = Ad,J) WME3A ♦ (X(I) - A(I,J)) 76 

77 

COMPUTE U»S ON MATERI AL I NTE RF ACE BETWEEN REGOONS A AND A1 78 

79 

IF(R2.NE.RA2) GO TO 101 80 

HGAl = HGA 81 

IF( ICOS.EQ.l) HGA1= HGA « DCOS( THETA( J1 ) » 82 

BIX=: HGAl * (TG-TAIU) ^ TAIKJ)) 83 

A(NRA,J)^ A(NRA,J) + WMEGA* ( ( AAINRA) ♦( AINRA, J1 )+ A( NR A, JP ) >/ 2 .0 84 

U A8(NRAM) I' A(NRAM,J) +BIX)/(A1X ^ HGAl/ CONA I( J ) ) -A( NRA, J ) ) 85 

TAUJ)= TAI(J) “ TAIKJ) i- A( NRA , J) /CONAK J) 86 

SUMAI = SUMAI * HGAl # (TG-TAKJ)) 87 

IF( J.EQ.2) SUMAI s SUMAI /2.0 88 

GO TO 4 89 

in DET = AIX/CONAIK JH- A2X/C0NAKJ) 90 

B1=AA{ NRA)*(A (NRA, JP) <-A(NRA, Jl) ) /2 • + AB ( NR AM) * A( NRAM, J I +A1A( D* 91 

1( Al( 1, JP )-»-Al{ 1, Jl) )*^A1B(1) *A1 (2, J) 92 

82 = TAI2(J) - TAIKJI 93 

XI = (-A2X * B2 Bl/CONAin J) )/OET 94 

X2 =(A1X =► B2 + 31/CONAI ( J) ) /DET 95 

A (N1A,J) = A (NRAfJ) WMEGA*( Xl-A(NRA,J)> 96 

A1(1,J) = A1(1,J)+ WMEGA1«(X2-A1 (1 , J) ) 97 

TAIU) = TAKJ) (A(NRA,J) /CONAK J) * A1 (1 , J )/CONAlI( J ) - (TAIKJ 98 

1) + TAI2(J)) )/2.0 99 

100 

COMPUTE U'S ON REGION A1 101 

102 

F(2) = A1B(1)*A1(1,J) +A1A(2) ♦(A1(2,JP) +A1(2,JD) 103 

DO 8 I = 3,NRA11 104 
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8 F(I) = AIA(I) *(Al(IfJP) + A1(I,J1)> 105 

F(N^AIXM) =A1B(NRA1XM) *A1(NRA1X,J) ♦ A1 A( NRAl X M) ♦( Al( NRAlX M, JP ) + 106 

lAKNRAlXM, Jl) ) 107 

CALL CLINE<A1B,A1C ,UAl ,FrX,NRAlXM) 108 

D3 21 I = 2»NRA1XM 109 

21 A1(I,J) = Ain,J) 4* WMEGAl * (X(U - Al ( I , JII 110 

HGAl = HGA 111 

IFdCOS .EQ. U HGAl = HGA ★ 0C0S(THETA( J1 U 112 

BIX = HGAl* <TG - TTA(J) 4- TTltJ)) 113 

AKMRAIX, J)=A1 (NR A1X,J) + WMEGAl* MAI A(NRA1X)*( Al( NR AIX , J1 > + Al( NRAIX 11 A 

1, )H-AlB(NRAlXM)*Al(NRAlXMt J) + BIX)/ ( AlC t NRAl X ) +HGA1/C0NA1 { J ) I- 115 

2AKNRA1X, J) J 116 

TTA(J) = TTA(J) “TTK JKAl (NRA1X,J)/C0NA1 ( J) 117 

SUMAI = SUMAl + HGAl *(TG-TTA(J)) 118 

IF(J.EQ. 2) SUMAI = SUMAI/2.0 119 

A CONTINUE 120 

121 

COMPUTE U'S ON BOUNDARY BETWEEN REGIONS A AND 0 122 

123 

BIX = HCAD * (TC - TBA(NTAP) + TBKNTAP)) 12A 

A(1,NTAP) =A(1,NTAP) WMEG A* M A A (1 ) *A( 1 ,NT A> + AD( 1 ) *D( 1 , 2 MBAD( 1) * 125 

1A(2.NTAPJ + BIX)/(CA0(1> 4- HCAD/COA A { NTAP M-A ( 1 , NT AP ») 126 

TBA(NTAP) = TBA(NTAPI - TBKNTAP) 4- A { 1 ,NT AP) /COAA( NT AP ) 127 

SUMAO = SUMAO + HCAD * (TC- TB A( NTAP) I 128 

DO 70 I=2tNRAM 129 

II = I-l 130 

70 A( UNTAP) =A(£,NTAP)+WMEGA*((AA(n*A(I,NTA)+ADm*On,2)+BAOUl )♦ 131 

1A(I1 ,NTAP) 4-BAD(I )*A(l4-l,NTAP)) /CADU) -A(I,NTAP)) 132 

IF(R2*EQ. RA2) GO TO 110 133 

DET = CAD(NRA)/C0NA1I(NTAP) 4- Cl AO (I I /CONAI(NTAP) 13A 

B2 = AA(NRA)*A(NRA,NTA)/2.04-AD(NRA)*0(NRA,2)4-8A0(NRAM)*A(NRAM,NTAP 135 

1)4-A1A( 1)*A1( 1,NTA)4-A1D(1)*DI (1 ,2 )4- B1 ADU ) *A1 ( 2 , NT AP ) 136 

B3 = TAI2 (NTAP) - TAIUNTAP) 137 

XI = (-ClAD(l) * B3 4- B2/C0NAIKNTAP) ) f DET 138 

X2 = (CAO(NRA) ♦ B3 4- B2 /CONAU NTAP ) ) / DET 139 

A(NRA,NTAP) = A(NRA,NTAP) 4- WMEGA * (X1-A{NRA,NTAP)» lAO 

Aid, NTAP) = Aid, NTAP) + WMEGAl * ( X2-A1 (1 , NT API I 141 

TAKNTAP) = TAI(NTAP) 4- ( A ( NRA ,NTAP) /CONAHNT AP) 4- Aid, NTAP) 142 

2/CONAlKNTAP) - (TAIKNTAP) 4- TAl 2 1 NTAP) I) / 2.0 143 

144 

COMPUTE U'S ON BOUNDARY BETWEEN REGIONS A1 AND 01 145 

146 

DO 71 I ^ 2,NRA1XM 147 

II = 1-1 148 

71 Aid, NTAP) =A1U,NTAP)+WM£GA1*( ( Al A( I)*A1 ( I , NT A) 4-A10( I) *Dl ( I, 2)+ 149 

1B1AD( II )*AI(I1 ,NTAP)+BlADn)*Al(l4-l,NTAP)MClAD( I )-Ald, NT AP ) ) 150 

HGAOl = HGAO 151 

IFdCOS .EQ. 1) HGAOl = HGAO ♦ OCOS ( THETA (NT A I) 152 

BIX = HGAOl* (TG - TTA(NTAP) 4- TTKNTAPI) 153 

AI(NRA1X,NTAP )= AK NR A1 X, NTAP) 4> WMEG A1 *M A1 A( NRAIX ) *A1( NRAIX, NT A ) 4- 154 

1B1AD(NRA1XM)*A1(NRA1XM,NTAP)4-BIX)/(HGA01/CONA1(NTAP) 4-ClAOtNRAlX) ) 155 

2-Al(NRAlX,NTAPn 156 

TTA(NITAP) = TTA(NTAP) - TTKNTAP) 4- Al (NRAIX, NTAP) /CONAKNTAP ) 157 

SU^^AI = SUMAI4.HGA01 *< TG -TTA(NTAP)) 158 

159 

COMPUTE U'S ON REGION D 160 

161 

GO TO 201 162 

110 HGADl = HGAO 163 

IFdC0S.6Q.il HGAOl * HGAO * DC OS ( TH£TA< NT A) ) 164 

BIX= HGAOl ♦ (TG-TAHNTAP) * TAIKNTAP)) 165 

A(NRA,NTAP)= A(NRA,NTAP) 4^ WMEGA * ( ( AA( NRA) * A( NRA, NT A) /2.0 ♦ 166 

lAD(NRA) * D(NRA,2) 4- B AO ( NRAM) *A (NRAM ,NT AP) 4-B IX )/ ( CAO( NRA) 167 

2 4. HGAOl/CONAI (NTAP)) - A(NRA,NTAP)) 168 

TAI(NTAP)= TAKNTAP) - TAIKNTAP) 4- A ( NRA ,NTAP )/ CONAK NT AP ) 169 

SUMAI = SUMAI ♦ HGADl ♦ ( TG-TAI ( NTAP) ) 170 

201 00 190 I =1,NRA 171 

190 0( 1,1) = AU ,NTAP) 172 

03 191 1=1, NRAIX 173 

191 01(1,1) = Aid, NTAP) 174 

CALL CONK TBO,TB1,NTD,CONBO) 175 

CALL CONI(TOI,TAI1 ,NTO,CONOn 176 

CALL C0N2( TDI,TAI2 ,NTD,CONDin 177 

DO 112 J=2,NTDM 178 

Jl = J -1 179 

JP = J 4- 1 180 

BIX = HCD*(TC-T80(J) + TB1( J)) 181 

0(1, J 1= 0(1, J)+WMEGO*((AOa)*(D(l, JP) +0(1, JD) +80(1)* 182 
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ID( 2f J )+BlX)/C:0(l)+HC0/C0NBD( Jn-0U»J)> 183 

TBD(J» ^ TBD(J» - TBl(J) ^ 0 (1 , J ) /CGNBD( J ) 184 

SUHA3 = SUMAO + HCO ( TC -TBO<J)) 185 

F(2) BD(1) *D(1,J1 i-AD(2» *CD(2,JP) +0(2,Jin 186 

D3 113 I = 3tNRAMl 187 

113 F(I> = A0(n ♦ (0(l»JP) ♦ 0«I,J1)> 188 

F(MRAM) = BO(NRAM) ♦D(NRA,J) A0( NRAMJ ♦«D(NRAM,JP) +D(NRAM,JD) 18S 

CALL CLINE(3D,CD,UD,F , X,NRAM) 19C 

03 26 I = 2,NRAM 191 

26 0( 1, = DUtJi + WMEGD*(X(I) - 192 

193 

COMPUTE U'S ON MATERIAL INTERFACE BETWEEN REGOGNS 0 AND D1 194 

195 

IF(R2.NE,RA2) GO TO 203 196 

HGAOl = A6X 197 

IFdCOS .EQ. 1) HGAOl A6X ♦ DCOS ( THETO ( J » + THETA(NTA)! 198 

BIX= HGAOl* (TG- TOMJ) * TAIKJM 199 

D(NRA,J) = D(NRAtJ) * WMEGD* ( ( AD ( NR A) ♦ ( D ( NRA, J1 H-0( NRA, J P ) 200 

1) + BO(NRAM) ♦ D(NRAM,J) + BIX)/(A5X + HGAOl/ CONDI C J ) ) - 201 

2D(NRA, J) ) 202 

TOHJ) = TDI(J) - TAIKJ) D ( NR A, Jl /COND U J ) 203 

SUMA1= SUMAI + HGAOl » (TG-TDUJII 2C4 

GO TO 112 205 

203 DET = A5X/C3N01I ( J) ^ A6X/C0NDUJ* 206 

Bl=AD(NRA)’*(0(NRA,JP)+D(NRA,Jin+BD(NRAM)*0(NRAM,JH-A10UI*(01(l, 207 

UP )+DUl,Jl) )+BiO(ll*01(2 , J) 208 

B2 = TA12IJ) - TAIKJ) 209 

XI = (- B2 * A6X * Bl/CONDIK J) l/DET 210 

X2 = (B2 ♦ A5X + Bl/CONOI ( J) I /DET 211 

D(NRA, J)= D(NRA,J)4-WMEG0*(X1“0<NRA,J) ) 212 

0K1,J) = 01(1, J) * WMEGOl (X2 - 01(1, J)) 213 

TDHJ) = TOKJ) + (D(NRA, J)/CQNDI( J) + 01 (1 , J ) /CONDII ( J ) - 214 

1 (TAIKJ) * TAI2( J))) /2.0 215 

216 

COMPUTE U»S ON REGION DI 217 

218 

K = NRAIX - J 219 

K1 * K-1 220 

F(2) = B1D(1)«DK1,J) +A1D(2) ♦(01(2, JP) +01(2, JD) 221 

03 114 I = 3,K1 222 

114 F(I) = AID(I) ♦ {OKI,JP) + D1U,JD) 223 

F(<) = BIO(IO*DKK+1,J) +A1D(K) ♦(01(K,JP) +01(K,JD) 224 

CALL CLINE(B1D,C10,UD1 ,F ,X,K) 225 

03 23 I = 2,K 226 

23 DKI,J) = DKIfJ) + WMEGOl* (X(I) - 01(I,J)) 227 

+1 226 

OKK, J )=DKK,J) + WMEGOl* ( (AID (K) *01 ( K , J-1 ) + B1D ( K-1 ) *DK K-1, J ) )/( AID 229 
1(< )+B1D(K-1) )-Dl(K,J) ) 230 

112 CONTINUE 231 

232 

COMPUTE U'S ON BOUNDARY BETWEEN REGIONS 0 AND 8 233 

234 

BIX = HCAB * (TC - TBD(NTO) + T6KNT0)) 235 

0(1, NTD ) =3(1, NTO ) +WMEG0* ( ( AD (1) *0(1 , NTOM) +BA( 1 )*B( 1 , 2) +8B0( 1 )* 236 

10(2, NTD ) + BTX)/(CB3(1) + HC AB/CONBO( NTD ))-0(l,NT0 )) 237 

TBO(NTO) = TBD(NTO) - TBKNTD) ♦ 0 ( 1 , NTO) / CON 80( NT 0) 238 

SUMAO = SUMAO + HCAB*(TC -TBD(NTO)) 239 

D3 170 I = 2,NRAM 240 

170 0( 1, NTO) = on ,NTD) + WMEG0* ((A0(n*0(I,NTDM)+BA( I)*B( I,2)+BBD( I-D* 241 

10( I“1,NTD)+3BD(I ) *0(1 + 1, NTD)) /CBD( D-D (I, NTO) ) 242 

TDD(l) = TDK NTO) 243 

244 

COMPUTE U AT MATERIAL INTERFACE POINT BETWEEN REGIONS D,D1,B1 245 

246 

CALL C0N3( TDD, TAI3,NRA1 , CONBC) 247 

IF(R2.NE,RA2) GO TO 204 248 

X3 = 01X+02X/CONOI (NTD) 249 

HGAOl = 03X +01BX(1)/CONBC(1) 250 

DET = X3/C0NBC(1) +HG AOl /C ONDI ( NTD ) 251 

BIX= D2X * (TG - TOKNTO) + TAIKNTO)) 252 

83 = D1BX(1)*(TG - TOKNTO) +TAI3(D) 253 

B1 = AD(NRA) * 0(NRA,NT0M) + BA(NRA) * B(NRA,2) + BBD( NRAH ) * 254 

1D(NRAM,NT0)+BIX + BlA(l) * BZ(1,2) + B1 B ( 1 ) *BZ ( 2,1 ) /2 ,0 + B3 255 

02 = TAI3(1) - TAIKNTO) 256 

XI = (B1/C0N8C(1) -B2 *HGAD1)/DET 257 

X2 = (82 *X3 + Bl/CONDI (NTD) ) /OET 258 

D(NRA,NTD)= D(NRA,NTD) + WMEGD* (Xl-O ( NRA , NTO ) ) 259 

BZ(i,l)= BZ(1,1) + WMEGBl ♦ ( X2 - 0Z(1,1)) 260 

TOKNTO) = TOKNTO) + ( 0 ( NR A ,NTO) /CONDK NTO) + BZ ( 1 , 1) / CONBC( 1 ) 261 
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1- (TAIKNTl) + TAI3(1))) /2.0 262 

SUMAI = SUMAI (D2X +01BXan*(TG -TDKNTDU 263 

03 205 I=2,NRA1M 264 

BfX= OieXd)* (TS + TAI3(I> - TDO(D) 265 

BZ(Itl)= BZntl) * WMEGBl * ((BIA(I) ♦ BZ(I,2) ♦ ( B1B< I-ll ♦BZC I-l. 266 

in + B ) /2.0 BIX)/<02BX(n + 01BX( n/CONBCUn-BZ( 267 

2I,l)i 268 

TDD«I)= TDDtn - TAI3(II + SZ ( I , i> /C0N8C ( 1 1 26S 

205 SU'IAI^: SUMAI + 013XU) * (TG-TOOCD) 27G 

8IX= DlBX(NRAl) ♦ ( TG * TAI3INRAI) - TOOINRAl) ) 271 

BZ(NRAl,l)=s BZINRAlfl) + WMEGBl ♦ ( ( B1B( NRAIM) ♦BZI NRAIM, I > /2 .0 ♦ 272 

lBIX)/(81BINRAlM)/2,0 ^ DIB X( NRAU / CDNBCI NRAl M -BZ( NRAl, 1) ) 273 

T00(>JRA1 )* TDDINRAl )- TAI3(NRA1 ) * BZ ( NRAl , 1 l/CONBCI NRAl ) 274 

SUMAI = SUMAI + OIBXINRAI ) * ( TG*TOO ( NRAl ) ) 275 

03 206 1=1, NRA 276 

206 B( I,l)= DU,NTD) 277 

G3 TO 207 278 

20A OET=D3X/(CONOI (NT0)*C0N01I (NTD) ) 01 X/CONOl I( NT 0) 4-D2X/ CONO H NTO )!/ 275 

1C3’^BC( 1) 280 

B1=ADINRA)*D(NRA,NT0M) + BBD (NRAMI NRAM,NTO ) +BA( NRAM^BINRA, 2 ) 281 

l+AlOd)* Dm,NTOMnBlO(l)*Dl(2,NTD)/2.D ^ B1 A( 1)*BZ( 1, 2 BIBIIJ *8 282 

2ZI2,nZ2.0 283 

B2= TAI2(NT0)~TAinNTD) 284 

B3 = TAI3(1) - TAIKNTOI 285 

Xl= (-B3 ♦ 33X/C0ND1I {NTD)+<B1/CCN01I (NT0)-B2*02XI/C0NBC( 1)1 286 

1/OET 287 

X2= ID3X*(B2“B3)/C0NDI (NT0)+(01X*B2 + Bl/COND I ( NTD ) )/ CONBCU ) ) 288 

1/ DET 289 

X3 = (- ID2X* 82-Bl/C0N01I(NT0n/C0N0I (NTD) * 03*( OlX/CON D1 1 280 

KMTOt 4- 02X/SONOI(NTOn)/OET 291 

D(NRA,NTO)= D(NRAfNTD) 4- WMEGO* ( XI -D ( NRA ,NTO )) 292 

01(1, NTD) = 31(1, NTD) t- WMEG01*(X2- D1(1,NT0)) 293 

0Z(1,1) = BZ(1,1)4- WMEGB1*(X3- BZ(1,D) 294 

TDI(MTD)= TOI(NT3)4-0 (NRA, NTOJ /CONDI (NTD) + D1 <1 , NT 01/ CONDI I( NTO) 295 

14- BZ( l,l)/C3NBCm- (TAIKNTD) 4- TAI2(NT0)+ TAI3{l)))/3,0 296 

TDD(1 ) = TDI (NTD) 297 

DO 90 I =1,NRA 298 

90 B( 1,1) = D(I,NTO) 299 

300 

COMPUTE U»S ON MATERIAL INTERFACE BETWEEN REGOCNS D1 AND BX 301 

302 

CALL C0N2 (TDD, TAI 2 , NRAl ,CONOD) 303 

D3 200 I =2,NRA1 304 

II = I-l 3C5 

DET = 01BX(I)/C0N3C(n 4- D20 X ( I ) /CONDO ( I) 306 

X3= 1.0 307 

IF (I.EQ. NRAl) X3=0. 0 308 

B1 = AIO(I) ♦ D1(I,NT0M) + BlOdl )♦ D1 d 1 ,NTD)/2.0 4- B1B( II )* 309 

1 BZ (II ,l)/2.0 + X3«^(BlDd )*Dl(l4-l,NTD)/2.04- B1 A( I) * BZ(I,2) + 310 

2 B1B( I ) ♦ 8Zd + l,l)/2.0) 311 

B2 = TAI3U) - TAI2(I) 312 

XI = (-B2*D23Xd) 4- Bl/CONBCd))/DET 313 

X2= (B2* D18X (1)4- Bl/CONDOd) ) /DET 314 

OK I, NTD) =01(1 , NTD )4-WMESDl *( XI -01 (I ,NTO) ) 315 

8Z(I,1) = BZ (1,1)4- WMEGBl* ( X2- 8Z(I,1)) 316 

200 TDO(I) = TDO(I) + (DZ d ,NTO ) /CONOOd > 4- 8Z d ,1 )/ CONBCU )- (TAI2 317 

1( I) 4- TAI3(I )) )/2.0 318 

319 

COMPUTE U»S ON REGION B 32C 

207 CALL C0Nl(TBB,TBltNT3,C0N8B) 321 

CALL CQN3( TBI ,TAI2 ,NTB ,C0NB1I ) 322 

CALL CONK TBI, TAIl, NTS ,C0N8I) 323 

DO 12 J=2,NTBM 324 

JP = J 4- 1 325 

J1 = J -1 326 

BIX = HCB*(T:-T8B( J)4-TB1( J) ) 327 

8( 1, J )= B( 1, J)4-WMEGB*n8Aa)*(8d, JP) 4-8(1, JD) 4-88(1)* 320 

1B( 2, J )+BIX)/(BC(1)4-HC8/C0NBQ( J) )-B(l, J) ) 329 

T8BU) = T8B<J) - T8KJ) 4- B ( 1 , J) /CONBB( J) 330 

SUMAO = SUMAO 4- HCB*(TC -T88(J)) 331 

F(2) = BB(1) *B(1,J) 4^BA(2) *(8(2,JPI 4-8(2,11)1 332 

03 13 I = 3,NRAM1 333 

13 F(I) = BAU) * (8d,JP) 4- BdtJl)) 334 

F(NRAM) =BB(NRAM)*8(NRA.J) 4-BA(NRAM) * ( 8 ( NRAM , JP) 4-B( NR AM, J I ) ) 335 

CALL CLINE(BB,6C,UB,F,X,NRAM) 336 

DO 24 I = 2, NRAM 337 

24 B( I,J) = Bd,J) 4^ WMEGB *(Xd) - 8d,J)) 338 

C 339 
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C COMPUTE U»S ON MATERIAL INTERFACE BETWEEN REGIONS 6 AND Bl 340 

OET * A3X/C3NBII ( J] + A4X/C0NBKJ) 341 

B1=BA( NRA)*(B (NRA , JP)+B(NRA, Jl) ) +BB( NR AM ) ♦BI NR AM, J ) +81 A( 1 > ♦( BZ ( 1, 342 

IJ’ )+BZ(l,Jl) J) 343 

B2 = TAI2U) TAIKJ) 344 

XI = (- B2 ♦ A4X + Bl/CONBin jn/DET 345 

X2 = (B2 * A3X 4. Bl/CONBKJUADET 346 

B(N1RA, J)= B(NRA, J)4-WMEGB*(XI“8<NRA,J) I 341 

BZ(1,J) = BZU,J) * WMEGBl * ( X2 - BZU,J)) 348 

TBKJ) = TBKJ) 4- {B(NRA, JI /CONBK J) 4- BZ ( 1 , J )/ CONBl M J ) - 349 

I (TAIKJ) 4- TAI2(J))>/2.0 350 

351 

COMPUTE U'S DN REGION Bl 352 

353 

K = NRAl-J 354 

K1 = K-1 355 

F(2) = BlB(l) *BZ(1,J) 4-BlA(2) *(BZ(2,JP) 4.BZ(2,Jin 356 

DO 14 I = 3,K1 357 

14 F( n = BlAd) ^(BZdfJP) 4- 8Z(I,J1)» 350 

F(X) = 81B(K>*BZ{K4-1, J) +BlA(Kl *<BZ(K,JPI +BZ<K,41U 359 

CALL CLINEIBIB ,B1C tUBl ,F ,X,K) 360 

DO 25 I = 2,K 361 

2 5 BZ( I, J ) = BZ(I ,J) + WMEGBl *(xn ) - BZtKJM 362 

<1 = K 363 

K=K+1 364 

8Z(K, J J=BZ(K, J)+ WME3Al*( (BIAUI *BZtK,Jl KBIBIKI )*8Z(K1 ,4K 365 

1/ (BIA(K) 4- 31B(K1 n-BZ(K,J)) 366 

12 CONTINUE 367 

368 

COMPUTE U»S ON BOUNDARY BETWEEN REGIONS B AND P 369 

370 

BIX = EC*( TC-TBB(NTB) 4- T8KNT8M 371 

BIlfNTB) = 3(1,NT3) 4 - WMEGB* ( (BA (1 » *8(1 ,NTBM» 4 - BE*P(2, 1 1/2 .0 372 

1+ AP(D* 6(2, NTB) 4 - BIX»/(BP(1I 4- EC/CONBB(NT 8) ) - BdtNTBH 373 

P( 1,1) = B(1,NTB) 374 

TBB(NTB) = T88{NTB)-TBl(NTB)4-B(lfNTB»/C0NB8(NT8) 375 

SUMAO = SUMAO 4- EC ♦( TC -TBB(NTBI) 376 

03 4C [=2,NRAM 377 

II = I-l 378 

B(I,NT0) = B(I ,NT8)4-WMEGB*U0A(n*B{I ,NT8M)4-BE*P<2, 1)4- 379 

lAPdl H'BUl ,NTB)+ APU )*8(l4.1,NTBn/8P(n- B(I,NTBM 380 

4C P( 1, I ) = B(I ,NTB) 381 

OET = DEX/CONBIHNTB) 4- DE Y/CONBU NTB) 382 

Bl= BA(NRA) * B(NRA,NTBM)4- BE*P (2 , NRA) /2*0 + AP(NRAM)* 383 

1B(NRAM,NTB) + B 1 A ( 1) <'8 Z (1 ,NTB «) 4-DE NB Z’S'Pl (2 ,1 ) / 2.04- AIP( 1 )*BZ( 2, NT B ) 384 

B2 ^ TAI2(NT3)-TAI1(NTB) 385 

XI = (-B2*DEY 4- Bl/CONBII (NTB) ) /OET 386 

X2 = (B2*DEX 4- B1/C0NBUNT8) )/OET 387 

B(NRA,NTB) = B ( NRA ,NTB ) 4 - WMEGB*( Xl-B( NR A,NTB) ) 388 

BZ(1,NTB) = BZ(i,NTB) 4 - WMEGBl ♦( X2-BZ(1 , NT 81) 389 

P(1,NRA) =3(NRA,NT8) 390 

PKl, 1) = BZ(1,NTB) 391 

TBKNTB) = T8I(NTB)4.(B{NRA,NTB)/C0NBUNTB)4- BZ ( l^NTB) /CONB 1 U 392 

INTB) - (TAIKNTB) 4 - TAT2(NT8) ) ) /2.0 393 

394 

COMPUTE U'S ON BOUNDARY BETWEEN REGIONS Bl AND PI 395 

396 

00 41 I = 2,NR8M 397 

II = I-l 398 

8Z(I,NT8)= BZ(I,NT8)4-WMEGBl*((BlA(n*BZ(I,NTBM) 4 - DENBZ+Pl ( 2, I) 4 - 399 

lAlPIIl )*BZ(I1 ,NTB)4-A1PU )*BZ<I4-1 ,NTB))/81P( n- BZ (I,NTB) ) 400 

41 Pl( 1, T )= BZ( I ,NTB) 401 

BZ(NRB,NTB)=BZ(NRB,NTB)4- WMEGBl^^tlBlAINRBl+BZ (NRB, NT8M)4-0EN8Z 402 

l*Pl(2,NRB)/2.0 4 - A1P(NRBM)*BZ(NRBM,NTB))/B1P(NRB)-BZ(NRB,NT8)) 403 

P1(1,NRB)= 8Z(NRB,NTB) 404 

RETURN 405 

END 406 


SUBROUTINE COMPB 1 

COMMON /FGRAN /NRA, NRA M, NTA, NTA P, NRAl ,NRA1M,NRAIX, NR AIXM, 2 

1NR8,NRBM,NT3,NTBM,NTD,NT0M 3 

COMMON /FORCN/ NXP ,NXP1 , NTQ, NTQl ,NZ ,NZ1 ,NZX,NZX1 4 

COMMON /for:/ P(20, 15) , pi (20,5) tQ(15 ,20) ,01(15,20) ,ZL(15,20), 5 
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III (20*15) ,TP3( 20) *TPI (20) ,TB0(20) ,TQI(20) *TZ B( 20 ) * T£ BR( 20 ) 6 

DOUBLE PRECISION P * PI , Q*Ql *Z L* ZR ,T PB *TPI ,TB0* T QI *T Z B, TZ BR 7 

CQV!M0N/F0RCl/PX22*HCP,PX2fPXl,PXA,P2X,PlX4,PlX3*HCPg*SlP*S2P* 8 

lHCa,QlX,Q2X,HCQZ,QZl,QZ2,HCZ,ZR3,HCRL,ZRl * ZR2 , ZR4, Z R5, HCZR * PX 3, 9 

1P1X2,P 1X1 10 

DOUBLE PRECISION PX22 ,HC P , PX2 , PXl * PX4 , P2 X , P1X4 , PIX 3* HCPQ* S IP* S2P* 11 

lHC3,QlXtQ2X,HCQZ*QZl,QZ2*HCZ*ZR3*HCRL,2Rl*ZR2*ZR4,ZR5*HCZR*PX3 12 

2*P1X2*P1X1 13 

COMMON /E0RC2/AQ( 15) *BQ (1 5 ) ,C Q(15 ) * A1 Q(5 ) , B1Q( 5 ) *C1 Q( 5 ),PQAI 15)* 14 

1P08(15)*PQA1< 15) »PQB1 (15) *AQT(2l ) >B0T(21 ) *C0r ( 21 ) * Z A( 15 ) *Z B ( 15 )» 15 

2ZC(15)* AQZ (15)*CQZ(15) ,ZARLU5)*ZBRL(15) 16 

DOUBLE PRECISION AQ*3Q *C Q* Al Q*8l Q*CI Q* PQA* PQB * PQAl * PQ61* 17 

1AQT,BQT,CQT* ZA , ZB * ZC , AQZ *CQZ , ZARL, ZBRL 18 

COMMON /F0UTC/THETQ(20) * XP(20) *RQA(15) ,RQ(5I *THETZ(20 )* 19 

1 XZ (20)*QT (20)*YP(20) ,YP1(5) 20 

DOUBLE PRECISION THETO , XP *RQA ,RQ *THETZ * XZ ,QT , YP *YP1 21 

DIMENSION T3K20) *TAI1 (20) *TAI2(20) *CONP( 20 ) * CONPI ( 2D ).C0NP1I( 20)* 22 

lC3Ng(20),CONQI (20) *C0NQ1I (20) *C0NZ(20 ) *CONZR( 20 ) *X ( 15 ) , F( 1 5 ) 23 

DOUBLE PRECISION TBl * TAI 1 * TAI 2 .CONP *C0NPI *C0NP1 1 ,CONQ* CONQ I* 24 

ICONQII *CONZ,CONZR*X*F 25 

COMMON /FORMA /R1,R2*RA2 * Y1 , RAH , PT*HC * TC ,HG *TG, PX * R2 B, BY 1, PT L* SUM AO, 26 

ISUMAI 27 

DOUBLE PRECISION R1*R2 *RA2 ,RAH,PT*HC ,TC,HG,TG*PX,R2B, PTL,SUMAO, 28 

ISUMAI* Y1,BY1 29 

COMMON /FOR WC/ WMEGP.WMEGPl *WMEGQ*WMEGQl 30 

DOUBLE PRECISION WMEOP * WMEGPl * WMEGO* WMEGQl 31 

DOUBLE PRECI SION B2,DET*X1 *61 X *B1 * X2 * X4,X5 32 

COMMON /CLB /UQ( 15 )*UZ( 15) *CP(15) *UP(15) *CP1(5) *UP1( 5),CZ( 15 )* 33 

1UZR(15) 34 

DOUBLE PRECISION UQ*U Z ,CP , UP *C PI *UP1 *C Z*UZR 35 

COMMON /CLA3/ PXA *PXAP ,PXA1 • PXAl P* ZRA * ZRP 36 

DOUBLE PRECISION PXA , PXAP , PXAl * PXAl P , ZRA , ZRP 37 

COMMON /FORANN/NRAll ,NRAM1 ,NRBM1 38 

39 

COMPUTE U»S IN RECTANGULAR REGION P 40 

41 

CALL CONK TP8 ,TB1 *NXP ,CONP) 42 

CALL CONI (TPI ,TAI1 *NXP, CONPI) 43 

CALL C0N3 (TPI ,TAI2*NXP,C0NPin 44 

DO 4 J=2*NXP1 45 

JP = J+1 46 

J1 = J-1 47 

8IX= HCP*(TC-TPB(J) TBKJ)) 48 

P( J*l)s P( J,l) 4- WM£GP^((PX22*(P(J1 ,1) + PUP ,1) )+PXl*P( J,2) ♦ 49 

1BIX)/(PX3 + HCP/C0NP( J))-P( J*l) ) 50 

TPB(J)« TPBU) - TBKJ) 4* P ( J ,1 ) /CONP U) 51 

F(2) = -PXA* P(J,1I 4- PXAP * (P(JP*2) 4* P(J1,2)) 52 

DO 5 I = 3*NRAM1 53 

5 F( n =: PXAP *(P(JP*I) 4- PU1,I)) 54 

F(NRAM) = -PXA * P(J*NRA ) +PXAP ♦(P(JP,NRAM) + P(J1,NRAM)) 55 

CALL RLINE(PXA,CP*UP*F ,X,NRAM) 56 

03 25 I = 2,NRAM 57 

25 P(J,I) = P(J*n 4- WMEGP * (X(I) - P(J,in 58 

59 

COMPUTE U'S ON MATERIAL INTERFACE BETWEEN REGIONS P*P1 60 

61 

DET = PX3/C0NP1I ( J) 4* P2X/C0NPKJ) 62 

PXl * P(J,NRAM) 4 . PX22*(P(J1 ,NRA) 4 - P(JP ,NRA)) + PlX 2 ♦ 63 

1P1(J*2) 4 - P1X4 * (PKJl ,1) 4 . PKJP ,1)) 64 

B2 == TAI2U) - TAIKJ) 65 

Xl= (81/CQNP1I(J) - B2*P2X)/OET 66 

X2=: (PX3*B2 4 - B1/C0NPI( J) ) /GET 67 

P(J,NRA)= P(J,NRA) + WMEGP * ( Xl-P( J , NRA)) 68 

f>l(J*l)= P1(J*1) 4 - WMEGPl * (X2-P1U*!)) 69 

TPI(J)s TPI(J) 4 - (P(J,NRA) /CONPI (J) + PI ( J ,1 )/ CONPI I( J ) - 70 

1 (TAIUJ) + TAI2(J)))/2.0 71 

72 

COMPUTE U»S ON RECTANGULAR REGION PI 73 

74 

F(2) = -PXAl *P1( J*1 ) 4.PXA1P*(P1 (JP,2) 4.PUJ1,2)) 75 

DO 6 I = 3,NRBM1 76 

6 F(I) = PXAIP *(P1(JP,I) 4^ P1(J1,I)) 77 

F(NRBM) = -PXAl* Pl(J,NRB ) + PXAIP* ( Pi ( JP * NRB M) 4.pi( J X, NRBM ) ) 78 

CALL RLINE(RXA1*CP1*UP1*F,X*NRBM) 79 

03 26 I = 2, NRBM 80 

26 P1(J,I) ^ Pl(J*n 4. WMEGPl * (X(I) - P1(J,I)) 81 

P1(J,NRB) = P1(J,NRB) 4-WMEGPl ♦ ( ( PI X2 *P1 ( J * NR BM) +P1X 4*( PI ( J 1* NRB) 82 

14-P1(JP,NRB)))/P2X - P1(J,NR8)) 03 

4 SUMAO s SUMAO +HCP ♦ ( TC - TPB(J)) 84 
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85 

C3MPUTE C*S ON BOUNDARY BETWEEN REGIONS P AND Q 86 

87 

BIX = HCPQ*(TC- TPB(NXP) + TBKNXP)) 88 

P(^^XP,1)::: P( N XP , 1 > + WMESP*({PX22 *P(NXPi»l)+AQm ♦0(1,2)+ PQA(l) 89 

1* P(NXP,2) BIX>/(PQB(1)+ HCPQ/CONPCNXP) )- P(NXP,1M 90 

TPB(NXP)= TP3(NXP|-TBl (NXP ) + P ( NXP,1 ) /CONP ( NXP ) 91 

SUMAD = SUMAO + HCPQ ♦(TC - TPB(NXP)) 92 

3(1, U= P(NXP,1) 93 

DD 50 I=2,NRAM 94 

II = I-l 95 

PINXP, l)= P(NXP,n + WMEGP^((PX2+P(NXPl,n + AQU)^Q(l,2)+ 96 

1P3A(I1 )*P(NXP,U » + PQA(n*P<NXP,I + in/PQBCn-P(NXP,n) 97 

50 Q( I,U= P(NXP,I ) 98 

GET = S1P/C3NP1I (NXP) 4- S2P/CONPI( NXP) 99 

B1 = PX22 *P(NXP1,NRA)+ AQ ( NRAI ♦OI NRA ,2 ) /2.0+ PQAC NRAM ) ♦ 100 

1P( NXP,NRAM)+P1X4*P1(NXP1 ,1)<-A1QU) *01 (1,2 )+PQAl(l) ♦Pl(NXP,2) 101 

B2 = TAI2(NXP)- TAIUNXP) 102 

XI = (-S2P*B2 + B1/C0NP1I(NXP))/DET 103 

X2 = (SIP *B2 + Bl/CONPl (NXP))/OET 104 

P(NXP,NRA)= P(NXP,NRA)+ WMEGP^( Xl-P( NXP,NRA) ) 105 

Pl(MXP,n= P1(NXP,1) + WMEGP1«(X2-P1 (NXP,in 106 

Q(NRA,n = P(NXP,NRA) 107 

TPI(NXP) = TPI(NXP) *■ (P(NXP,NRA) /CONPKNXP) + PI ( NXP, 1 )/ CONP II 108 

l('JXP) - (TAIKNXP) + TAI2(NXPn ) /2.0 109 

110 

C3MPUTE U»S ON BOUNDARY BETWEEN i^EGIONS PI AND Q1 111 

112 

31(1,1)= P1(NXP,1) 113 

D3 51 I=2,NRBM 114 

II = I-l 115 

PKNXP ,I )= Pl(NXP,nf WMEGP1*( (PI Xl*Pl(NXPl,n+AlQ( U+OK I, 2) 116 

1+P3AUI1 )*P1(NXP,I1 )+ PQAl (n^Pl (NXP,I + 1))/PQB1( I) - PKNXP,!)) 117 

51 Dl( It n= PKNXP, n 118 

PKNXP ,NRB )=P1( NXPtNRB ) + WMEGPl*( ( PI X4^P1 ( NXPl , NRB) +PQAKNRBM ) ♦ 119 

IP K NXP ,NRBM)+AQT(2I* Q1 (NRB,2M /CQT( K -PI ( NXP , NRB) ) 120 

Q1(NRB,1)= PKNXP, NRB) 121 

122 

COMPUTE U*S ON REGION Q 123 

124 

CALL CaNl(TBQ,TBl,NTQ,CONQ) 125 

CALL CONI (T3I ,TAI1 ,NTQ,CONQI) 126 

CALL C0N3(TQI,TAI2,NTQ,CQNQin 127 

03 1 J«2,NTQ1 128 

JP = J + 1 129 

J1 = J -1 130 

BIX= HCQ^(TC-TBQ( J) + TB1( J) ) 131 

3( ItJ )= Q(1,J)+WMESQ4^( (AQ(1)*(Q(1, J1 )+QU,JP )H“BQ(1)^ 132 

13(2,J)+BIX)/(CQ(l)+HCQ/CONQ(Jn-Q(l,J)) 133 

TBO(J)= TBQ(J)-TBl(J)+Qa,J)/CONQCJ) 134 

SUMAO = SUMAO + HCQ * ( TC - TBQ(J)) 135 

F(2) = BQ(1) *Q(1,JH'AQ(2)*(Q(2,JP) + Q(2,JD) 136 

D3 2 I = 3,NRAM1 137 

2 F( n = AQ(I) ♦ {Q(I,JP) + QUfJU) 138 

F(NRAM) =BO(NRAM) *Q( NRA, J) «^AQ( NRAM) ♦{ Q(NRAM, JP) +Q( NRAM, Jl) • 139 

CALL CLINE(BQ,CQ,UQ,F,X,NRAM) 140 

DO 20 I = 2,NRAM 141 

20 0( I,J ) = Q(It J) + WM=SQ ♦(XU) - QU ,JM 142 

143 

COMPUTE U«S ON MATERIAL INTERFACE BETWEEN REGIONS Q AND Ql 144 

145 

DET = QIX/CONOIK J)+Q2X/C0NQI < J) 146 

Bl= AQ(NRA)*(Q(NRA,JP )+Q(NRA,Jl ) ) /2,0+BQ(NRAM) *Q( NRAM, 147 

IJ) +A1Q(1)*(Q1(1,JP )+01(l,Jl n+BlQ(K+Ql(2, J) 148 

82= TAI2(J>-TAIKJ) 149 

Xl= (-Q2X*B2+B1/C0NQ1I ( JM /DET 150 

X2= (Q1X«B2^B1/C0NQK jn /DET 151 

0(NRA,J)* Q(NRA, J)4-WMEGQ«(Xl-(3(NRA,jn 152 

31(1, J)= QKl, J)«^WMESQ1«(X2-Q1(1,J)) 153 

TO I( J) = TQK J) + CQ(NRA, J) /CONQI ( J » + Q1 ( 1 , J) /CONQl I ( J ) - 154 

1( TAIK J)*TA!2( J)) ) /2, 0 155 

156 

COMPUTE U»S ON REGION Ql 157 

158 

00 3 I = 2,NRBM1 159 

II = I-l 160 

3 3KI,J)= Ql(I,J) + WME3Ql+((A10(n«(Ql(I,Jl )+Ql(I,JP ))+ 161 

IBlQdl )*QKI1 ,JK^BlO(n*Ql(K-l ,J))/C10(I)-Q1(I,J) ) 162 
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J)=Q1(NRBM, J»+WMEGQ1*( (AI0(NRBM»*(01( NRBM, J1 l + QK NRBM»JP 163 

1) )+BlQ(NRBMl)*Ql(NRBMl , J)^'BQT( *01 (NRB,J) » / ( A1Q( NRBM ) ♦a.O ♦81Q( NRB 164 

2'il )+BQT( J) )-Ql(NRBM,J) ) 165 

Ql(NRB,J)s 31(NRB,J) + WMEGQl* ( ( AQT ( J) *Q1 ( NRB , J1 ) + A0T(JP )* 166 

IQltNRBfJP ) + BQTIJ)^ QHNRBM,J))/CQTU) -QUNRB,J)) 167 

1 CONTINUE 168 

169 

COMPUTE U«S ON BOUNDARY BETWEEN REGIONS Q AND ZL 170 

171 

BIX = HCQZ*(TC-TBQ(NTQ) * TBKNTCI) 172 

QdfMTQ) =: 3(l,NTQdWM£GQ*MAQa )*Qa,NTQl)+ZAm*ZL(l»2) 173 

U AQZ(l) * QiZfHTQ) ^ BIXI/ICQZU) + HCQZ/CONQ( NTQ) )-Q( 1,NTQ) ) 174 

ZL( 1, I ) = Q(1,NTQ) 175 

TBOINTQ) - T8Q(NTQ» - TBKNTQI + Q( 1 » NTQ) /CON Q( NT Q ) 176 

SUMAD = SUMAO + HCQZ *(TC - TBQ( NTQ) ) 177 

DO 200 I=2fNRAM 178 

II = 1-1 179 

Q{ I,NTQ) » QU jNTQM-WMEGQ* U AQCI )*G(I , NT Q1 dZ A { I ) *ZU I,2J+ 180 

lAOZUl )*Q(I1 ,NTQ)MQZ(I )*Q(dl ,NTQHZCQZm-Qn,NTQn 181 

20C ZL( I, n = Q(I »NTO) 182 

DET = QZI/CONQIKNTQ) QZ2 /CONQM NTQ) 183 

B1 = AQ(NRA) ♦ Q(NRA,NTQ1) /2.0+ AQZ( NRAM) *QCN R AM,NT Q) +Z A( NRA) 184 

I* ZL(NRA,2)+ A1Q(1)*Q1 (1,NTQ1)+ B1 Q(1 ) *Q1 ( 2 , NTQ)/2 .0 185 

B2 = TA12(NTQ)-TAIUNTQ) 184 

XI = (-QZ2*B2 + 81/CQNQinNTQ))/0ET 187 

X2 = (QZ1*B2 + Bl/CONQUNTQ) ) /DET 188 

Q(NRA,NTQ)=3(NRA,NTQ) + WMEGQ*(X1-Q(NRA,NTQ) I 189 

31(1, NTQ) = Qia,NTQ) + WMEGQl*(X2-Ql(l,NTQ)) 190 

T3UVITQ) =TQI(NTQ) + {Q(NRA,NTQ) /CONQI ( NT Q) +Q1 ( 1 , NT Q) / CONQl U NTQ) 191 

1-(TAI1(NTQ) + TAI2(NTQ) ) 1/2.0 192 

ZL(NRA,l)sQ(NRA,NTQ) 193 

194 

COMPUTE U*S ON BOUNDARY OF REGION Q1 195 

196 

DO 201 I - 2,NRBMl 197 

II = 1-1 198 

201 Ql( I,NTQ)=Q1(I ,NTQ)^‘WMEGQ1»UA1Q(I)*Q1 (I ,NTQ1 ) + 01Q( II )*Q1(I1 ,NTQ 199 

l)/2.0+BlQ( I) *31(1+1, NTQ) /2.0) / ( AIQ ( I ) + ( B1 Q( I ) +B1QI II ) )/2.0l- 200 

2Ql(I,NTQ)) 201 

Q1(NRBM,NTQ) =Q1(NRBM,NTQ) + WMEGQ1*( ( A1 Q ( NRBM) *Q1 ( NR BM, NTQl) +81Q(NRB 202 

1M1)*Q1 (NRBMltNTQ) /2. 0+8QT( NTQ) *Q1 ( NRB , NTQ) ) / ( A1Q( NR8M ) +B1Q( NRBM 1 )/ 203 

22.0+BQT(NTQ) )-Ql(NRBM,NTQ) ) 204 

01(NR8,NTQ)=Q1(NRB,NTQ) + WMEGQ1*( (AQT { NTQ) *01 ( NRB, NT Q1 ) + 205 

IB3T(NTQ)*Q1CNR8M,NTQ) ) / (A QTt NTQ) +8 QT ( NTQ) ) -Q1 ( NRB, NTQ) ) 206 

207 

COMPUTE U*S ON REGION ZL 208 

209 

CALL CONI (TZB,TB1,NZ,C0NZ) 210 

D3 30 J=2,NZ1 211 

JP = J + 1 212 

J1 = J -1 213 

BIX == HCZ*{T:-TZB( J)+ TBl(J)) 214 

ZLll,J)= ZL(l,J) + WMEGQ*(tZA(l)*IZLtl ,J1 )+ZL(l,JP ) )+ZB< 1) *ZLI 2, J ) 215 

1+BIX)/(ZC(1)+HCZ/C0NZ(J)) - ZLU.JH 216 

TZ8(J)= TZ8(J) - TBK J) + ZL(I, J)/CONZ(J) 217 

SUMAO = SUMAO + HC Z * ( TC - TZ8(JM 218 

F(2) = ZB(1) *ZL(1 ,J) + ZA(2) * ( ZL (2 , J P) +ZL (2 , J 1 ) ) 219 

03 31 I = 3,NRAM1 220 

31 F(I) = ZA(I) *(ZLU,JP) + ZLdtJin 221 

F(NRAM) = Z3(NRAM)*ZL(NRA,J) +ZA(NRAM) *( ZL ( NR AM, J P) +ZL(NRAM, JD) 222 

CALL CLINE(ZB,ZC,UZ,F,X,NRAM) 223 

09 22 I = 2,NRAM 224 

22 ZL(I,J) = ZL(I,J) +WMEGQ *{X(I) - ZL(I,J)) 225 

ZL(NRA,J) = ZL(NRA,J) + WMEGQ *((ZA(NRA) * (ZL(NRA,J1 ) + ZL(NRA, 226 

UP )) + ZBINRAM) * ZL(NRAM,J) )/ZC(NRA) - ZL(NRA,J)) 227 

30 CONTINUE 228 

229 

COMPUTE- U*S ON BOUNDARY BETWEEN REGIONS ZL AND ZR 230 

231 

BIX = HCRL *(TC-TZB(NZ)+ TBKNZ)) 232 

ZL(1,NZ) ^ ZL(1,NZ ) + WMEGQ ♦ ( (Z ARL( I ) *ZL (2 , NZ ) +Z A( 1 ) * 233 

IZL(1,NZ1)+ ZR1*ZR(2,1) + 6 IXI / ( ZBRL( 1 ) +HCRL/CQNZ( NZ U - 234 

2ZL(1 ,NZ)) 235 

TZB(NZ)= TZB(NZ) - TBl (NZ) + ZL( 1 , NZ) /CONZ ( NZ ) 236 

SUMAO = SUMAO + HCRL * ( TC - TZB(NZ)) 237 

03 34 I=2,NRAM 238 

II = I-l 239 
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34 ZL( I,NZ)= ZL ( I fNZH'WMEGQ+C (ZARL( I) »NZ I ♦ 240 

IZARKIl ) *ZL(I1 ,NZ) * ZR3*ZR(2 ,n + ZA(I»*ZL( I,NZ1 nZ 241 

2ZBRL< I ) - ZL(I,NZ») 242 

ZL(MRA,NZ) = ZL<NRA,NZ) + MMEGQ ♦ ((ZA(NRA) * ZL(NRA,NZ1I + 243 

IZRl * ZR(2,NRA) * ZARL(NRAM) ♦ ZU NR AM ^NZ ) ) / Z BRLCN RA) -ZLCNRAiNZI 244 

2) 245 

D3 35 1=1, NRA 246 

35 ZR(1, [)= ZL(I,NZ) 247 

248 

COMPUTE U'S ON REGION ZR 24S 

250 

CALL CONI (TZBR,TB1 ,NZX,CONZR) 251 

no 32 J=2,NZX 252 

Jl = J~1 253 

JP = J+1 254 

IF(J .LT. NZX) GO TO 28 255 

DO 29 I = 1,NRA 256 

29 ZR( JP, I) = ZR( Jl,n 257 

28 BIX = HCZR *(TC - TZBR(J) ^ TBl < J) J 258 

ZR(J,1)= ZR(J,U + WME3Q*({ZR1*(ZRU1 tU + ZRUP ,1)) + 259 

1 IRZ^lRiJfZ) + BIX)/(ZR4 + HCZR/CONZRUU -ZR( J , 1 ) ) 260 

TZ0R(J)= TZ3R( J)-TB1( J) ZR ( J,U /CONZRI J) 261 

SUMAQ = SUMAO + HCZR *(TC - TZ8R(JM 262 

F(2) = ~ZRA ♦ZRIJ,!) +ZRP ♦(ZR(JP,2) ♦ZRU1,2M 263 

DO 33 I = 3,NRAM1 264 

33 F(I) = ZRP ♦ (ZR(JP,n + ZR{J1,IU 265 

F(NRAM> = -ZRA ^ZR(J,NRA ) +ZRP ZR ( JP, NRAM) +ZR( Jl, NRAM) > 266 

CALL RLINEIZRA ,CZ,UZR ,F,X,NRAM) 267 

DO 27 I = 2,NRAM 268 

27 ZR(J,I) = ZR(J,n WMEGQ ♦(X(I) - ZRU,IM 269 

ZR(J,NRA) = ZR(J,NRAJ +WMEGQ* { ( ZR2*ZR( J,NRAM) +ZR1 ♦(ZRUl, NRA I + 270 

1 ZR( JP,NRA)) ) /ZR4 - ZRU,NRAH 271 

32 CONTINUE 272 

RETURN 273 

ENO 274 


SUBROUTINE CONI ( TOLD , U ,NTI P,C OND ) 1 

2 

CONDUCTIVITY FOR STAINLESS STEEL 3 

TREFl = 70.0 FOR SAMPLE PROBLEM 4 

5 

DIMENSION TOLD (1) ,U(l ) ,C0ND(50> ,B(6,5) ,TKI6) ,TEG(6 ) 6 

DOUBLE PRECISION TOLD , U,COND ,B , TK, TEG , SUM ,TNEW 7 

DATA (B< 1, J>,J =1, 5) Z. 803803070-4, .637514990-6, -.142016130-8,. 1083 8 

171240-11,0. DOZ 9 

DATA (B(2,J), J=1,5)Z.115204890-3,.19479961D-6, -.1240770-9, .624877 10 

11D-13,Q.00Z 11 

DATA (B(3,J), J =1 , 5) Z.40113088D-3 ,-.384859870-6, .339124520-9,-. 99 12 

1203914D-13,. 1D638424D-16Z 13 

DATA (TKU),J = 1,3) Z70.00,400.00, 1200. DOZ 14 

DATA IN/IZ 15 

IF ( IN .EQ. 1) GO TO 6 16 

4 00 1 J =1,NTIM 17 

IF (TOLOU) .GT. 70.D0) GO TO 20 18 

TDLD( J ) = 70.00 19 

C3ND( J ) = .11841930-3 20 

U(J) = 0.0 21 

GO TO 1 22 

20 K= 1 23 

DD 2 L = 1,2 24 

LI = L + 1 25 

2 IF(TOLOU) .GT. TK( LI) ) K = LI 26 

8 C3ND(J)= B<K,1) 4- TOLDU) ♦ CB(K,2) ♦ TOLO(J) ♦ (B(K,3) ♦ TOLD( J )♦ 27 

1(B(K,4) + TOLDU) * B(K,5)))) 28 

N = K 29 

SUM = 0.0 30 

DO 7 I = 1,N 31 

7 SUM = SUM + TEGU ) 32 

TNEW = SLM 4- TOLD(J) * (B(K,1) + TOLDU) ★(B<K,2)Z2.0 ♦ TOLDU) ♦ 33 

1 (B(K, 31/3.0 + TOLDU) ♦ <B(K,4)/4.3 ♦ TOLOU) * B(K, 5 )/5.0 ) ) ) ) - 34 

2 TK(N) ♦ (B(K,1) 4 - TK(N) (B(K,2)Z2.0 ♦ TK(N) ♦ (B(K, 31/3.0 ♦ TK 35 

3 U) ♦ (B(K, 41/4.0 4- TK(N) * B ( K ,5 ) Z5 .0 ) ) ) ) 36 

U( J ) = TNEW/CONDU) 37 

1 CONTINUE 38 

RETURN 39 
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6 TEGdl = 0.0 40 

DD 5 I = 1,2 41 

II = I + 1 42 

5 TEGim = TK(in ★ (B(I , 1 ) ^ TKUn ♦ (0(1 ,21/2.0 ♦ TKII 1 1 * ( 43 

1 B( I ,3)/3.0 + TK(in ♦ (8(1 ,41/4.0 ^ TKIII) * 8(1 ,5)/3.0))n - 44 

2TK(II * (8(1 ,1) + TK(n ♦ (B(l ,21/2.0 * TK( II * (8(1 ,31/3.0 > 45 

3 TK(I) ♦ (B(I ,41/4.0 ^ TK(II * 8(1 ,51/5.0)))) 46 

IM = 2 47 

GO TO 4 48 

END 49 


SUBROUTINE C 0N2( TOLO , U ,NTI )«,C0ND ) 1 

2 

CONDUCTIVITY NICKEL 3 

TREF2 = 230.0 FOR SAMPLE PROBLEM 4 

5 

DIMENSION TOLDdl ,U(U ,COND(50I ,8(6,5) ,TK(6) ,TEG(6) 6 

DOUBLE PRECISION TOLO , U,COND ,B ,T K, TEG ,SUM ,TNEW 7 

DATA (B( 1,J),J=1,5)/ .4045C27D-2 ,-.134431050-4,. 19014014D-7,-. 8813 8 

141260- n,-. 63 8609670-1 5/ 9 

DATA(B(2,J),J=1,5)/ . 165812 570-2, -. 272994110-5, .3 44585970-8,- .238 10 

1504690-11,. 67843 9540-1 5/ 11 

DATA ( B( 3, J) ,J = 1,5)/. 432941010-3, .22290446 D-6 ,-.17 147530-9 ,. 120048 12 

178D- 12, -.236501 90-16/ 13 

DATA ( TK(J),J = 1, 31/230.00, 660.00, 1360.00/ 14 

DATA IN/1/ 15 

IF ( IN .EQ. 1) GO TO 6 16 

4 D3 1 J=1,NTIM 17 

20 K= 1 18 

D3 2 L=l,2 19 

LI = L + I 20 

2 IF (TOLD( J).GT.TK(Lin K=L1 21 

8 CONO(J)= B(K,1) I- TOLD(J) « (B(K,2) + TOLO(J) # (B(K,3) + TOLOU)* 22 

1(B(K,4) + TOLD(J) ♦ 8(K,5)m 23 

N = K 24 

SUM == 0.0 25 

00 7 I = 1,N 26 

7 SUM - SUM + TEGU) 27 

TNEW ^ SUM + TOLD(J) ♦ (B(K,1) * TOLO(J) *(B(K,2)/2.0 TOLD(J) * 28 

1 (B(K, 31/3.0 TQLO(J) * (8(K, 41/4.0 + TOLO( J ) * B(K, 5 1/5.0 > ) H - 29 

2 TK(N) * (B(K,1) 4- TK(N) ♦ (8(K, 21/2.0 + TK(N) ♦ (0(K,3)/3.O + TK 30 

3 (N) ♦ (8(K,4)/4.C + TK(N) ^ B ( K ,5 ) /5.0 ) )) ) 31 

U( J) = TNEW/COND(J) 32 

1 CONTINUE 33 

RETURN 34 

6 TEG(l) = 0.0 35 

00 5 1=1,2 36 

II = I 4- 1 37 

5 TEG(II) = TK(in * (B(I ,11 + TK(II) ♦ (B(I ,21/2.0 4- TK( I U ♦ ( 38 

1B( I ,3)/ 3.0 4- TKdlJ * (B(I ,41/4.0 4- TK(II) * B( I ,5)/5.0)))) - 39 

2TK( I) ♦ (B(I ,1) + TK(U ♦ (B(I ,2)/2.0 + TK( I ) * ( B( I ,3)/3.0 + 40 

3 TK(I) 4^ (B(I ,4)/4.0 4. TK(I) * 8(1 ,5)/5.0)> )) 41 

IN = 2 42 

GO TO 4 43 

END 44 


SUBROUTINE C0N3( TOLD , U ,NTI M,C QNO ) 1 

2 

CONDUCTIVITY COPPER 3 

TREF3 = 120.0 FOR SAMPLE PROBLEM 4 

5 

DIMENSION TOLDd) ,U(1) ,C0ND(50) , B( 6 , 5 ) ,TK ( 6) ,T EG(6 ) 6 

DOUBLE PRECISION TOLD , U,COND ,B , TK, TEG , SUM ,TNEW 7 

DATA (8( 1,J) , J=l,5)/1. 269232D0, -.66746332 0-1, .14229340-2, -.1399347 8 

19D-4, .524714240-7/ 9 

DATA (8( 2, J) , J=l, 5) /. 44088491 00 ,-.12795 5330-1 ,.137 510660-3, -.60595 10 

19480-6,. 823575130-9/ 11 

DATA ( B( 3, J) , Jr=l,5)/.4 9401865D-1 ,-. 667958 510-3 ,. 37 7709000-5, -. 9341 12 

138690- 8, .849680070-11 / 13 

DATA (B(4, J) , J = l,5)/. 70182 2020-2 1002119 0-4 ,. 216036330-7, -.21689 14 
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1429D- 10, . 81662 7650-14/ 15 

DATA ( B( 5, J) ,J:=1, 5)/. 520293820-2 191 831800-5 554337990-9, . 14533 16 

17480-12,0.00/ 17 

DATA (B(6,J), J = 1,5)/ .460-2,0.00,0.00,0.03,0.00/ 18 

DATA (TK(J),J = l,5)/80.00, 120.00, 350. DO, 863. 00, 1660. 00/ 19 

DATA IN/1/ 20 

IF ( IN .EO. 1) GO TO 6 21 

4 DD 1 J=1,NTIM 22 

20 K= X 23 

DO 2 L = 1,5 24 

LI = L 4- 1 25 

2 IF(TOLO(J) .GT. TK( L ) ) K = LI 26 

8 CONO(J)= B{K,1) + TOLD(J) * (B{K,2) ^ TOLD(J) ♦ (B(K,3) + TOLO(J)* 27 

1(B(K,4) + T3L0<J) * B(K,5)))) 28 

N = K-1 29 

IF(TOLD<J) .LE. 120.00) N = 2 30 

SUM = 0.0 31 

03 7 I = 1,N 32 

7 SUM = SUM + TEGU ) 33 

TNEW = SUM + TOLO(J) * <8(K,1) ♦ TOLO(J) *(B(K,2)/2.0 + TOLDCJ) * 34 

1 (8(K,3)/3.0 * TOLD(J) ♦ (BIK,4)/4.0 ^ TOLO(J) * B (K, 5 )/5.0 )) > ) * 35 

2 TK(N) * (B(K,1) + TK(N) * (B(K,2)/2.0 TK<N) * (B<K, 31/3.0 + TK 36 

3 <N) ♦ (8(K,4)/4.0 4- TK(N) ♦ 8( K,5 ) /5.0) ) ) ) 37 

U( J ) = TMEW/CONDU) 38 

1 CONTINUE 39 

RETURN 40 

6 TEG(l) = 0.0 41 

TEG(2) = 0.0 42 

03 5 I = 2,4 43 

II = I + 1 44 

5 TEG(II) = TKtin * (B(II,1) TKUI) * (B(II,2)/2.0 TKU I ) ♦ ( 45 

lB(II,3)/3.0 TKUn ♦ (8(II,4)/4.0 TK(I I) ♦ BU 1 ,5 ) /5.0 M ) ) - 46 

2TK(n ♦ (Bni,l) + TK(I) * (B(II,2)/2.0 + TKU) * (BUI,3)/3.0 + 47 

3 TKU) * (Bdl, 41/4.0 ^ TKU) * B( 1 1 ,5 ) /5 .0 ) ) ) ) 48 

TEGI6) = 0.0 49 

IN = 2 50 

GO TO 4 51 

END 52 


SUBROUTINE RL I NE < A ,EL , U,F , X, N) 1 

2 

ROUTINE 70 SOLVE AX = F WHERE A IS TRI-OIAGQNAL WITH CONSTANT OFF 3 

DIAGONAL ELEMENTS AND A = LU WHERE L AND U ARE LOWER AND UPPER 4 

TRIANGULAR MATRICES 5 

6 

DIMENSION EL( 1) ,U<1) ,F (1 ) ,xm ,Y(15) 7 

DOUBLE PRECISION A,EL,U,F,X,Y 8 

Nl = N-1 9 

Y(2) = F(2)/EL(2) 10 

DO I I = 3,N 11 

1 Y(I) = CF(I) - A * YU-in/ELU) 12 

X(N ) = YIN) 13 

DO 2 J = 2,N1 14 

I == N + 1 - J 15 

2 X( I) = YU) - Um ♦ X(H*1) 16 

RETURN 17 

END 18 


47 



oooo oooo 


SUBROUTINE CL I NE (E «,E L »U tF » X , N) 1 

2 

ROUTINE TO SOLVE AX = F WHERE A IS TRI-DI AGONAL AND A = LU,L AND U 3 

ARE LOWER AND UPPER TRIANGULAR MATRICES 4 

5 

DIMENSION EMlDfELd) fU(l) ,F(l),X(l) ,Y(35) 6 

DOUBLE PRECISION E M,E L ,U»F , X, Y 7 

N1 = N-l 8 

Y( 2) = F(2)/EL(2) 9 

03 1 I = 3,N 1C 

1 Y( n = <F(I) + EM(I-l) ♦ Y(I-in/EL(I) U 

X(N) = YIN) 12 

DO 2 J = 2,N1 13 

I = N-J +1 14 

2 X( n = YU) - UU) * X(H-l) 15 

RETURN 16 

END 17 


subroutine OUTPUT! I PUNCH) 1 

COMMON /F0RMA/R1,R2,RA2 ,YFINT,FI NTHK, ANGTOP ,HCOOL,T COOL , HGAS, TGAS* 2 

1XL3NGM ,R2B,YFIN8,ANGBDT,SUMA0,SUMAI 3 

DOUBLE PRECISION Rl ,R2 iRA2 , YF I NT ,F I NTHK , ANGTOP, HCOOLf TCOOL t HGASi 4 

ITGAS, XL0NGM,R2B,YFINB ,ANGBOT , SUMAC , SUMAI ,TOL 5 

COMMON /FORA/ A (15,25) ,A1 (35 ,25) ,8(15 ,15) ,BZ(20 ,15) ,D( 15, 15 ) , 6 

1 D1(35,15),TBA(25) ,TAI(25) ,TTA{25) ,TBB(15) ,T8 I ( 15 ) ,T80( 15) , 7 

2TDK 15),TD0(15) 8 

DOUBLE PRECISION A , A1 , B ,BZ ,D ,D1 , T8 A ,T AI ,TT A,T 8B,TB I ,T 8D, 9 

1T01,T0D 10 

COMMON /FORC/ P ( 20 ,1 5) , PI ( 20 ,5 ) , C( 15 ,20 ) , Q1 (1 5 , 20 ) , ZU 15, 20 ) , 11 

IZR ( 20,15) ,TP8(20) , TP I (20) ,TBQ(20) ,TQI (20) ,TZ B(20 ) ,TZ BR( 20 ) 12 

DOUBLE PRECISION P , PI , Q, Q1 ,ZL ,ZR ,T FB ,TPI ,T8Q, TCI ,TZ B, TZ BR 13 

COMMON /FORCN/ NXP , NXPl , NTQ , NTQl ,NZ, NZl ,N ZX, NZX 1 14 

COMMQN/FOUTA/THETA(25) ,RA(15) ,RA1X(35) ,R Al (35 ) ,T HET D( 15 ), 15 

ITHETBtlS) 16 

DOUBLE PRECISION THETA ,RA ,RA1 X ,RA1 ,THETD,THET B 17 

COMMON /F0UTC/THETQ(20) , XP (20) ,RQ A ( 15 ) , RQ(5 ) ,THET Z ( 20 ) , 18 

1 XZ (20), QT (20 ,YP(20) ,YP1 (5) 19 

DOUBLE PRECISION THE TQ , XP ,RQA , RQ ,THETZ ,XZ ,QT , YP,YP1 20 

COMM3N/FORAN/NRA,NRAM,NTA,NTAP,NRA1 ,NRA1M,NRA1X,NRA1XM, 21 

1NRB,NRBM,NT3,NT8M,NT0 ,NTDM 22 

DATA CON V/57. 2957795/ 23 

DIMENSION U( 80) , TEMP (80) 24 

25 

IPUNCH 1 PUNCH BCD CAROS WITH U DISTRIBUTION FOR ALL 26 

SUBREGIONS AND TEMPS ON BOUNDARIES AND INTERFACES 27 

28 

IF( IPUNCH .NE. 1) GO TO 600 29 

CALL BCOUMP < A ( 1 , 1 ) ,A ( 1 5 ,25 ) ,0) 30 

CALL BCDUMP ( A 1 ( 1 , 1) , Al (3 5 ,25 ) ,0 ) 31 

CALL BCDUMP ( D (1 , 1 ) ,0 ( 15 ,15 ) ,0) 32 

CALL BCDUMP ( D I (1 , 1) ,D1 (3 5 ,15 ) ,0 ) 33 

CALL BCOUMP ( B ( 1 , 1 ) ,B ( 15 , 15 ) ,0) 34 

CALL BCOUMP ( B Z( 1 , 1) ,8 Z(20 ,1 5 ) ,0 ) 35 

CALL BCOUMP ( T8 A { 1 ) ,TB A ( 2 5 ) ,0) 36 

CALL BCDUMP ( TA I( 1) , TA I { 2 5 ) ,0) 37 

call BCDUMP (TTA(l) ,TTA(25) ,0) 38 

CALL BCOUMP ( TBD ( 1 ) ,TBD ( 15 ) , 0 ) 39 

CALL BCOUMP ( TOI { 1 ) ,TD I ( 15 ) ,0 ) 40 

CALL BCDUMP ( TOD ( 1 ) ,TDD( 1 5 ) ,0 ) 41 

CALL BCDUMP ( TBB ( 1 ) , TBB (1 5 ) , 0) 42 

CALL BCDUMP ( T8I ( 1 ) ,TB I ( 1 5 ) ,0 ) 43 

CALL BCDUMP ( P ( 1 ,1 ) ,P ( 20 ,1 5) ,0 ) 44 

CALL BCOUMP ( P 1 ( 1 , 1) , PI (20 ,5 ) ,0 ) 45 

CALL BCDUMP ( Q (1 , 1 ) ,Q ( 15 ,20) ,0 ) 46 

CALL BCDUMP ( Q 1 ( 1 , 1 ) , Q1 (1 5 ,20) ,0 ) 47 

CALL BCDUMP ( ZL ( 1 , 1) , ZL( 1 5 ,20 ) ,0 ) 48 

CALL BCDUMP ( ZR ( 1 , 1 ) , ZR (20 ,l 5 ) ,0 ) 49 

CALL BCOUMP ( TPB ( 1 ) , TPB( 2 0) ,0) 50 

CALL BCOUMP ( TPI { 1 ) , TP I ( 20) ,0) 51 

CALL BCOUMP (TBQ(l) ,TBQ(20) ,0) 52 

CALL BCDUMP ( TQI ( 1 ) ,TQI (20) ,0) 53 

CALL BCOUMP ( TZB ( 1 ) , TZB ( 2 0) ,0 ) 54 

CALL BCDUMP ( TZ8R ( 1 ) , T2BR (20) ,0 ) 55 

6CC CONTINUE 56 

C 57 


48 



ooo non ooo ooo ono ooo ooo 


C REGION A 58 

C 59 

T1 = TBA(l) 60 

WRITE (6, 100) 61 

DD 1 I=2,NTAP 62 

THETAII) = THETAd) CONV 63 

WRITE! 6, 101) THETAU-1) 64 

DD 2 J=1,NRA 65 

2 U(J) = AUtl) 66 

CALL SLBT1(U,T1,TEMP,1 ,NRA) 67 

T1 = TEMPIl) 68 

T2 = TEMP(NRA) 69 

WR ITEI 6» 102) (RAU) ,TEMP( J) ,J=1 ,NRA) 70 

71 

REGION A1 72 

73 

IF(R2 .EQ. RA2) GO TO 1 74 

03 4 J=1,NRA1X 75 

4 IJ( J) = AKJfl) *^6 

CALL $UBT2(U,T2,TEMP,I ,NRA1X) 77 

WR ITEI 6, 102) (RAl ( J) ,TEMP( J) ,J==1 ,NRA1X) 78 

1 CONTINUE 79 

80 

REGION D 81 

82 

00 14 I = 2,NTD 83 

THB = THETOm ♦ CONV + THETA(NTA) 84 

WRITE (6,101) THB 85 

DO 15 J = 1,NRA 86 

15 U(J) = D(J,I) 87 

CALL SUBT1(U,T1,TEMP,1,NRA) 88 

T1 = TEMP(l) 89 

T2 = TEMP(NRA) 90 

WRITE(6,1C2) (RA (J) , TEMPI J) , J= 1,NRA) 91 

IF(R2 .EQ, RA2) GO TO 14 92 

93 

REGION 01 94 

95 

JJ = NRAIX +1-1 96 

03 16 J = 1,JJ 97 

16 U( J ) = OK J, I ) 98 

CALL SUBT2(U,T2,TEMP,1 ,JJ) 99 

WRITE! 6,102) (RAl U) , TEMP ( J) ,J=1 ,JJ) 100 

lA CONTINUE 101 

102 

REGION B 103 

104 

THBL = THB 105 

03 5 I = 1,NTB 106 

THB = THETBU) ♦ CONV + THBL 107 

WRITE (6,101) THB 108 

03 6 J = 1,NRA 109 

6 U(J ) = B(J,I ) 110 

CALL SUBTKU,T1,TEMP,1,NRA) 111 

T1 = TEMP(l) 112 

T2 = TEMP(NRA) 113 

WRITE(6,102) (RA( J) , TEMPI J) , J=1,NRA) 114 

115 

REGION B1 116 

117 

JJ = NRAl +1-1 118 

03 7 J = 1,JJ 119 

7 U(J) = BZ(J,n 120 

CALL SUBT3(U,T2,TEMP,1 ,JJ) 121 

MRITE(6,102) IRAK J),TEMPU) , J= 1 ,JJ) 122 

5 CONTINUE 123 

124 

REGION P 125 

126 

WRITE(6,120) 127 

03 30 1=1, NXP 128 

WRITE(6,IG5) XPU) 129 

03 31 J=1,NRA 130 

31 U( J )= P( !,J) 131 

CALL SUBTUU,T1,TEMP,1 ,NRA) 132 

WRITE(6,102) (YP(J), TEMP(J), J=1,NRA) 133 

134 

REGION PI 135 

136 
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T1 = TEMP(l) 137 

T2 = TEMP(NRA) 138 

DG 32 J=1,NR8 139 

32 U( J )= P1(I,J1 lAO 

CALL SL8T3(U,T2,TEHP,1 ,NRBI lAl 

WRITE(6,102» (YPKJI, TEMPCJ), J=1,NRB) 142 

3P CONTINUE 143 

C 144 

C REGION Q 145 

C 146 

WRITE (6,122) 147 

DO 33 I:=1,NT0 148 

THB = THETQd)* CONV 149 

WRITE( 6, 101) THB 150 

00 34 J=1,NRA 151 

34 U( J )= Q( J,I) 152 

CALL SUBT1(U,T1,TEMP,1 ,NRA) 153 

WR I TE( 6, 1(>2) (RQA(J), TEMP(J), J=1,NRA) 154 

155 

REGION Q1 156 

157 

T1 = TEMP(l) 158 

T2 = TEMP(NRA) 159 

03 35 J=1,NR3 160 

35 U(J )= Ql(J,n 161 

CALL SUB T3( U,T2, TEMP, 1 ,NRB) 162 

WRITE(6,102) (RQ{J), TEMPI J), J»1 , NRBM) ,QT ( I) ,T EMP( NRB » 163 

33 CONTINUE 164 

165 

REGION ZL 166 

167 

THBL = THB 168 

DO 38 1=2, NZ 169 

THB = THETZtl) *CONV + THBL 170 

WRITE( 6,101) THB 171 

00 39 J=l, NRA 172 

39 U(J)= ZL{J,I) 173 

CALL SL8T1(U,T1, TEMP, 1, NRA) 174 

T1 = TEMPI 1) 175 

WRITE(6,102) (RQA(J), TEMPU), J=1,NRA) 176 

38 CONTINUE 177 

178 

REGION ZR 179 

180 

WRITE(6,123) 181 

DO 40 I = 1,NZX 182 

WRITE(6,105) XZ(I) 183 

DO 41 J=1,NRA 184 

41 U(J)= ZR(I,J) 185 

CALL SUBT1(U,T1,TEMP,1 ,NRA) 186 

T1 = TEMP(l) 187 

WRITE(6,102) (YP(4), TEMPU), J=1,NRA) 188 

40 CONTINUE 189 

ICO F0RMAT(1H1,37X46HTEMPERATURES ON THE TOP RADIAL SECTION OF TU8E//1 190 

IH ) 191 

101 FORMAT(1HK,50X7HTHET = El 3# 6//9X1HR ,15 XIHT , 15 XIHR, 15X IHT, 192 

115X1HR ,15X1HT,15X1HR,15X1HT/1H ) 193 

102 F0RMAT(8( 3XE13.6) ) 194 

105 F3RMAT(1HK,53X4HX = El 3. 6//9XIH Y ,1 5 XIHT , 1 5XIH Y , 15X IHT , 195 

U5X1HY,15X1HT, 15X1HY,15X1HT/1H ) 196 

120 F0RMAT(1HK,37X46HTEMPERATURES ON THE MIDDLE RECTANGULAR SECTION//! 197 

IH ) 198 

122 F0RMAT(1HK,37X48HTEMPERATURES ON THE LOWER RADIAL SECTION OF TUBE/ 199 

1/lH ) 200 

123 F0RMATaHK,37X45HTEMPeRATUR6S ON THE LOWER RECTANGULAR SECTI0N//1H 201 

1 ) 202 

RETURN 203 

END 204 


SUBROUTINE SIJBTUU,T1 ,TEMP,N1 ,N2) 1 

2 

SUBROUTINE TO TRANSFORM U FOR STAINLESS TO TEMPERATURE 3 

TREFl = 70.0 FOR SAMPLE PROBLEM 4 

5 

DATA ( 8( 1, J) , J =1,5)/. 803 80307D-4,. 63751499 D-6,-. 1420 1613D-8, .1083 6 
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171240-11,0.00/ 7 

DATA J=1,5)/.11520489D-3».194799610-6»-.1240770-9, .624877 8 

UD-13f 0.00/ 9 

DATA (B(3,J), J =1, 5)/. 40il3088D-3f-. 384859870-6, .33912452D-9,-. 99 10 

12039140- 13 t. 106384240-16/ 11 

DATA (TK(J),J = 1,3»/70.00,400.DO,12 00 . 00/ 12 

DIMENSION U(80) , TEMP {80) ,B <6 ,5 ) t TK (5 ) ,TEG(6) 13 

DOUBLE PRECISION B ,TK 14 

DATA IN/1/ 15 

IF ( IN .EQ. 1) GO TO 6 16 

4 DO 1 I=N1,N2 17 

3 < = 1 18 

DO 2 L = 1,2 19 

L 1= L + 1 20 

2 IF(T1 .GT. TKI LU ) K = LI 21 

C3ND = B(K,1) + T1 ♦ (B(K,2) + T1 ♦ (B(K,3) + T1 ♦ 22 

1(B(K,4) + T1 * B{K,5) )} ) 23 

N == K 24 

SUM = 0.0 25 

D3 7 M = 1,N 26 

7 SUM = SUM + TEG(M) 27 

TNEW = SUM 4- T1 ★ (B(K,n + T1 *(8(K,2)/2.0 + TI ♦ 28 

1 (B(K,3)/3.0 + Tl * (B(K,4)/4.0 + Tl ♦ B ( K, 5 ) /5.0 ) H ) - 29 

2 TK(N) ♦ (B(K,1) 4- TKINI * (B(K, 21/2.0 4- TK(N) ♦ (0(K, 31/3.0 4- TK 30 

3 (N) * (B(K, 41/4.0 4- TK(N) * B(K,5)/5.0nM 31 

TNEW = Tl - {-U(n ^ TNEW)/COND 32 

REL = ABS< (Tl-TNEWI/TNEW) 33 

T1=TNEW 34 

IF (REL.LE.l.OE-6) GO TO 1 35 

GO TO 3 36 

1 TEMP( I ) = TNEW 37 

RETURN 38 

6 TEGCl) = 0.0 39 

DO 5 I = 1,2 40 

II I + 1 41 

5 TEG(in = TXCin * (BU ,U + TK(in * (811 ,2)/2.0 + TK( I 1) * C 42 

1B( I ,3)/3.0 4 - TK(II) * (B(I ,41/4.0 4 - TKlII) ♦ B(I ,5)/5.Dn)) - 43 

2T<(I) * (B(I ,n + TK(U ♦ (B(I ,2)/2.0 + TK( I) * ( B( I ,31/3.0 4 - 44 

3 TK(I ) ♦ (B( I ,41 /4.0 4- TK(I) ♦ B(I ,5)/5.0),) )) 45 

IN 2 46 

G3 TO 4 47 

END 48 


SUBROUTINE SUBT2(U,T1 ,TEMP,N1,N2) 1 

2 

SUBROUTINE TO TRANSFORM U FOR NICKEL TO TEMPERATURE 3 

TREF2 = 230.0 FOR SAMPLE PROBLEM 4 

5 

DATA { 8( 1, J), J=l,5> / .40450270-2 ,-.134431050-4, .190140140-7, -.8813 6 

14126D-ll,-.63 8 609670-15/ 7 

0ATA(B(2,J), J = l,5)/ .16581 2 570-2, -.272994110-5,. 344585970-8, -.-238 8 

1504690-11, ,678439540-15/ 9 

DATA (B( 3, J) ,J=1,5)/. 432941010-3 ,. 222904460-6 , -.17 147530-9 ,. 120048 10 

1780- 12, -.23650 193-16/ 11 

DATA ( TK{ J ), J = l,3) /230.00, 660.00, 1360.00/ 12 

dimension U( 80) , TEMP (80) ,B (6 ,5 ) , TK (5 ) ,TEG(6) 13 

DOUBLE PRECISION B,TK 14 

DATA IN/1/ 15 

IF ( IN .EQ. 1) GO TO 6 16 

4 DO 1 I=N1,N2 17 

3 K = 1 18 

DO 2 L=l,2 19 

Ll= L+1 20 

2 IF(T1 .GT. TK(L1M K * LI 21 

CONO = B(K,1) + Tl ♦ (B(K,2) + Tl ♦ (B(K,3) + Tl ‘ ♦ 22 

1(B(K, 4) + Tl ♦ 8(K,5) )) ) 23 

N ^ K 24 

SUM =5 r.O 25 

DO 7 M = 1,N 26 

7 SUM =: SUM + TEG(M) 27 

TNEW = SUM + Tl * (81K,1) ^ Tl ♦(B(K,2)/2.0 + Tl ♦ 28 

1 (B(K,3)/3.0 ^ Tl ♦ (B(K, 4)74.0 ♦ Tl ♦ B (K, 5 )/5.0 ) > )) - 29 

2 TK(N) ♦ (8(K,1) TK(N) ♦ (8(K, 21/2.0 TK(N) ♦ <B(K,3)/3.0 ♦ TK 30 

3 (N) * (B(K,4)/4.0 + TK(N) ♦ B(K,5)/5.0n)) 31 

TNEW ^ Tl - (-U(I) ♦ TNEW)/COND 32 
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noon 


REL = ABS({T1*-TN6W)/TNEWI 33 

T1=TNEV< 3^ 

IF (REL.LE.l.OE-6) GO TO I 35 

GO TO 3 36 

I TEMP( r ) = TNEW 37 

DEBUG (U(I I = Nl,M2» 38 

RETURN 39 

6 TEG( 1) =0.0 40 

DO 5 1=1,2 41 

11=1+1 42 

5 TEGdl I = TKCin ♦ iBil ,1) + TKtlH ♦ (B(I ,2)/2.0 + TK( 1 1) ♦ ( 43 

18(1 ,3)/3.0 + TK(in ♦ (B(I , 41/4.0 + TK(I I) ♦ BU ,51/5.0)))) - 44 

2TK(ri ♦ (B(I ,1) + TK(U ♦ (6(1 ,21/2.0 + TK( H ♦ (6(1 ,3)/3.0 + 45 

3 TKU) * (B(I ,4i/4.0 + TK(I ) ♦ B« I ,5)/5.0)))) 46 

IN = 2 47 

GO TO ^ 48 

END 49 


SUBROUTINE SUB T3 ( U,Tl , TE MP ,N1 ,N2) 1 

2 

SUBROUTINE TO TRANSFORM U FOR COPPER TO TEMPERATURE 3 

TREF3 = 120.0’ FOR SAMPLE PROBLEM 4 

5 

DATA ( B( 1, J) , J = l,5)/l. 26923200, -.66746332 0-1, .14229 340-2, -.1399347 6 

190-4, .52^714240-7/ 7 

DATA (S( 2, J),J = 1, 5)/. 44088491 DO, -.127955330-1 137 51066D-3, -.60595 8 

19480- 6,. 8235 75130-9/ 9 

DATA CB(3,J),J = 1, 5)/. 494018650-1 ,-.6679585 lD-3 , .37 7709000- 5, -. 9341 10 

13869D- 8, .849680070-11 / 11 

DATA (B(4,J),J=1,5)/.701 822020-2 ,-.10021190-^,. 2160 3633D-7,-. 21689 12 

14290-10, .81662765D-14/ 13 

DATA ( 8(5, J) ,J=l,5)/.52029382D-2.,.19l83180D-6 ,-.554337990-9, .14533 lA 

17480-12,0.00/ 15 

DATA (B(6,J), J = 1,5)/ .460-2,0.00,0.00,0.03,0.00/ 16 

DATA (TK(J),J = l,5)/80.00, 120.00, 360. DO, 860. 00,1660.00/ 17 

DIMENS ION U( 80) , TEMPI 80) ,B (6 ,5 ) , TK ( 5 ) ,TEG(6) 18 

DDUBLE PRECISION B,TK 19 

DATA IN/1/ 20 

IF ( IN .EQ. 1) GO TO 6 21 

4 00 1 I=N1,N2 22 

3 K = 1 23 

03 2 L = 1,5 24 

Ll= L+1 25 

2 IF (Tl.G7.TKiL )) K=L1 26 

C3ND = B(K,1) + T1 ^ (B(K,2) + T1 ♦ (B(K,3) + T1 ♦ 27 

1(B(K,4) + T1 * 8(K,5)>) ) 28 

N = K-1 29 

SUM = 0.0 ’ 30 

IF( T1 .LE. 120.00) N = 2 31 

00 7 M = 1,N 32 

7 SUM = SUM + TEGiM) 33 

TNEW = SUM + T1 ♦ (8(K,1) + T1 ♦ ( B( K, 2 )/ 2 .0 + T1 * 34 

1 (BCK,3)/3.0 + T1 * (B(K,4)/4.0 + T1 ♦ B (K, 5 )/5.0 ) ) ) ) - 35 

2 T<(N) * (B(K,1) + TK(N) ♦ (B(K,2)/2.C t TK(N) * {B(K,3)/3.0 + TK 36 

3 (N) * (8(K,4)/4.0 + TK(N) ♦ B ( K ,5 ) /5 .0 ) ) ) ) 37 

TNEW = T1 - (-Um + TNEW)/COND 38 

REL = ABS( (Tl- TNEW) /TNEW) 39 

T1 = TNE)^ 40 

IF iREL.LE.l.OE-6) GO TO 1 41 

G3 T3 3 42 

1 TEMPI I) = TNEW 43 

DEBUG (Um, I = N1,N2) 44 

RETURN 45 

6 TEG( 1 ) = 0.0 46 

TEG(2) = 0.0 47 

03 5 I = 2,4 48 

II = I + 1 49 

5 TEGdl ) = TK(II) * (B(II ,1) + TKIIIJ ♦ <B(II,2)/2.0 + TKt I I) * ( 50 

1B( II,3)/3.0 + TKiin * (Bdl ,4)/4.0 + TK(II) * B( I 1 , 5 ) /5.0 )) ) ) - 51 

2T<(I) ♦ (Bdl,l) + TK(I) ♦ (B(II,2)/2.0 + TK( I) ♦ (B(II,3)/3.0 + 52 

3 TK{ n * (B( II ,4) /4.0 + TKU ) ♦ Bdl ,5)/5.0)) ) ) ^ 53 

TEG(6) 0.0 54 

IN = 2 55 

G3 TO 4 56 

END 57 
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